Acoustic FDTD Solver
|
Simulate the propagation of the elastic wave using the finite-difference time-domain method. More...
#include <ELFDTD22.hpp>
Public Member Functions | |
ELFDTD22 () | |
void | calcCoefficients (Model &model, const Grid &grid, Function *signal, const double &cfl_number) |
See notes on 2012-05-14 for the expression of the coefficients. | |
void | execute (const Grid &grid, const CartPosition &sig_loc, Function *signal, const double &time_span) |
See notes on 2012-05-14 for the expression of the update equations. | |
void | attach (const unsigned &j, const unsigned &k) |
std::vector< std::vector < double > > | detach () |
~ELFDTD22 () |
Private Attributes | |
double ** | yya |
Coefficients of the stress components. | |
double ** | yyb |
double ** | zza |
double ** | zzb |
double ** | yza |
double ** | yzb |
double ** | uya |
Coefficients of the velocity components. | |
double ** | uyb |
double ** | uza |
double ** | uzb |
std::vector< unsigned > | sensor_j |
std::vector< unsigned > | sensor_k |
std::vector< std::vector < double > > | sensor |
double | dt |
temporal step |
Simulate the propagation of the elastic wave using the finite-difference time-domain method.
Use the velocity-stress formulation.
ELFDTD22::ELFDTD22 | ( | ) |
ELFDTD22::~ELFDTD22 | ( | ) |
void ELFDTD22::attach | ( | const unsigned & | j, |
const unsigned & | k | ||
) |
void ELFDTD22::calcCoefficients | ( | Model & | model, |
const Grid & | grid, | ||
Function * | signal, | ||
const double & | cfl_number | ||
) |
See notes on 2012-05-14 for the expression of the coefficients.
std::vector< std::vector< double > > ELFDTD22::detach | ( | ) |
void ELFDTD22::execute | ( | const Grid & | grid, |
const CartPosition & | sig_loc, | ||
Function * | signal, | ||
const double & | time_span | ||
) |
See notes on 2012-05-14 for the expression of the update equations.
From the grid and the cartposition to calculate the location of the hard source.
Pre-allocate the memory space for the field quantities
Current field - refers to the current (n for p and n+1/2 for u, )
Precedent field - refers to the one occuring immediately before the current (n-1 for p and n-1/2 for u )
|
private |
temporal step
|
private |
|
private |
|
private |
|
private |
Coefficients of the velocity components.
|
private |
|
private |
|
private |
|
private |
Coefficients of the stress components.
|
private |
|
private |
|
private |
|
private |
|
private |