Nella ricerca sui problemi inversi, è comune costruire un set di dati sintetici da un set noto di parametri e quindi verificare se la tecnica di inversione può ricostruire tali parametri. Nel fare ciò, è importante aggiungere livelli adeguati di rumore casuale ai dati sintetici. Inoltre, se il metodo utilizzato …
Ho un set di dati che sta lentamente cambiando e devo tenere traccia degli autovettori / autovalori della sua matrice di covarianza. Ho usato scipy.linalg.eigh, ma è troppo costoso e non usa il fatto che ho già una decomposizione che è solo leggermente errata. Qualcuno può suggerire un approccio migliore …
Nel documento http://www.sciencedirect.com/science/article/pii/S0045782509003521 , un'equazione locale-elemento HDG è descritta a pagina 584 equazione (4), con una delle equazioni che assume la forma seguente - ( uh, ∇ q)K= - ⟨ u^h⋅n , q-q¯⟩∂K-(uh,∇q)K=-⟨u^h⋅n,q-q¯⟩∂K-(u_h,\nabla q)_K = -\left\langle\hat{u}_h \cdot n, q - \bar{q}\right\rangle_{\partial K} Qual è l'approssimazione variazionale all'equazione continua , con …
Sto cercando metodi che consentano di stimare l'entropia delle informazioni di una distribuzione quando le uniche modalità pratiche di campionamento da quella distribuzione sono i metodi Monte Carlo. Il mio problema non è diverso dal modello Ising standard che viene generalmente utilizzato come esempio introduttivo per il campionamento di Metropolis-Hastings. …
Ho un sistema di equazioni non lineari che voglio risolvere numericamente:nnn f = ( f 1 , … , f n )f( x ) = af(X)=un'\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) Questo sistema ha una serie di caratteristiche che lo rendono particolarmente difficile da …
So che i solutori di flusso incomprimibili e comprimibili sono specificamente progettati per risolvere diversi tipi di problemi con diverse proprietà del fluido / condizioni di flusso. Chiaramente, tra i vantaggi dell'utilizzo di solutori di flusso incomprimibili per modellare problemi con fluidi incomprimibili è che l'equazione dell'energia può essere trascurata, …
Ho un modello considerevole (~ 5000 righe) scritto in C. È un programma seriale, senza generazione casuale di numeri da nessuna parte. Fa uso della libreria FFTW per le funzioni che usano FFT - Non conosco i dettagli dell'implementazione di FFTW, ma presumo che anche le funzioni ivi contenute siano …
Le cosiddette quadtree e octree griglie sono abbastanza attraenti per le applicazioni che richiedono adaptive mesh refinement. Sono ad esempio utilizzati in Gerris e Paramesh . Qualcuno è a conoscenza di un buon formato di file per tali griglie e supporta il software di visualizzazione? Vedi anche questa risposta nelle …
Ho cercato di realizzare un modello spaziale 2D biologicamente accurato di strati di tessuto, in cui si verificano diversi processi fisiologici. Ciò include principalmente reazioni chimiche, diffusione e flussi oltre i confini. Sto realizzando questo modello in COMSOL Multiphysics, un pacchetto software ad elementi finiti che risolve diversi tipi di …
Nella domanda recente Risolvendo simbolicamente il sistema di 7 equazioni algebriche non lineari , Brian Borchers ha confermato sperimentalmente che Maple può risolvere un sistema polinomiale che Matlab / Mupad non può gestire. Ho sentito in passato da persone che lavorano sul campo che Maple ha un'implementazione di alta qualità …
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. …
Esistono configurazioni del metodo degli elementi finiti che forniscono stime di errore nella norma (ovvero, limiti su )? Quali famiglie di elementi possono essere utilizzate per implementarle?W1 , ∞W1,∞W^{1,\infty}∥ u'h- u'∥∞‖uh'-u'‖∞\|u'_h - u'\|_\infty ( Crossposted da MathOverflow, dove ha riscontrato scarso interesse, ma probabilmente qui posso trovare più persone con …
L'ho letto (1) Le operazioni mal condizionate devono essere eseguite prima di quelle ben condizionate. Ad esempio, si dovrebbe calcolare come ( x - y ) z poiché la sottrazione è mal condizionata mentre la moltiplicazione no.x z- yzxz−yzxz-yz( x - y) z(x−y)z(x-y)z Tuttavia, un'analisi degli errori del primo ordine …
Lo spazio iperbolico nel modello del semispazio superiore di Poincaré sembra un normale RnRn\Bbb R^n ma con la nozione di angolo e distanza distorta in modo relativamente semplice. Nello spazio euclideo posso campionare un punto casuale uniformemente in una palla in diversi modi, ad esempio generando nnn campioni gaussiani indipendenti …
Qualche tempo fa avevo scritto un codice che tentava di calcolare il senza usare le funzioni di libreria. Ieri stavo rivedendo il vecchio codice e ho cercato di renderlo il più veloce possibile (e corretto). Ecco il mio tentativo finora:l o g( x )log(x)log(x) const double ee = exp(1); double …
We use cookies and other tracking technologies to improve your browsing experience on our website,
to show you personalized content and targeted ads, to analyze our website traffic,
and to understand where our visitors are coming from.
By continuing, you consent to our use of cookies and other tracking technologies and
affirm you're at least 16 years old or have consent from a parent or guardian.