Pertanyaan yang diberi tag «stability»

Saya memiliki beberapa masalah optimisasi global non-cembung yang menantang untuk dipecahkan. Saat ini saya menggunakan MATLAB's Optimization Toolbox (khusus, fmincon()dengan algoritma = 'sqp'), yang cukup efektif . Namun, sebagian besar kode saya menggunakan Python, dan saya ingin melakukan optimasi dengan Python juga. Apakah ada pemecah NLP dengan binding Python yang …

77 python  optimization  nonlinear-programming  linear-algebra  linear-solver  sparse-matrix  ode  automatic-differentiation  efficiency  matrix  symbolic-computation  pde  multiphysics  operator-splitting  fluid-dynamics  pde  multigrid  pde  hyperbolic-pde  quadrature  precision  numerical-analysis  fluid-dynamics  interpolation  c  ode  testing  data-sets  pde  fluid-dynamics  hyperbolic-pde  pde  finite-element  fluid-dynamics  stability  incompressible  pde  hyperbolic-pde  fluid-dynamics  pde  precision  floating-point  convex-optimization  conditioning  pde  hyperbolic-pde  stability  implicit-methods  explicit-methods  fluid-dynamics  accuracy  pde  hyperbolic-pde  petsc  linear-algebra 

The Lax ekuivalen teorema menyatakan bahwa konsistensi dan stabilitas skema numerik untuk masalah nilai awal linear adalah kondisi yang diperlukan dan cukup untuk konvergensi. Tetapi untuk masalah nonlinier, metode numerik dapat menyatu sangat masuk akal untuk hasil yang salah, meskipun konsisten dan stabil. Sebagai contoh, makalah ini menunjukkan bagaimana metode …

19 pde  stability  convergence 

Saya mencoba menerapkan fungsi berikut dalam floating point presisi ganda dengan kesalahan relatif rendah : logsum(x,y)=log(exp(x)+exp(y))logsum(x,y)=log⁡(exp⁡(x)+exp⁡(y))\mathrm{logsum}(x,y) = \log(\exp(x) + \exp(y)) Ini digunakan secara luas dalam aplikasi statistik untuk menambahkan probabilitas atau kepadatan probabilitas yang diwakili dalam ruang log. Tentu saja, baik atau dapat dengan mudah overflow atau underflow, yang akan …

18 floating-point  stability  numerics 

Saya ingin tahu apakah ada cara cepat untuk menghitung jarak Euclidean dari dua vektor dalam Oktaf. Sepertinya tidak ada fungsi khusus untuk itu, jadi haruskah saya menggunakan rumus saja sqrt?

16 octave  discretization  nonlinear-equations  newton-method  visualization  fluid-dynamics  mesh-generation  finite-element  finite-volume  optimization  algorithms  approximation  fluid-dynamics  navier-stokes  comsol  modeling  optimization  sparse-matrix  matrix  condition-number  visualization  matlab  quadrature  blas  intel-mkl  finite-element  gpu  discontinuous-galerkin  mathematica  optimization  convex-optimization  algorithms  reference-request  matlab  statistics  finite-element  numerical-analysis  petsc  molecular-dynamics  machine-learning  statistics  visualization  open-source  statistics  image-processing  visualization  python  petsc  finite-element  fluid-dynamics  stability  navier-stokes  incompressible 

Metode numerik untuk memecahkan PDE (atau ODE) terbagi dalam dua kategori besar: metode eksplisit dan implisit. Metode implisit memungkinkan timesteps stabil yang lebih besar tetapi membutuhkan lebih banyak pekerjaan per langkah. Untuk PDE hiperbolik, kebijaksanaan umum adalah bahwa metode implisit biasanya tidak membuahkan hasil karena penggunaan timesteps lebih besar dari …

16 pde  hyperbolic-pde  stability  implicit-methods  explicit-methods 

Setelah melakukan beberapa matematika yang berkaitan dengan stabilitas elemen dalam masalah 3D Stokes saya sedikit terkejut menyadari bahwa P2- P1P2-P1P_2-P_1 tidak stabil untuk mesh tetrahedral acak. Lebih tepatnya, jika Anda memiliki elemen di mana semua node dan tiga dari empat aspek terletak pada batas domain dengan kondisi Dirichlet, Anda akhirnya …

15 stability  navier-stokes  unstructured-mesh 

Saya menemukan komentar membingungkan di koran PJ van der Houwen, Pengembangan metode Runge-Kutta untuk persamaan diferensial parsial, Appl. Tidak. Matematika 20: 261, 1996 Pada baris 8ff di halaman 264, van der Houwen menulis: "Untuk polinomial Taylor ini menyiratkan bahwa interval stabilitas imajiner kosong untuk p=1,2,5,6,9,10,⋯p=1,2,5,6,9,10,⋯p = 1, 2, 5, 6, …

15 numerical-analysis  stability  runge-kutta 

Saya telah menerapkan pemecah Euler mundur dalam python 3 (menggunakan numpy). Untuk kenyamanan saya sendiri dan sebagai latihan, saya juga menulis fungsi kecil yang menghitung perkiraan perbedaan hingga dari gradien sehingga saya tidak selalu harus menentukan Jacobian secara analitis (jika itu mungkin!). Menggunakan deskripsi yang disediakan dalam Ascher dan Petzold …

13 pde  stability  numpy  scipy  newton-method 

Saya sedang mengerjakan pemecahan persamaan poroelastisitas satu dimensi (model biot), yang diberikan sebagai: −(λ+2μ)∂2u∂x2+∂p∂x=0−(λ+2μ)∂2u∂x2+∂p∂x=0-(\lambda+ 2\mu) \frac{\partial^2 u}{\partial x^2} + \frac{\partial p}{\partial x} = 0 ∂∂t[γp+∂u∂x]−κη[∂2p∂x2]=q(x,t)∂∂t[γp+∂u∂x]−κη[∂2p∂x2]=q(x,t)\frac{\partial}{\partial t} \left[ \gamma p + \frac{\partial u}{\partial x}\right] -\frac{\kappa}{\eta}\left[\frac{\partial^2 p}{\partial x^2}\right] =q(x,t) pada domain dan dengan syarat batas: Ω=(0,1)Ω=(0,1)\Omega=(0,1) p=0,(λ+2μ)∂u∂x=−u0p=0,(λ+2μ)∂u∂x=−u0p=0, (\lambda + 2\mu)\frac{\partial u}{\partial x}=-u_0 …

13 pde  finite-difference  stability  numerical-analysis 

Asumsikan saya memiliki fungsi bernilai nyata dari beberapa variabel x_i yang ingin saya evaluasi secara numerik. Secara umum rumus untuk f dapat mengandung produk, rasional, fungsi trancendental dll. Dan akan lama untuk menyelidiki stabilitas numerik secara analitik. Atau setidaknya akan memakan waktu untuk melakukannya dalam praktik. Asumsikan saya tidak memiliki …

12 stability 

Aliran bilangan Reynolds yang tinggi menghasilkan lapisan batas yang sangat tipis. Jika resolusi dinding digunakan dalam Simulasi Eddy Besar, rasio aspek mungkin berada di urutan . Banyak metode menjadi tidak stabil dalam rezim ini karena konstanta inf-sup menurun sebagai akar kuadrat dari rasio aspek atau lebih buruk. Konstanta inf-sup penting …

12 pde  finite-element  fluid-dynamics  stability  incompressible 

Saya sudah membaca itu (1) Operasi yang terkondisi harus dilakukan sebelum yang terkondisi dengan baik. Sebagai contoh, seseorang harus menghitung xz−yzxz−yzxz-yz sebagai(x−y)z(x−y)z(x-y)z karena pengurangan tidak dikondisikan sedangkan multiplikasi tidak. Namun, analisis kesalahan orde pertama dari kedua algoritma mengungkapkan bahwa mereka hanya berbeda oleh faktor tiga (*), dan saya tidak melihat …

10 stability  floating-point 

Saat ini saya memiliki kode yang menggunakan algoritma Crank-Nicholson, tapi saya pikir saya ingin pindah ke algoritma tingkat tinggi untuk mencatat waktu. Saya tahu bahwa algoritma Crank-Nicholson stabil di domain yang ingin saya gunakan, tetapi saya khawatir bahwa beberapa algoritma lain mungkin tidak. Saya tahu bagaimana cara menghitung kestabilan wilayah …

10 pde  stability  parabolic-pde 

Saya telah menggunakan Intel MKL SVD ( dgesvdmelalui SciPy) dan memperhatikan bahwa hasilnya sangat berbeda ketika saya mengubah presisi di antara float32dan float64ketika matriks saya dikondisikan secara buruk / bukan peringkat penuh. Apakah ada panduan tentang jumlah minimum regularisasi yang harus saya tambahkan untuk membuat hasil tidak sensitif terhadap float32-> …

10 numerical-analysis  stability  svd  intel-mkl  numerical-limitations 

Saya membaca makalah [1] di mana mereka memecahkan persamaan non-linear berikut menggunakan metode beda hingga. Mereka juga menganalisis stabilitas skema menggunakan analisis stabilitas Von Neumann. Namun, seperti yang disadari oleh penulis, ini hanya berlaku untuk PDE linier. Jadi penulis menyiasatinya dengan "pembekuan" istilah non-linear, yaitu mereka mengganti u u x …

9 pde  finite-difference  numerical-analysis  nonlinear-equations  stability