Khoa học tính toán

Tôi biết rằng ATLAS có thể tự tối ưu hóa cho máy được biên dịch và do đó lợi ích tối đa được tìm thấy bằng cách biên dịch từ nguồn. Có bất kỳ lợi ích nào khi biên dịch LAPACK từ nguồn không? Sẽ dễ dàng hơn nhiều nếu …

13 performance  lapack 

Thông thường, tôi viết mã nối tiếp và khi tôi làm, tôi viết các bài kiểm tra đơn vị với một số khung kiểm tra kiểu xUnit (MATLAB xUnit, PyUnit / mũi hoặc khung kiểm tra C ++ của Google). Dựa trên một tìm kiếm đáng chú ý của Google, …

13 parallel-computing  petsc  mpi  testing  reproducibility 

Tôi đang cố gắng giải phương trình MHD lý tưởng bằng các phương pháp bán rời, tái tạo không gian ENO và bước thời gian TVD RK. Tôi đang nhận được các giải pháp trạng thái ổn định khác nhau với thứ tự thời gian khác nhau. Nó có đúng …

13 fluid-dynamics 

Hàm có số ít gần . Điểm kỳ dị đó có thể được nâng lên, mặc dù: với , người ta nên có , vì Và do đó Tuy nhiên, hình thức không chỉ không được xác định tại , nó cũng không ổn định về số lượng trong vùng …

13 c++  c 

Tiêu đề là câu hỏi. Kỹ thuật này liên quan đến việc sử dụng "ma trận các đồng yếu tố" hoặc "ma trận liên hợp" và đưa ra các công thức rõ ràng cho các thành phần nghịch đảo của ma trận vuông. Thật không dễ dàng để làm bằng …

13 linear-algebra  matrix  complexity 

Gần đây, tôi đã thúc đẩy nhóm của mình bao gồm nhiều thử nghiệm hơn khi viết mã của họ. Có một số lỗi lớn cần nhiều thời gian để bắt hơn là có thể nói là cần thiết, bởi vì chúng tôi không có chế độ kiểm tra tốt. …

13 testing 

Tôi quan tâm đến việc thực hiện một lưới di chuyển cho một vấn đề khuếch tán. Các phương pháp lưới di chuyển thích ứng đưa ra một ví dụ tốt về cách thực hiện điều này cho phương trình của Burgers trong 1D bằng cách sử dụng sai phân …

13 finite-difference  adaptive-mesh-refinement  advection-diffusion 

Tôi phải giải các bài toán eigenvalue tổng quát trong đó và đều là tam giác, là đối xứng dương xác định và thực, nhưng chỉ là đối xứng phức tạp (không xác định hoặc Hermiti). Hơn nữa, tôi cần sự xuất tinh đầy đủ. Tôi hiện chỉ đang gọi …

13 linear-algebra  eigensystem 

Tôi đọc một số tài liệu tham khảo bao gồm cả điều này . Tôi hơi bối rối không biết vấn đề tối ưu hóa đã được xây dựng và cố gắng giải quyết. Là nó minimizesubject to∥x∥1Ax=bminimize‖x‖1subject toAx=b\begin{array}{ll} \text{minimize} & \|x\|_1\\ \text{subject to} & Ax=b\end{array} hoặc và minimizesubject to∥x∥0Ax=bminimize‖x‖0subject …

13 optimization 

Tôi có sự nhầm lẫn này về quy tắc Armijo được sử dụng trong tìm kiếm dòng. Tôi đã đọc lại tìm kiếm dòng theo dõi nhưng không hiểu quy tắc Armijo này là gì. Bất cứ ai cũng có thể xây dựng quy tắc Armijo là gì? Các wikipedia …

13 optimization 

Câu hỏi: Những phương pháp nào có sẵn để tính toán chính xác và hiệu quả cấu trúc thưa thớt của ma trận phần tử hữu hạn? Thông tin: Tôi đang làm việc trên một bộ giải phương trình áp suất Poisson, sử dụng phương pháp của Galerkin với cơ …

13 matrix  finite-element  petsc  performance 

Theo Wikipedia , tốc độ hội tụ được biểu thị bằng tỷ lệ cụ thể của các chỉ tiêu vectơ. Tôi đang cố gắng tìm hiểu sự khác biệt giữa tốc độ "tuyến tính" và "bậc hai", tại các thời điểm khác nhau (về cơ bản, "ở đầu" của lần …

13 iterative-method  convergence 

Trong chương trình mẫu này, tôi đang làm điều tương tự (ít nhất là tôi nghĩ vậy) theo hai cách khác nhau. Tôi đang chạy cái này trên máy tính Linux của tôi và theo dõi việc sử dụng bộ nhớ hàng đầu. Sử dụng gfortran tôi thấy rằng theo …

13 performance  fortran 

Tôi đang làm việc để giải các phương trình poroelasticity một chiều (mô hình của biot), được đưa ra như sau: trên miềnvà với các điều kiện biên: - ( λ + 2 μ ) ∂2bạn∂x2+ ∂p∂x= 0-(λ+2μ)∂2bạn∂x2+∂p∂x= =0-(\lambda+ 2\mu) \frac{\partial^2 u}{\partial x^2} + \frac{\partial p}{\partial x} = 0 Ω=(0,1)∂∂t[ …

13 pde  finite-difference  stability  numerical-analysis 

\newcommand{\v}[1]{\boldsymbol{#1}} Giả sử chúng ta có phương trình mô hình dòng Stokes: {−div(ν∇u)+∇pdivu=f=0{−div(ν∇u)+∇p=fdivu=0 \tag{1} \left\{ \begin{aligned} -\mathrm{div}(\nu \nabla \v{u}) + \nabla p &= \v{f} \\ \mathrm{div} \v{u} &= 0 \end{aligned} \right. Trong đó độ nhớt ν(x)ν(x)\nu(x) là một hàm, đối với phần tử hữu hạn hỗn hợp tiêu chuẩn, giả …

13 finite-element