Khoa học tính toán

Đặt là ma trận thực, vuông, dày đặc. và đối xứng nhau. Để choG QA , G , QA,G,QA, G, QGGGQQQ H= [ A- Q- G- AT]H=[A−G−Q−AT]H = \begin{bmatrix} A & -G \\ -Q &-A^T \end{bmatrix} là một ma trận Hamilton. Tôi muốn để tính hàm mũ ma trận . …

10 linear-algebra  matrix  dense-matrix 

Trong bài báo http://www.scTHERirect.com/science/article/pii/S0045782509003521 , một phương trình phần tử-cục bộ HDG được mô tả trên trang 584 phương trình (4), với một trong các phương trình có dạng sau - ( bạnh, ∇ q)K= - ⟨ u^h⋅ n , q- q¯⟩∂K-(bạnh,∇q)K= =-⟨bạn^h⋅n,q-q¯⟩∂K-(u_h,\nabla q)_K = -\left\langle\hat{u}_h \cdot n, q - …

10 finite-element  discontinuous-galerkin 

Tôi đang tìm kiếm các phương pháp cho phép ước tính entropy thông tin của phân phối khi các cách lấy mẫu thực tế duy nhất từ ​​phân phối đó là phương pháp Monte Carlo. Vấn đề của tôi không giống như mô hình Ising tiêu chuẩn thường được sử …

10 monte-carlo  random-sampling 

Tôi có một hệ thống phương trình phi tuyến tính mà tôi muốn giải bằng số:nnn f = ( f 1 , Mạnh , f n )f(x)=af(x)=a\mathbf{f}(\mathbf{x})=\mathbf{a} f=(f1,…,fn)x=(x1,…,xn)f=(f1,…,fn)x=(x1,…,xn)\mathbf{f}=(f_1,\dots,f_n)\quad\mathbf{x}=(x_1,\dots,x_n) Hệ thống này có một số đặc điểm khiến nó đặc biệt khó xử lý. Tôi đang tìm kiếm ý tưởng về …

10 numerical-analysis  nonlinear-equations 

Tôi có một bộ dữ liệu đang dần thay đổi và tôi cần theo dõi các hàm riêng / giá trị riêng của ma trận hiệp phương sai của nó. Tôi đã sử dụng scipy.linalg.eigh, nhưng nó quá đắt, và nó không sử dụng thực tế là tôi đã có …

10 linear-algebra  optimization  python  eigenvalues 

Cái gọi là lưới quadtree và octree khá hấp dẫn cho các ứng dụng đòi hỏi phải tinh chỉnh lưới thích ứng. Chúng là ví dụ được sử dụng trong Gerris và Paramesh . Có ai biết một định dạng tệp tốt cho các lưới như vậy và hỗ trợ …

10 visualization 

Tôi có một mô hình khá lớn (~ 5000 dòng) được viết bằng C. Đây là một chương trình nối tiếp, không có tạo số ngẫu nhiên ở bất cứ đâu. Nó sử dụng thư viện FFTW cho các chức năng bằng cách sử dụng FFT - Tôi không biết …

10 floating-point  precision  computer-arithmetic 

Tôi biết rằng các bộ giải dòng chảy không nén và nén được thiết kế đặc biệt để giải quyết các loại vấn đề khác nhau với các tính chất / điều kiện lưu lượng khác nhau. Rõ ràng, trong số những lợi thế của việc sử dụng bộ giải …

10 fluid-dynamics 

Trong câu hỏi gần đây Giải hệ thống 7 phương trình đại số phi tuyến một cách tượng trưng , Brian Borchers đã xác nhận bằng thực nghiệm rằng Maple có thể giải một hệ đa thức mà Matlab / Mupad không thể xử lý. Trước đây tôi đã nghe …

10 reference-request  polynomials  symbolic-computation  maple 

Tôi đã cố gắng tạo ra một mô hình không gian 2D chính xác về mặt sinh học của các lớp mô, nơi xảy ra các quá trình sinh lý khác nhau. Điều này bao gồm chủ yếu là các phản ứng hóa học, khuếch tán và thông lượng trên …

10 finite-element  comsol  computational-biology 

Tôi có chương trình con nhỏ của riêng mình để tích hợp số (bậc hai), đó là bản chuyển thể C ++ của chương trình ALGOL do Bulirsch & Stoer xuất bản năm 1967 (Numerische Mathematik, 9, 271-278). Tôi muốn nâng cấp lên một thuật toán (thích ứng) hiện đại …

10 c++  quadrature 

Không gian Hyperbolic trong mô hình không gian nửa trên Poincaré trông giống như RnRn\Bbb R^n thông thường nhưng với khái niệm góc và khoảng cách bị bóp méo một cách tương đối đơn giản. Trong không gian Euclide tôi có thể lấy mẫu một điểm ngẫu nhiên thống nhất …

10 computational-geometry  monte-carlo  geometry  random-sampling 

Trong câu trả lời của tôi cho một câu hỏi về MSE liên quan đến mô phỏng vật lý 2D Hamilton, tôi đã đề nghị sử dụng một bộ tích hợp đối xứng bậc cao . Sau đó, tôi nghĩ rằng có thể là một ý tưởng tốt để chứng …

10 python  time-integration 

Tôi đã đọc nó (1) Các hoạt động có điều kiện nên được thực hiện trước khi hoạt động tốt. Ví dụ, người ta nên tính xz−yzxz-yzxz-yz là vì phép trừ được điều hòa trong khi phép nhân không.(x−y)z(x-y)z(x-y)z Tuy nhiên, phân tích lỗi thứ nhất của cả hai thuật …

10 stability  floating-point 

Có các thiết lập phương thức phần tử hữu hạn cung cấp các ước tính lỗi trong định mức (nghĩa là giới hạn trên ) không? Những họ các yếu tố có thể được sử dụng để thực hiện chúng?W1 , ∞W1,∞W^{1,\infty}∥ u'h- bạn'∥∞‖bạnh'-bạn'‖∞\|u'_h - u'\|_\infty ( Crossposted từ MathOverflow, …

10 finite-element  error-estimation