«c++» etiketlenmiş sorular

C ++ için MATLAB mevcut bir kodu taşıma ve çözmek için bir lineer sistem bilgisi am (yerine daha tipik bir şekilde daha A x = b )x A = bxbir=bxA=bA x = bbirx=bAx=b matrisi yoğun ve genel formda olmakla birlikte 1000x1000'den daha büyük değildir. MATLAB'da çözüm, işlev veya eğik çizgi …

12 linear-solver  c++  lapack 

SIAM Bilimsel Hesaplama Dergisi'nde yer alan "İfade Şablonları Yeniden Ziyaret Edildi: Güncel Metodolojilerin Performans Analizi" makalesi "Blaze" doğrusal cebir kütüphanesine atıfta bulunmaktadır. Daha önce duymadım ve çevrimiçi referanslar bulamıyorum. (Açıkça görülen google aramaları yukarıdaki makaleyi geri veriyor.) Peki bu kütüphane nedir ve nereden daha fazla bilgi edinebilirim?

12 linear-algebra  reference-request  c++ 

QuasiMonteCarloYöntemi kullanarak 3 veya 4 boyutlu bazı integraller gerçekleştiren bir Mathematica program var . Sorun şu ki, bu hesaplamaların bazılarının HPC kümemizde bulunan maksimum iş süresinde tamamlanamayacağı noktaya kadar can sıkıcı bir şekilde uzun zaman alıyor. Bu yüzden programı büyük bir faktörle hızlandıracağından şüphelendiğim C ++ 'da yeniden yazmayı düşünüyorum. …

12 libraries  monte-carlo  c++  mathematica 

Fiziksel bir problemi örtük sayısal şema kullanarak çözüyorum. Bu beni tridiagonal matris ile doğrusal bir denklemin çözülmesine götürüyor. Ben kodlu ettik bu algoritma Wikipedia'dan. Bu tür bir denklemin optimize edilmiş bir şekilde çözülmesini sağlayan etkili bir kütüphane olup olmadığını merak ediyorum. Önemli bir not, matrisin sadece sistem parametreleri değiştiğinde değişmesidir, …

12 linear-algebra  libraries  c++ 

Sabit yerçekimi büyüklüğünde olmayan serbest düşme hareketini belirlemek için Euler yöntemini Runge-Kutta 4. sıra ile nasıl değiştirebilirim (örn. Yerden 10 000 km yukarıda serbest düşüş)? Şimdiye kadar Euler yöntemi ile basit entegrasyon yazdım: while() { v += getMagnitude(x) * dt; x += v * dt; time += dt; } x …

11 numerical-analysis  c++  runge-kutta  integral-equations 

Bana 12×1212×1212 \times 12 matris QQQ , bu simetrik, ters çevrilebilir, pozitif kesin ve yoğun. Eğer test mi det(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1) burada JJJ tamamen birli bir matristir. Şu anda armadillo kütüphanesi ile yapıyorum ama çok yavaş olduğu ortaya çıkıyor. Mesele şu ki, bunu bir trilyon matris için …

11 linear-algebra  matrix  c++  graph-theory  c 

1967'de Bulirsch & Stoer tarafından yayınlanan bir ALGOL programının C ++ uyarlaması olan sayısal entegrasyon (kareleme) için kendi küçük alt programım var (Numerische Mathematik, 9, 271-278). Daha modern (uyarlanabilir) bir algoritmaya yükseltmek ve böyle bir (ücretsiz) C ++ kütüphanesi olup olmadığını merak ediyorum. GSL (ki bu C) olarak bir göz …

10 c++  quadrature 

Runge-Kutta 8. dereceden yöntemini (89) C ++ ile yazılmış bir gök mekaniği / astrodinamik uygulamasında bir Windows makinesi kullanarak kullanmak istiyorum. Bu yüzden herkes belgelenmiş ve kullanımı ücretsiz iyi bir kütüphane / uygulama biliyor mu merak ediyorum? Beklenilen herhangi bir derleme problemi olmadığı sürece C ile yazılmışsa sorun yoktur. Şimdiye …

10 c++  reference-request  libraries  c  runge-kutta 

C ++ 11, örneğin C ++ 03'ün bir kopya oluşturma veya kopyalama ataması gerçekleştirmesi gereken durumlarda kod performansını artırabilecek taşıma semantiği sunar. Bu makalede , aşağıdaki kod C + 11 ile derlendiğinde 5 kat hızlanma yaşadığı bildirilmektedir: vector<vector<int> > V; for(int k = 0; k < 100000; ++k) { vector<int> …

10 finite-element  performance  c++ 

Şu anda temel-spline kollokasyonu adı verilen bir diferansiyel denklem çözme yöntemi üzerinde çalışıyorum. Ne ile sorun yaşıyorum ile ilişkisi ile keyfi bir sipariş spline oluşturmak için bir yöntem inşa etmektir ile başlangıç ​​koşulu ve bu sorunla başlamakta bile sorun yaşıyorum, çünkü özyinelemeli "üst" veya "alt" ile başlayabilir ve genel bir …

10 c++  b-spline 

Temel olarak, FEM hemen hemen "çözülmüş" bir sorun gibi görünüyor. Trilinos, PETSc, FEniCS, Libmesh veya MOOSE gibi çok sayıda güçlü çerçeve var. Ortak noktalarından biri: Son derece “ağır”. İlk olarak, kurulum normalde süper ağrılıdır. İkincisi, arayüzleri / API'ları kalın ve ağırdır - tüm fikrinizi ilgili kütüphanenin düşüncesine çevirmeniz gerekir. Bu, …

9 finite-element  python  c++  petsc  fenics 

double trap(double func(double), double b, double a, double N) { double j; double s; double h = (b-a)/(N-1.0); //Width of trapezia double func1 = func(a); double func2; for (s=0,j=a;j<b;j+=h){ func2 = func(j+h); s = s + 0.5*(func1+func2)*h; func1 = func2; } return s; } Yukarıdaki trapezi kullanarak func()sınırları arasındaki 1D …

9 c++  performance 

Büyük bir kübik özdeğer problemim var: (bir0+ λbir1+λ2bir2+λ3bir3) x =0.(bir0+λbir1+λ2bir2+λ3bir3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Bunu lineer bir özdeğer problemine dönüştürerek çözebilirdim ama bir sistemle sonuçlanırdı 32323^2 büyük gibi: ⎡⎣⎢-bir0000ben000ben⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥= λ⎡⎣⎢bir1ben0bir20benbir300⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥,[-bir0000ben000ben][xyz]=λ[bir1bir2bir3ben000ben0][xyz],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 \\ 0 & \mathbf{I} & 0 \\ 0 & 0 & …

9 c++  eigensystem  eigenvalues 

Şu anda matlab "fmincon" fonksiyonunda uygulandığı gibi doğrusal olmayan kısıtlı minimizasyon problemini çözmeye çalışıyorum. Beklentilerim, minimize etmek (fun1, x0, uB, lB, fun2), burada x0 başlangıç ​​durumudur, fun1 en aza indirilmesi gereken işlevdir, uB üst sınırlardır, lB daha düşük sınırlardır ve fun2, doğrusal olmayan eşitliklerin vektörlerini sağlayan işlevdir / http://www.mathworks.com/help/optim/ug/fmincon.html adresinde …

9 c++  nonlinear-programming 

Bu yüzden Python öğrenmeye zahmet etmem gerekip gerekmediğini tartışıyorum. Profesörlerimle konuştuğumdan, Matlab uygulamalı matematik / bilgisayar bilimlerinde akademi söz konusu olduğunda kullanılan ortak dil gibi görünüyor ; ederken sanayi , benim profesör (özellikle. endüstrisinde çalışmış olanlar) söylediler öğrenme c ++ en güvenli yoldur. Hem akademi hem de endüstride, Python ile …

9 matlab  python  c++