Khoa học tính toán

Tôi vừa đọc Chương 3 trong "Hướng dẫn đa điểm" của Briggs / Henson / McCormick, liên kết . Văn bản nói về các chu trình Multigrid như chu kỳ V, chu trình mu, FMG. Điều khiến tôi chú ý: Trong hầu hết các quy trình lặp, người ta sẽ …

12 multigrid 

Đây có thể là một câu hỏi phù hợp hơn cho phía Khuyến nghị phần mềm của SE, tuy nhiên tôi tin rằng những người thường xuyên sử dụng phần này của SE có nhiều khả năng có thể trả lời câu hỏi này. Tôi đang tìm kiếm một sự …

12 finite-element  comsol  octave  openmodelica  modeling 

Giả sử tôi có các phiên bản chậm và nhanh của một số mã và muốn báo cáo số tăng tốc so sánh hai mã. Tôi chạy phiên bản chậm lần và phiên bản nhanh lần, tạo ra thời gian và . Cách đơn giản nhất để tạo tốc độ …

12 statistics  publications 

Chúng tôi đang nghiên cứu mô hình Bayes cho quy trình không gian thời gian và đang sử dụng bộ lấy mẫu No-U-Turn (NUTS) yêu cầu một mô hình cho xác suất log và độ dốc của nó đối với các tham số mô hình. Ngắn gọn hơn, chúng ta …

12 monte-carlo  automatic-differentiation 

Hãy và f ( → x ) : [ 0 , 1 ] n → C là một chức năng trong các biến này.x⃗ =(x1,x2,…,xn)∈[0,1]nx→=(x1,x2,…,xn)∈[0,1]n\vec{x} = (x_1, x_2, \dots, x_n) \in [0,1]^nf(x⃗ ):[0,1]n→Cf(x→):[0,1]n→Cf(\vec{x}): [0,1]^n \to \mathbb{C} Có một sơ đồ đệ quy cho tích phân lặp này không? ∫[0,1]n∏dxif(x⃗ )∫[0,1]n∏dxif(x→)\int_{[0,1]^n} …

12 numerical-analysis  fourier-analysis 

Trong các phương trình Navier-Stokes không thể nén, thuật ngữ áp lực thường được gọi là hệ số nhân Lagrange thi hành điều kiện không thể nén.ρ ( bạnt+ ( u ⋅ ∇ ) u )∇ ⋅ u= - ∇ p + μ Δ u + f=0ρ(ut+(u⋅∇)u)=−∇p+μΔu+f∇⋅u=0\begin{align*} \rho\left(\mathbf{u}_t + (\mathbf{u} …

12 optimization  fluid-dynamics  navier-stokes  incompressible 

Tôi không quen với định dạng tệp được sử dụng trong HDF5, nhưng tôi tự hỏi liệu các tệp HDF5 có phù hợp để kiểm soát sửa đổi với git (hoặc ví dụ Mercurial hoặc Subversion) không? Tôi đoán điều tôi muốn nói là: các tệp HDF5 phù hợp với …

12 data-management  data-storage 

Tôi đang cố gắng tìm hiểu xem có cách nào nhanh hơn để tính toán tất cả các giá trị riêng và hàm riêng của một ma trận kề kề rất lớn và thưa thớt hơn là sử dụng scipy.spude.linalg.eigsh Theo tôi biết, phương pháp này chỉ sử dụng độ …

12 linear-algebra  python  performance  eigensystem  scipy 

Hệ thống của tôi là một vấn đề FE đối xứng với các số nhân trễ (ví dụ như luồng Stokes không thể nén): ( MộtBBTC)(MộtBTBC)\begin{pmatrix}A & B^T \\ B & C\end{pmatrix} Trong đó là trường hợp điển hình (tôi thậm chí đã chắc chắn rằng các phương trình được …

12 finite-element  petsc  krylov-method  preconditioning 

Khi FEM-discretizing và giải quyết một vấn đề phản ứng-khuếch tán, ví dụ: với 0 < ε « 1 (số ít nhiễu loạn), các giải pháp của vấn đề rời rạc thường sẽ triển lãm các lớp dao động gần với ranh giới. Với Ω = ( 0 , 1 …

12 finite-element  discretization  numerical-analysis  singular-perturbation 

Này và cảm ơn vì đã dành thời gian để xem xét câu hỏi của tôi. Đây là phiên bản cập nhật của câu hỏi của tôi mà tôi đã đăng trước đó trong vật lý.stackexchange.com Tôi hiện đang nghiên cứu một máy kéo sợi exciton 2D Bose-Einstein Condensate và …

12 matlab  ode  quantum-mechanics 

Một số họ phương thức tồn tại cho phương trình số. Nếu tôi có một lớp tích phân cụ thể, làm thế nào để tôi chọn phương thức lý tưởng? Các câu hỏi liên quan để hỏi cả về tích phân (ví dụ: nó có trơn tru không? Nó có …

12 quadrature 

Giả sử tôi có một vấn đề giá trị biên: d2bạndx2+ dvdx= f trong Ωd2bạndx2+dvdx= =f trong Ω\frac{d^2u}{dx^2} + \frac{dv}{dx}=f \text{ in } \Omega dbạndx+ d2vdx2= g trong Ωdbạndx+d2vdx2= =g trong Ω\frac{du}{dx} +\frac{d^2v}{dx^2} =g \text{ in } \Omega u = h trong ∂Ωbạn= =h trong ∂Ωu=h \text{ in } \partial\Omega …

12 pde  iterative-method 

Tôi đang chuyển một mã hiện có từ MATLAB sang C ++ và có một hệ thống tuyến tính để giải (chứ không phải là dạng A x = b điển hình hơnxA=bxA=bxA=bAx=bAx=bAx=b ) Ma trận dày đặc và có dạng chung, nhưng không lớn hơn 1000x1000. Vì vậy, trong …

12 linear-solver  c++  lapack 

Tôi đang mô phỏng dòng chảy không thể nén trên một hình trụ ở số Reynold là 500. Tôi đang giải phương trình stokes điều hướng bằng phương pháp hiệu chỉnh áp suất. Giải pháp của tôi trở nên không ổn định sau một thời gian nhất định (khoảng 5s). …

12 fluid-dynamics