Scienza computazionale c++

2

risolvi

Sto trasferendo un codice esistente da MATLAB a C ++ e ho un sistema lineare per risolvere (piuttosto che la forma più tipica A x = b )xA=bxA=bxA=bAx=bAx=bAx=b La matrice è densa e di forma generale, ma non è più grande di 1000x1000. Quindi in MATLAB, la soluzione è trovata …

12 linear-solver c++ lapack

3

Libreria di algebra lineare Blaze?

L'articolo "Modelli di espressione rivisitati: un'analisi delle prestazioni delle attuali metodologie" nel SIAM Journal of Scientific Computing fa riferimento alla libreria di algebra lineare "Blaze". Non ne ho mai sentito parlare prima e non riesco a trovare riferimenti online. (Le ovvie ricerche su Google stanno restituendo il documento sopra.) Cos'è …

12 linear-algebra reference-request c++

1

Sostituendo l'integrazione QuasiMonteCarlo di Mathematica in C ++

Ho un programma Mathematica che esegue alcuni integrali in 3 o 4 dimensioni usando il QuasiMonteCarlometodo. Il problema è che richiede un tempo fastidiosamente lungo per essere eseguito, al punto che alcuni di questi calcoli non possono essere completati nel tempo di lavoro massimo disponibile sul nostro cluster HPC. Quindi …

12 libraries monte-carlo c++ mathematica

3

Implementazione efficiente dell'algoritmo a matrice tridiagonale

Sto risolvendo un problema fisico usando uno schema numerico implicito. Questo mi porta a risolvere un'equazione lineare con matrice tridiagonale. Ho codificato questo algoritmo da Wikipedia. Mi chiedo se esiste una libreria efficiente che consenta di risolvere questo tipo di equazione in modo ottimizzato. Una nota importante è che la …

12 linear-algebra libraries c++

1

Applicazione del metodo Runge-Kutta a ODE di secondo ordine

Come posso sostituire il metodo Euler con il Runge-Kutta del 4 ° ordine per determinare il movimento di caduta libera in magnitudo gravitazionale non costante (es. Caduta libera da 10 000 km dal suolo)? Finora ho scritto una semplice integrazione con il metodo Euler: while() { v += getMagnitude(x) * …

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

3

Verifica se due matrici 12x12 hanno lo stesso determinante

12×1212×1212 \times 12QQQdet(Q)=det(12I−Q−J)(1)det(Q)=det(12I−Q−J)(1)\det(Q) = \det(12I-Q-J) \; \; (1)JJJ Attualmente lo sto facendo con la libreria di armadillo ma risulta essere troppo lento. Il fatto è che devo farlo per un trilione di matrici e si scopre che il calcolo dei due determinanti è il collo di bottiglia del mio programma. …

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

5

Libreria C ++ per intergrazione numerica (quadratura)

Ho la mia piccola subroutine per l'integrazione numerica (quadratura), che è un adattamento C ++ di un programma ALGOL pubblicato da Bulirsch & Stoer nel 1967 (Numerische Mathematik, 9, 271-278). Vorrei passare a un algoritmo più moderno (adattivo) e chiedermi se ci sono librerie C ++ (gratuite) che forniscono tale. …

10 c++ quadrature

4

Alla ricerca di Runge-Kutta 8 ° ordine in C / C ++

Vorrei usare il metodo Runge-Kutta dell'8 ° ordine (89) in un'applicazione di meccanica / astrodinamica celeste, scritta in C ++, usando una macchina Windows. Quindi mi chiedo se qualcuno conosce una buona libreria / implementazione documentata e libera da usare? Va bene se è scritto in C, purché non ci …

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

1

Qual è l'impatto della semantica del movimento C ++ 11 nel contesto dell'informatica scientifica?

C ++ 11 introduce la semantica di spostamento che può, ad esempio, migliorare le prestazioni del codice in situazioni in cui C ++ 03 dovrebbe eseguire una costruzione o assegnazione di copie. Questo articolo riporta che il codice seguente presenta una velocità 5x quando viene compilato con C + 11: …

10 finite-element performance c++

3

Come creare una funzione spline ricorsiva in C ++

Al momento sto lavorando a un metodo di risoluzione delle equazioni differenziali chiamato collocazione base-spline. Ciò di cui sto avendo problemi è la costruzione di un metodo per costruire una spline di ordine arbitrario, con la relazione con la condizione iniziale B 1 i (x)={ 1Bk+1i(x)=x−xixk+i−xiBki+xk+i+1−xxk+i+1−xi+1Bki+1(x)Bik+1(x)=x−xixk+i−xiBik+xk+i+1−xxk+i+1−xi+1Bi+1k(x) B^{k+1}_{i}(x)= \frac{x-x_i}{x_{k+i}-x_i}B^k_i + \frac{x_{k+i+1}-x}{x_{k+i+1}-x_{i+1}}B^k_{i+1}(x) …

10 c++ b-spline

3

Ci sono pacchetti FEM “leggeri” in giro?

Fondamentalmente, la FEM sembra essere un problema praticamente "risolto". Esistono numerosi potenti framework, come Trilinos, PETSc, FEniCS, Libmesh o MOOSE. Una cosa che hanno in comune: sono estremamente "pesanti". Innanzitutto, l'installazione è normalmente molto dolorosa. In secondo luogo, la loro interfaccia / API è spessa e pesante: devi tradurre tutta …

9 finite-element python c++ petsc fenics

1

Implementazione del metodo Jacobi-Davidson per il problema degli autovalori cubici

Ho un grosso problema di autovalori cubici: (A0+λA1+λ2A2+λ3A3)x=0.(A0+λA1+λ2A2+λ3A3)x=0.\left(\mathbf{A}_0 + \lambda\mathbf{A}_1 + \lambda^2\mathbf{A}_2 + \lambda^3\mathbf{A}_3\right)\mathbf{x} = 0. Potrei risolverlo convertendomi in un problema di autovalore lineare ma si tradurrebbe in un sistema di dimensioni:32323^2 ⎡⎣⎢−A0000I000I⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥=λ⎡⎣⎢A1I0A20IA300⎤⎦⎥⎡⎣⎢xyz⎤⎦⎥,[−A0000I000I][xyz]=λ[A1A2A3I000I0][xyz],\begin{bmatrix} -\mathbf{A}_0 & 0 & 0 \\ 0 & \mathbf{I} & 0 \\ 0 & 0 & …

9 c++ eigensystem eigenvalues

10

È possibile ottimizzare questo codice di integrazione in modo che funzioni più velocemente?

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; } Quanto sopra è il mio codice …

9 c++ performance

2

Libreria C ++ per minimizzazione vincolata non lineare

Attualmente sto cercando di risolvere il problema di minimizzazione vincolata non lineare come implementato nella funzione "fmincon" di matlab. Le mie aspettative sono minimizzare (fun1, x0, uB, lB, fun2) dove x0 è lo stato iniziale, fun1 è la funzione che deve essere minimizzata, uB sono i limiti superiori, lB sono …

9 c++ nonlinear-programming

6

Domanda delicata: dove si inserisce Python nell'immagine?

Quindi ho discusso se dovrei preoccuparmi o meno di imparare Python. Dal parlare con i miei professori, Matlab sembra essere il linguaggio comune usato nella matematica applicata / scienze computazionali per quanto riguarda il mondo accademico ; mentre nell'industria , i miei professori (specialmente quelli che hanno lavorato nell'industria) hanno …

9 matlab python c++

Domande taggate «c++»