Scienza computazionale

Domande e risposte per gli scienziati che utilizzano i computer per risolvere problemi scientifici

4
Quadratura numerica con derivati
La maggior parte dei metodi numerici per la quadratura considera l'integrando come una funzione di scatola nera. E se avessimo maggiori informazioni? In particolare, quali benefici possiamo derivare dalla conoscenza dei primi derivati ​​dell'integrando? Quali altre informazioni potrebbero essere utili? In particolare per i derivati: le stime di errore per …
6
Qual è il modo migliore per effettuare ricerche riproducibili se hai bisogno di librerie proprietarie?
La ricerca riproducibile nel calcolo mira a rendere il codice necessario per generare i risultati in un documento di calcolo disponibile per altri ricercatori in modo che possano eseguire questo codice per riprodurre i risultati in quel documento. Mi piacerebbe rendere riproducibili tutte le mie ricerche, ma sto incontrando un …
2
Che cos'è lo pseudo time-stepping?
Mentre leggevo della letteratura sui solutori della PDE, mi sono imbattuto oggi nel termine pseudo-passo temporale . Sembra essere un termine comune, tuttavia non sono riuscito a trovare una buona definizione o un articolo introduttivo per esso. Quindi: che cos'è lo pseudo time-stepping e come viene solitamente utilizzato?
1
Come si possono applicare le wavelet al PDE?
Vorrei imparare come applicare i metodi wavelet al PDE, ma sfortunatamente non conosco una buona risorsa per conoscere questo argomento. Sembra che molte introduzioni alle wavelet si concentrino sulla teoria dell'interpolazione, ad esempio assemblando un segnale mediante una sovrapposizione di preferibilmente poche wavelet. Le domande di PDE sono talvolta menzionate, …
18 pde  wavelet 
3
Risoluzione di problemi di ottimizzazione non lineare non vincolati su GPU
Sto cercando di risolvere alcuni problemi di ottimizzazione non lineare non vincolati su GPU (CUDA). La funzione obiettivo è una funzione non lineare regolare, e il suo gradiente è relativamente economico da calcolare analiticamente, quindi non ho bisogno di preoccuparmi dell'approssimazione numerica. Voglio risolvere questo problema con la maggior parte …
4
Esiste una libreria per scopi generici per il perfezionamento della mesh adattativa a griglia strutturata?
Vuoi migliorare questo post? Fornisci risposte dettagliate a questa domanda, comprese le citazioni e una spiegazione del perché la tua risposta è corretta. Le risposte senza dettagli sufficienti possono essere modificate o eliminate. Il raffinamento di maglie adattative (AMR) è una tecnica comune per affrontare il problema delle scale spaziali …
2
F2Py con array di forme allocabili e ipotizzate
Mi piacerebbe usare f2pycon il moderno Fortran. In particolare, sto cercando di far funzionare il seguente esempio di base. Questo è il più piccolo esempio utile che potrei generare. ! alloc_test.f90 subroutine f(x, z) implicit none ! Argument Declarations ! real*8, intent(in) :: x(:) real*8, intent(out) :: z(:) ! Variable …
18 python  fortran 
1
Cancellazione catastrofica nel logum
Sto cercando di implementare la seguente funzione in virgola mobile a precisione doppia con errore relativo basso : 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)) Questo è ampiamente utilizzato nelle applicazioni statistiche per aggiungere probabilità o densità di probabilità rappresentate nello spazio del registro. Ovviamente, o potrebbero facilmente overflow o underflow, il …
2
Esiste un algoritmo efficiente per le frazioni continue valutate a matrice?
Supponiamo di avere un'equazione di matrice definita ricorsivamente come A[n] = inverse([1 - b[n]A[n+1]]) * a[n] Quindi l'equazione di A [1] sembra simile a una frazione continua, per la quale esistono alcuni metodi altamente efficienti che evitano il noioso ricalcolo (vedere "Ricette numeriche" per alcuni esempi). Tuttavia, mi chiedo se …
18 algorithms 
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.