Ilmu Komputasi

Saya ingin tahu bagaimana kondisi Dirichlet biasanya diterapkan ketika menggunakan metode volume hingga pada grid non-seragam berpusat sel, Implementasi saya saat ini hanya memaksakan kondisi batas saya memperbaiki nilai sel pertama, ϕ1= gD( xL.)ϕ1=gD(xL.) \phi_1 = g_D(x_L) di mana adalah variabel solusi dan adalah nilai kondisi batas Dirichlet di lhs …

10 boundary-conditions  finite-volume  poisson 

Algoritma Parareal, PITA, dan PFASST semuanya adalah teknik lintas domain untuk memparalelkan solusi dari masalah yang bergantung pada waktu. Apa prinsip panduan di balik metode ini? Apa perbedaan utama di antara mereka? Dapatkah saya mengatakan bahwa satu didasarkan pada yang lain? Bagaimana? Bagaimana dengan aplikasi mereka? Saya tahu tidak akan …

10 parallel-computing  time-integration 

Untuk reduksi model, saya ingin menghitung vektor singular kiri terkait dengan - katakanlah 20 - nilai singular terbesar dari matriks , di mana N ≈ 10 6 dan k ≈ 10 3 . Sayangnya, matriks A saya akan padat tanpa struktur apa pun.A∈RN,kA∈RN,kA \in \mathbb R^{N,k}N≈106N≈106N\approx 10^6k≈103k≈103k\approx 10^3AAA Jika saya …

10 linear-algebra  python  reference-request  numpy  svd 

Saya mencari untuk port beberapa kode yang menyelesaikan satu set persamaan diferensial parsial (PDE) dengan metode volume hingga dalam bentuk IMPLICIT (untuk diskritisasi waktu). Akibatnya ada sistem persamaan tridiagonal dalam arah x, y, z yang ditangani oleh skema ADI / TDMA. Saya tidak bisa menemukan apa pun mengenai solusi implisit …

10 parallel-computing  implicit-methods  cuda 

Apa perbedaan antara FEM eksplisit dan FEM implisit? Menurut posting di sini , tampaknya satu-satunya perbedaan adalah apakah integrasi waktu implisit atau eksplisit digunakan. Seperti yang saya ingat dari satu buku yang saya baca, FEM implisit adalah di mana massa tidak disamakan dengan node. Apa definisi pasti FEM eksplisit dan …

10 finite-element  implicit-methods  explicit-methods 

Saya baru mengenal komputasi ilmiah. Saya mencari implementasi Fortran (lebih disukai f90) dari Octree. Masalah saya memerlukan Oktree yang membagi domain saya sampai tidak ada lebih dari beberapa partikel N (atau sumber di mana saya tahu nilai kerapatan yang dapat dicolokkan dalam metode persamaan integral untuk menyelesaikan persamaan seperti persamaan …

10 fortran  integral-equations 

Saya mencoba menyelesaikan persamaan 2D Poisson menggunakan metode Discontinuous Galerkin (DG) dan diskritisasi berikut (saya punya file png tapi saya tidak diizinkan mengunggahnya, maaf): Persamaan: ∇⋅(κ∇T)+f=0∇⋅(κ∇T)+f=0\nabla \cdot( \kappa \nabla T) + f = 0 Persamaan baru: q=κ∇T∇⋅q=−fq=κ∇T∇⋅q=−fq = \kappa \nabla T\\\nabla \cdot q = -f Bentuk lemah dengan numerik fluks …

10 pde  finite-element  fenics 

Dalam metode solusi yang diproduksi (MMS) seseorang mendalilkan solusi yang tepat, menggantikannya dalam persamaan dan menghitung istilah sumber yang sesuai. Solusinya kemudian digunakan untuk verifikasi kode. Untuk persamaan Navier-Stokes yang tidak dapat dimampatkan, MMS dengan mudah menyebabkan istilah sumber (bukan nol) dalam persamaan kontinuitas. Tetapi tidak semua kode memungkinkan istilah …

10 navier-stokes  incompressible  verification 

Sebagaimana diketahui bahwa aliran yang tidak dapat dimampatkan tidak ada dalam kenyataan, ini merupakan asumsi yang diperkenalkan untuk menyederhanakan persamaan yang mengatur. Kami tidak dapat menerapkan asumsi ini secara langsung. Secara umum bilangan Mach (M <0,3 untuk aliran yang tidak dapat dimampatkan), variasi kepadatan (variasi kepadatan nol) dan divergensi kecepatan …

10 fluid-dynamics 

Saya menggunakan skema beda hingga Crank-Nicolson untuk menyelesaikan persamaan panas 1D. Saya bertanya-tanya apakah prinsip maksimum / minimum dari persamaan panas (yaitu bahwa maksimum / minimum terjadi pada kondisi awal atau pada batas-batas) juga berlaku untuk solusi diskritisasi. Ini mungkin tersirat oleh fakta bahwa Crank-Nicolson adalah skema yang stabil dan …

10 linear-algebra  pde  finite-difference  crank-nicolson 

C ++ 11 memperkenalkan langkah semantik yang dapat, misalnya, meningkatkan kinerja kode dalam situasi di mana C ++ 03 perlu melakukan konstruksi salin atau tugas penyalinan. Artikel ini melaporkan bahwa kode berikut mengalami kecepatan 5x saat dikompilasi dengan C + 11: vector<vector<int> > V; for(int k = 0; k < …

10 finite-element  performance  c++ 

Kita semua akrab dengan banyak metode komputasi untuk menyelesaikan sistem linear standar Ax=b.Ax=b. Ax=b. Namun, saya ingin tahu apakah ada metode komputasi "standar" untuk memecahkan sistem linier yang lebih umum (dimensi terbatas) dari formulir tersebut. A m 1 × n 1 B m 2 × n 2 L m 1 …

10 linear-algebra 

Tutup. Pertanyaan ini di luar topik . Saat ini tidak menerima jawaban. Ingin meningkatkan pertanyaan ini? Perbarui pertanyaan agar sesuai dengan topik untuk Computational Science Stack Exchange. Ditutup 2 tahun yang lalu . Seperti judulnya, saya ingin dapat menghentikan kode di lokasi tertentu dan menghentikannya dengan aman. Saya tidak dapat …

10 matlab 

dalam matlab, baik linsolve dan mldivide digunakan untuk menyelesaikan sistem persamaan linear, dalam semua kasus yang ditentukan, ditentukan secara berlebihan dan tidak ditentukan. Membaca dokumen mereka, saya bertanya-tanya apa perbedaan di antara mereka? Apakah mereka menggunakan algoritma faktorisasi matriks dan triangularisasi yang hampir sama dalam tiga kasus? Jika A memiliki …

10 matlab  linear-solver 

Saya memiliki masalah pengoptimalan yang terlihat seperti berikut ini minJ,Bs.t.∑ij|Jij|MJ+BY=XminJ,B∑ij|Jij|s.t.MJ+BY=X \begin{array}{rl} \min_{J,B} & \sum_{ij} |J_{ij}|\\ \textrm{s.t.} & MJ + BY =X \end{array} Di sini, variabel saya adalah matriks JJJ dan BBB , tetapi seluruh masalahnya masih berupa program linier; variabel yang tersisa diperbaiki. Ketika saya mencoba memasukkan program ini ke …

10 optimization  convex-optimization  linear-programming