Câu hỏi được gắn thẻ «linear-solver»

Tôi muốn giải quyết một vấn đề phi tuyến tính với các ràng buộc đẳng thức phi tuyến tính và tôi đang sử dụng Lagrangian tăng cường với một thuật ngữ chính quy hóa hình phạt, như được biết đến, làm hỏng số điều kiện của các hệ thống tuyến …

12 linear-solver  preconditioning  constrained-optimization 

Các phương pháp Newton-Krylov (JFNK) của Jacobian và các phương pháp Krylov nói chung, có thể rất hữu ích vì chúng không yêu cầu lưu trữ rõ ràng hoặc xây dựng ma trận, chỉ kết quả của các sản phẩm vector ma trận. Nếu bạn thực sự hình thành hệ …

11 fluid-dynamics  linear-solver  preconditioning 

Có bất kỳ hy vọng nào trong việc giải quyết hệ thống tuyến tính sau một cách hiệu quả bằng phương pháp lặp không? A∈Rn×n,x∈Rn,b∈Rn, with n>106A∈Rn×n,x∈Rn,b∈Rn, with n>106A \in \mathbb{R}^{n \times n}, x \in \mathbb{R}^n, b \in \mathbb{R}^n \text{, with } n > 10^6 Ax=bAx=bAx=b với , nơi Δ …

11 linear-solver  dense-matrix  linear-system 

Cho hệ thống trong đó , tôi đọc rằng, trong trường hợp Jacobi lặp được sử dụng làm bộ giải, phương thức sẽ không hội tụ nếu có giá trị khác không thành phần trong null-không gian của . Vì vậy, làm thế nào người ta có thể chính thức …

11 linear-algebra  optimization  matrix  iterative-method  linear-solver 

Vui lòng liệt kê gói Python (petc4py, v.v ...) và các bộ giải trực tiếp thưa thớt mà nó hỗ trợ. Một câu trả lời (cộng đồng-wiki) cho mỗi gói, xin vui lòng.

11 python  linear-solver  sparse-matrix 

trong matlab, cả linsolve và mldivide đều được sử dụng để giải một hệ phương trình tuyến tính, trong tất cả các trường hợp xác định, quá hạn và không xác định. Đọc tài liệu của họ, tôi đã tự hỏi sự khác biệt giữa chúng là gì? Có phải …

10 matlab  linear-solver 

Tôi hiện đang làm việc để giải quyết các hệ thống đối xứng rất lớn (nhưng không tích cực xác định), được tạo ra bởi một số thuật toán nhất định. Các ma trận này có độ thưa thớt khối đẹp có thể được sử dụng để giải song song. …

10 parallel-computing  linear-solver  sparse-matrix  linear-system  gmres 

Tôi muốn biết bộ giải tuyến tính cổ điển nào (ví dụ Gauss-Seidel, Jacobi, SOR) được đảm bảo hội tụ cho bài toán trong đó A là bán xác định dương và dĩ nhiên b ∈ i m ( A )Ax=bAx=bAx=bAAAb∈im(A)b∈im(A)b \in im(A) (Lưu ý là bán xác định và …

10 linear-algebra  iterative-method  convergence  linear-solver 

Tôi đã phát triển một giải pháp hoạt động của Phương pháp phần tử hữu hạn để giải quyết các vấn đề truyền nhiệt bằng cách sử dụng GPU và OpenCL bằng phương pháp Conjugate Gradient. Nhược điểm chính của phương pháp này là nhu cầu cao về bộ nhớ. …

10 parallel-computing  opencl  linear-solver 

Tôi rất quan tâm đến việc tối ưu hóa địa ngục khỏi việc giải hệ thống tuyến tính cho các ma trận nhỏ (10 x 10), đôi khi được gọi là ma trận nhỏ . Có một giải pháp sẵn sàng cho việc này? Ma trận có thể được giả …

9 linear-algebra  linear-solver  compiling  open-source  assembly 

Tôi đang tìm kiếm một giải pháp rõ ràng (dám nói là tối ưu?), Vấn đề thực tuyến tính 3x3, , . Một ∈ R 3 × 3 , b ∈ R 3Ax=bAx=b\mathbf{A}\mathbf{x} = \mathbf{b}A∈R3×3,b∈R3A∈R3×3,b∈R3\mathbf{A} \in \mathbf{R}^{3 \times 3}, \mathbf{b} \in \mathbf{R}^{3} Ma trận là chung, nhưng gần với ma …

9 linear-solver  reference-request 

Có quy tắc đơn giản nào để nói nếu nó đáng để làm SOR thay vì Gauss-Seidel không? (và cách có thể để ước tính tham số realxation )ωω\omega Ý tôi là chỉ bằng cách nhìn vào ma trận , hoặc kiến ​​thức về vấn đề cụ thể mà ma …

9 linear-solver  iterative-method  heuristics 

Tôi dự định giải Ax = b trong đó A là ma trận vuông hoặc hình chữ nhật có điều kiện phức tạp, thưa thớt, không đối xứng và có điều kiện cao (số điều kiện ~ 1E + 20). Tôi đã có thể giải quyết hệ thống với ZGELSS …

9 linear-solver  sparse-matrix  lapack  condition-number 

Tôi cần phải giải Ax = b, nhưng tôi nhận ra rằng ngay cả khi nó thưa thớt, việc lưu trữ các hệ số ma trận của vấn đề của tôi sẽ chiếm quá nhiều bộ nhớ. Vì vậy, bây giờ tôi đang xem xét sử dụng phương pháp không …

9 linear-solver  petsc  sparse-operator 

Hầu hết các chương trình của tôi là mã nghiên cứu một lần trong C để sử dụng riêng. Tôi chưa bao giờ phân phối bất kỳ mã nào ngoài các cộng tác viên thân thiết. Tôi đã phát triển một thuật toán mà tôi đang xuất bản trong một …

9 linear-solver  libraries  c++  c