Khoa học tính toán

Các Lax tương đương lý khẳng định rằng sự nhất quán và ổn định của một chương trình số cho một vấn đề giá trị ban đầu tuyến tính là một điều kiện cần và đủ để hội tụ. Nhưng đối với các bài toán phi tuyến, các phương pháp …

19 pde  stability  convergence 

Một câu trả lời gần đây đã đề cập đến việc sử dụng Máy tạo số ngẫu nhiên ( RNG ) của Fortuna hoặc Mersenne Twister để tạo ra một mô phỏng Monte Carlo . Tôi chưa nghe nói về Fortuna trước đây vì vậy tôi đã tra cứu nó …

19 algorithms  random-number-generation 

y0y0y_0β∈ Rβ∈R\beta \in \mathbb Ryyybạnbạnu12∥ y- y0∥2+ β2∥ u ∥212‖y-y0‖2+β2‖bạn‖2\begin{equation} \frac{1}{2} \| y - y_0 \|^2 + \frac{\beta}{2} \| u \|^2 \end{equation}Một y= u .Mộty= =bạn.\begin{equation} Ay = u. \end{equation}y, y0, U ∈ Rny,y0,bạn∈Rn y, y_0, u \in \mathbb R^n Một ∈ Rn × nMột∈Rn×n A \in \mathbb R^{n \times …

19 optimization  reference-request  time-integration 

làma trận thưa thớt xác định đối xứng dương tính (SPD) n × n . G là một ma trận đường chéo thưa thớt. n lớn ( n > 10000) và số lượng khác không trong G thường là 100 ~ 1000.AAAn×nn×nn \times nGGGnnnnnnGGG đã được factorized theo hình thức …

19 linear-algebra 

Đối với giải pháp của các hệ tuyến tính lớn sử dụng các phương pháp lặp, thường được quan tâm để đưa ra điều kiện tiên quyết, ví dụ giải quyết thay vì , trong đó ở đây được sử dụng cho tiền điều kiện bên trái của hệ thống. …

19 linear-algebra  iterative-method  preconditioning 

Tôi đang cố gắng giải quyết một số vấn đề tối ưu hóa phi tuyến không giới hạn trên GPU (CUDA). Hàm mục tiêu là một hàm phi tuyến trơn tru và độ dốc của nó tương đối rẻ để tính toán phân tích, vì vậy tôi không cần bận …

18 optimization  cuda 

Trong khi đọc một số tài liệu về người giải quyết PDE, tôi đã bắt gặp thuật ngữ giả bước thời gian ngày hôm nay. Nó dường như là một thuật ngữ phổ biến, tuy nhiên tôi đã thất bại trong việc tìm ra một định nghĩa hay một bài …

18 pde  time-integration 

Trong hầu hết (tất cả?) Việc thực hiện Phương pháp Đa cực nhanh (FMM), octrees được sử dụng để phân tách miền liên quan. Về mặt lý thuyết, octrees cung cấp một giới hạn thể tích đơn giản, rất hữu ích để chứng minh thời gian chạy O (n) của …

18 algorithms 

Tôi muốn tìm hiểu làm thế nào các phương pháp wavelet có thể được áp dụng cho PDE, nhưng tiếc là tôi không biết một tài nguyên tốt để tìm hiểu về chủ đề này. Dường như nhiều giới thiệu về wavelet tập trung vào lý thuyết nội suy, ví …

18 pde  wavelet 

Các thuật toán tốt nhất (khả năng mở rộng và hiệu quả) để tạo ra các lưới tứ giác không cấu trúc trong 2D là gì? Tôi có thể tìm một máy phát điện lưới bốn cấu trúc tốt ở đâu? (ưu tiên nguồn mở)

18 mesh-generation  computational-geometry 

Bạn muốn cải thiện bài này? Cung cấp câu trả lời chi tiết cho câu hỏi này, bao gồm trích dẫn và giải thích lý do tại sao câu trả lời của bạn là chính xác. Câu trả lời mà không đủ chi tiết có thể được chỉnh sửa hoặc …

18 pde  libraries  parallel-computing  adaptive-mesh-refinement 

Các đại lượng vật lý như áp suất, mật độ, năng lượng, nhiệt độ và nồng độ phải luôn dương, nhưng phương pháp số đôi khi tính toán các giá trị âm trong quá trình giải. Điều này không ổn vì các phương trình sẽ tính toán các giá trị …

18 pde  fluid-dynamics  hyperbolic-pde 

Tôi đang cố gắng thực hiện chức năng sau trong dấu phẩy động có độ chính xác kép với sai số tương đối thấp : l o g s u m (x,y) = nhật ký( điểm kinh nghiệm( x ) + điểm kinh nghiệm( y) )tôiogSbạnm(x,y)= =đăng nhập⁡(điểm kinh nghiệm⁡(x)+điểm …

18 floating-point  stability  numerics 

Tôi muốn sử dụng f2pyvới Fortran hiện đại. Cụ thể, tôi đang cố gắng lấy ví dụ cơ bản sau để làm việc. Đây là ví dụ hữu ích nhỏ nhất tôi có thể tạo ra. ! alloc_test.f90 subroutine f(x, z) implicit none ! Argument Declarations ! real*8, intent(in) :: …

18 python  fortran 

Giả sử tôi có một phương trình ma trận được định nghĩa đệ quy là A[n] = inverse([1 - b[n]A[n+1]]) * a[n] Sau đó, phương trình của A [1] trông tương tự như một phân số tiếp tục, trong đó có một số phương pháp hiệu quả cao để tránh …

18 algorithms