Khoa học tính toán

Tôi nhận được cảnh báo trong chủ đề bài viết khi cố gắng tối ưu hóa một chức năng trong Python bằng hàm scipy.optizes.fmin_bfss . Đầu ra hoàn chỉnh: Cảnh báo: Lỗi mong muốn không nhất thiết phải đạt được do mất độ chính xác Current function value: nan Iterations: …

8 optimization  python 

Tôi chỉ cần tính toán một số lượng nhỏ các thành phần Fourier tần số thấp của mảng 2 chiều phức tạp. Tôi sẽ tính toán các thành phần Fourier tương tự lặp đi lặp lại khi mảng đầu vào thay đổi. Rõ ràng, trong giới hạn mà tôi chỉ …

8 matlab  fourier-analysis 

Tôi đang cố gắng làm MD trên hơi nước. Theo tôi biết có tồn tại một số mô hình nước cho nước lỏng, chẳng hạn như SPC, SPC / E, TIP3P, nhưng chúng cũng sẽ áp dụng cho trạng thái hơi của nước? Và sự khác biệt của mô phỏng …

8 molecular-dynamics 

Theo như tôi hiểu, thứ tự mà các tin nhắn MPI điểm-điểm không chặn (Isend và Irecv) được nhận phù hợp với thứ tự chúng được gửi. Có bất kỳ kỹ thuật để ưu tiên tin nhắn nhất định hơn những người khác? Ví dụ: chúng tôi có một thuật …

8 mpi 

Tôi đã triển khai CG trong FORTRAN bằng cách liên kết nó với Intel MKL. Khi có các câu như: ( Tham khảo Wikipedia ) p=r; x=x+alpha*p r=r-alpha*Ap; hoặc những cái tương tự trong QMR (với số lượng lớn hơn nhiều) v_tld = r; y = v_tld; rho = norm( …

8 performance  fortran  blas 

Tôi đã cố chạy các trường hợp thử nghiệm cho CG và tôi cần tạo: ma trận xác định dương đối xứng có kích thước> 10.000 ĐẦY ĐỦ Chỉ sử dụng các chỉ số ma trận và nếu cần 1 vectơ (Giống như )A(i,j)=x(i)−x(j)(i+j)A(i,j)=x(i)−x(j)(i+j)A(i,j) = \dfrac{x(i) - x(j)}{(i+j)} Với số …

8 linear-algebra 

Một số tài nguyên tốt (sách, bài viết, trang web) về giao diện đa giác và thuật toán hợp là gì?

8 computational-geometry  reference-request 

Tôi viết mã bằng Fortran và C cho các thuật toán ma trận khác nhau. Tuy nhiên, khi tôi cấu hình mã của mình bằng VTune, tôi thường chạy theo một số thuật ngữ mà tôi không thể đánh giá đầy đủ. Có một nguồn tài nguyên tốt để học …

8 parallel-computing  reference-request 

Giả sử bạn giải quyết một hệ thống tuyến tính bằng một phương pháp lặp, ví dụ: độ dốc liên hợp hoặc phép lặp Richardson. Sau đó, bạn cố gắng giải quyết một hệ thống tuyến tính hơi bị nhiễu trong ma trận và phía bên tay phải, giả sử, …

8 linear-algebra  finite-element  linear-solver 

Tôi đang tìm kiếm một thư viện C ++ miễn phí, thuận tiện cho các biểu đồ bao gồm một bộ giải cho bài toán đồ thị con chung (MCS) tối đa. Tôi biết về Thư viện đồ thị Boost và LEMON , nhưng không bao gồm bộ giải MCS.

8 libraries  c++  graph-theory 

Cách tốt, phổ biến để diễn đạt từ điển (= maps) trong mã giả là gì? Các cơ sở dữ liệu về cơ bản cho phép lưu trữ các giá trị cho các khóa, lặp lại trên tất cả các cặp khóa / giá trị, kiểm tra để bao gồm …

8 algorithms  publications 

Tôi đang chơi với các triển khai đơn giản của các thuật toán để tìm độ phủ đỉnh tối thiểu / khớp thẻ tối đa phù hợp trong các biểu đồ lưỡng cực. Tuy nhiên, tôi dường như gặp khó khăn khi googling cho một số bộ dữ liệu thử …

8 data-sets  graph-theory 

Vấn đề này phát sinh từ một dự án mô hình thống kê Bayes. Để tính toán với mô hình của tôi, tôi cần thực hiện tích hợp trong đó một phần của tích phân là Phân phối "Pólya" hoặc "Dirichlet-Multinomial", p ( n ∣ α ) = ( N! …

8 quadrature 

"Ưu và nhược điểm" của hai phương pháp sàng lọc lưới này là gì? Cả hai dường như là phương pháp phổ biến. Tôi có thể tự nhiên tưởng tượng rằng tinh chế màu đỏ toàn cầu tương đối dễ thực hiện và hấp dẫn cho nghiên cứu trong các …

8 hpc  adaptive-mesh-refinement 

Tôi đã tạo một spline để phù hợp với dữ liệu của mình trong python bằng cách sử dụng: spline=scipy.interpolate.UnivariateSpline(energy, fpp, k=4) Phương trình tôi muốn sử dụng bao gồm tổng giữa n = 2 và n = vô cùng, trong đó n là thứ tự của vi phân tại …

8 python  interpolation