Khoa học tính toán

Cách hiệu quả nhất để có được một ma trận thưa thớt phức tạp từ mã Fortran của tôi sang PETSc là gì? Tôi hiểu rằng đây là vấn đề phụ thuộc, vì vậy tôi đã cố gắng cung cấp càng nhiều chi tiết có liên quan càng tốt dưới …

8 matrix  petsc  performance  fortran  sparse-matrix 

Tôi đang tìm kiếm các ví dụ về phương trình Helmholtz và Biharmonic trong tọa độ Cartesian với các giải pháp chính xác, để so sánh các giải pháp số của tôi với nó. Tôi đã có thể tìm thấy khá nhiều ví dụ trên internet, nơi vấn đề với …

8 pde  finite-element  reference-request  elliptic-pde 

Tôi có hai biến kvà tlà hàm của hai biến khác p1và p2. Tôi cũng biết giá trị tối đa của họ. Tôi không có bất kỳ biểu hiện phân tích cho điều này. Tôi muốn tìm các giá trị của kvà tđó là gần nhất với giá trị tối …

8 optimization 

Tôi đang cố gắng tìm hiểu thêm về chia lưới âm lượng và đã quyết định thử thực hiện một mesher âm lượng đơn giản. Chiến lược tôi đã chọn là chia nhỏ không gian của mình bằng cách sử dụng một quãng tám, được tinh chỉnh dựa trên một …

8 finite-element  mesh-generation 

Có một triển khai mã nguồn mở, an toàn luồng của các bộ giải tối ưu hóa lồi trong C / C ++ không? Một số thư viện như NLopt , Ipopt , OPT ++ không đáp ứng yêu cầu của tôi. OPT ++ và Ipopt không an toàn cho …

8 optimization  c 

Tôi cần phải giải hệ hai phương trình vi phân từng phần được ghép bằng số. ∂x1∂t∂x2∂t= c1∇2x1+ f1( x1, x2)= c2∇2x2+ K∂x1∂t∂x1∂t=c1∇2x1+f1(x1,x2)∂x2∂t=c2∇2x2+K∂x1∂t\begin{align} \frac{\partial x_1}{\partial t} &= c_1\nabla ^2 x_1 + f_1(x_1,x_2)\\ \frac{\partial x_2}{\partial t} &= c_2\nabla ^2 x_2 + K\frac{\partial x_1}{\partial t} \end{align} Miền của hệ thống là một …

8 pde  fourier-analysis 

Giả sử rằng một mô phỏng số trực tiếp được thực hiện, một phương pháp tốt để ước tính tác động của quy mô nhỏ trên quy mô lớn trong động lực học chất lỏng là gì? Ví dụ, có thích hợp để so sánh hai lần chạy với kích …

8 fluid-dynamics  statistics  simulation 

Tôi đã đọc về việc thực hiện phép nhân ma trận dày đặc khi ma trận không phù hợp với bộ đệm. Một trong những biểu đồ tôi đã thấy (slide 9 từ các slide này ) cho thấy hiệu suất giảm đột ngột khi sử dụng thuật toán ngây …

8 linear-algebra  matrix  blas 

Tôi muốn làm cho chương trình python của mình nhanh bằng cách sử dụng cython, nhưng vòng lặp bên trong của tôi vẫn thực hiện các cuộc gọi python chậm đến trình tạo số ngẫu nhiên! Vài năm trước, vấn đề tương tự đã được đưa ra bởi một người …

8 python  random-number-generation  sage 

Gần đây tôi đã thừa hưởng một khối lớn mã di sản giải quyết một vấn đề rất khó khăn, thoáng qua. Tôi muốn chứng minh rằng kích thước bước không gian và thời gian đủ nhỏ để bản chất định tính của giải pháp tính toán sẽ không thay …

8 convergence  accuracy  time-integration 

Tôi có một loạt các quy trình MPI trao đổi thông điệp yêu cầu qua lại. Các quy trình không biết quy trình nào khác sẽ gửi tin nhắn cho họ hoặc bao nhiêu. Với tình huống này, tôi muốn một cách hiệu quả để biết liệu tất cả các …

8 mpi 

Phương trình nhiệt được rời rạc trong không gian với FV (hoặc FEM) và thu được phương trình bán rời (hệ thống ODE). Cách tiếp cận này, được gọi là phương pháp của các dòng , cho phép dễ dàng chuyển đổi từ sự rời rạc tạm thời này sang …

8 finite-element  pde  finite-volume  software  rothe-method 

Tôi đã thấy rất nhiều ấn phẩm trong các tạp chí Vật lý tính toán sử dụng các số liệu khác nhau để thực hiện mã của họ. Đặc biệt đối với mã GPGPU, dường như có rất nhiều kết quả về thời gian mà mọi người công bố. Cụ …

8 performance  publications  documentation 

Tôi xin lỗi trước nếu câu hỏi này là ngớ ngẩn. Tôi cần tính toán gốc của bạn - f( u ) = 0u−f(u)=0\begin{equation} u -f(u) =0 \end{equation} Trong đó là một vectơ thực và là một hàm có giá trị vectơ thực. Tôi đã bắt đầu với phương pháp …

8 iterative-method  nonlinear-equations  roots 

Tôi có một phần của một vấn đề được mô tả bởi phương trình bảo toàn động lượng: ∂ρ∂t+1sinθ∂∂θ(ρusinθ)=0∂ρ∂t+1sin⁡θ∂∂θ(ρusin⁡θ)=0\frac{\partial \rho}{\partial t} + \frac{1}{\sin\theta} \frac{\partial}{\partial \theta}(\rho u \sin \theta) =0 Trong đó và (vận tốc không đổi).u=f(θ)u=f(θ)u=f(\theta)ρ=f(θ,t)ρ=f(θ,t) \rho = f(\theta,t) Chắc chắn người ta có thể áp dụng một trong những …

8 finite-difference  fluid-dynamics  hyperbolic-pde