SW4  Git Revision 67ab458ccdb97763a04025b4bf840d7814d6d27f
SW4
 All Classes Pages
Public Types | Public Member Functions | Public Attributes | List of all members
EW Class Reference

Public Types

enum  InputMode {
  UNDEFINED, Efile, GaussianHill, GridFile,
  CartesianGrid, TopoImage, Rfile
}
 

Public Member Functions

 EW (const string &name, vector< Source * > &a_GlobalUniqueSources, vector< TimeSeries * > &a_GlobalTimeSeries, bool invproblem=false)
 
bool wasParsingSuccessful ()
 
bool isInitialized ()
 
void set_output_options (bool output_load, bool output_detailed_timing)
 
void setGMTOutput (string filename, string wppfilename)
 
void saveGMTFile (vector< Source * > &a_GlobalUniqueSources)
 
void allocateCartesianSolverArrays (double a_global_zmax)
 
void setGoalTime (double t)
 
void setNumberSteps (int steps)
 
int getNumberOfSteps () const
 
void setupRun (vector< Source * > &a_GlobalUniqueSources)
 
void solve (vector< Source * > &a_GlobalSources, vector< TimeSeries * > &a_GlobalTimeSeries)
 
void solve_backward (vector< Source * > &a_Sources, vector< TimeSeries * > &a_TimeSeries, double gradient[11], double hessian[121])
 
void solve_allpars (vector< Source * > &a_GlobalSources, vector< Sarray > &a_Rho, vector< Sarray > &a_Mu, vector< Sarray > &a_Lambda, vector< TimeSeries * > &a_GlobalTimeSeries, vector< Sarray > &a_U, vector< Sarray > &a_Um, vector< DataPatches * > &Upred_saved_sides, vector< DataPatches * > &Ucorr_saved_sides, bool save_sides)
 
void solve_backward_allpars (vector< Source * > &a_GlobalSources, vector< Sarray > &a_Rho, vector< Sarray > &a_Mu, vector< Sarray > &a_Lambda, vector< TimeSeries * > &a_GlobalTimeSeries, vector< Sarray > &a_U, vector< Sarray > &a_Um, vector< DataPatches * > &Upred_saved_sides, vector< DataPatches * > &Ucorr_saved_sides, double gradients[11], vector< Sarray > &gRho, vector< Sarray > &gMu, vector< Sarray > &gLambda)
 
bool parseInputFile (vector< Source * > &a_GlobalSources, vector< TimeSeries * > &a_GlobalTimeSeries)
 
void parsedate (char *datestr, int &year, int &month, int &day, int &hour, int &minute, int &second, int &msecond, int &fail)
 
void extractRecordData (TimeSeries::receiverMode mode, int i0, int j0, int k0, int grid0, vector< double > &uRec, vector< Sarray > &Um2, vector< Sarray > &U)
 
void badOption (string name, char *option) const
 
bool startswith (const char begin[], char *line)
 
void processGrid (char *buffer)
 
void processRefinement (char *buffer)
 
void deprecatedOption (const string &command, const string &oldone, const string &newone)
 
void processTime (char *buffer)
 
void processTwilight (char *buffer)
 
void processFileIO (char *buffer)
 
void processImage (char *buffer)
 
void processImage3D (char *buffer)
 
void deprecatedImageMode (int value, const char *name) const
 
void processTestPointSource (char *buffer)
 
void processTestRayleigh (char *buffer)
 
void processTestLamb (char *buffer)
 
void processTestEnergy (char *buffer)
 
bool checkTestEnergyPeriodic (char *buffer)
 
void processSource (char *buffer, vector< Source * > &a_GlobalUniqueSources)
 
void processRupture (char *buffer, vector< Source * > &a_GlobalUniqueSources)
 
void processMaterial (char *buffer)
 
void processMaterialIfile (char *buffer)
 
void processMaterialBlock (char *buffer, int &blockCount)
 
void processMaterialPfile (char *buffer)
 
void processMaterialEtree (char *buffer)
 
void processMaterialVimaterial (char *buffer)
 
void processMaterialInvtest (char *buffer)
 
void processMaterialRfile (char *buffer)
 
void processAnisotropicMaterialBlock (char *buffer, int &ablockCount)
 
void processReceiver (char *buffer, vector< TimeSeries * > &a_GlobalTimeSeries)
 
void processObservation (char *buffer, vector< TimeSeries * > &a_GlobalTimeSeries)
 
void processBoundaryConditions (char *buffer)
 
void processPrefilter (char *buffer)
 
void processGMT (char *buffer)
 
void processDeveloper (char *buffer)
 
void processGlobalMaterial (char *buffer)
 
void processTopography (char *buffer)
 
void processAttenuation (char *buffer)
 
void processRandomize (char *buffer)
 
void side_plane (int g, int side, int wind[6], int nGhost)
 
void setPrintCycle (int cycle)
 
void setVerbosity (int level)
 
void setQuiet (bool stealth)
 
bool getQuiet ()
 
int getVerbosity ()
 
int getRank ()
 
void setDebugIO (bool onoff)
 
void printTime (int cycle, double t, bool force=false) const
 
void printPreamble (vector< Source * > &a_Sources) const
 
void switch_on_checkfornan ()
 
void switch_on_error_log ()
 
void set_energylog (string logfile, bool print, bool elog)
 
void set_inner_loop (int loopnr)
 
void set_cflnumber (double cfl)
 
void set_testing_mode (bool a_testing)
 
bool get_testing_mode ()
 
void default_bcs ()
 
void update_curvilinear_cartesian_interface (vector< Sarray > &a_U)
 
void set_twilight_forcing (ForcingTwilight *a_forcing)
 
void initialData (double a_t, vector< Sarray > &a_U, vector< Sarray * > &a_AlphaVE)
 
bool exactSol (double a_t, vector< Sarray > &a_U, vector< Sarray * > &a_AlphaVE, vector< Source * > &source)
 
void exactRhsTwilight (double a_t, vector< Sarray > &a_F)
 
void exactAccTwilight (double a_t, vector< Sarray > &a_Uacc)
 
void Force (double a_t, vector< Sarray > &a_F, vector< GridPointSource * > point_sources)
 
void Force_tt (double a_t, vector< Sarray > &a_F, vector< GridPointSource * > point_sources)
 
void normOfDifference (vector< Sarray > &a_Uex, vector< Sarray > &a_U, double &diffInf, double &diffL2, double &xInf, vector< Source * > &a_globalSources)
 
void normOfDifferenceGhostPoints (vector< Sarray > &a_Uex, vector< Sarray > &a_U, double &diffInf, double &diffL2)
 
void normOfSurfaceDifference (vector< Sarray > &a_Uex, vector< Sarray > &a_U, double &diffInf, double &diffL2, double &solInf, double &solL2, vector< Source * > &a_globalSources)
 
void test_sources (vector< GridPointSource * > &a_point_sources, vector< Source * > &a_global_unique_sources, vector< Sarray > &F)
 
void testSourceDiscretization (int kx[3], int ky[3], int kz[3], double moments[3], vector< GridPointSource * > &point_sources, vector< Sarray > &F)
 
void setupSBPCoeff ()
 
void simpleAttenuation (vector< Sarray > &a_Up)
 
void enforceBC (vector< Sarray > &a_U, vector< Sarray > &a_Mu, vector< Sarray > &a_Lambda, double t, vector< double ** > &a_BCForcing)
 
void enforceBCfreeAtt (vector< Sarray > &a_Up, vector< Sarray > &a_U, vector< Sarray > &a_Um, vector< Sarray > &a_Mu, vector< Sarray > &a_Lambda, vector< Sarray * > &a_AlphaVEp, vector< Sarray * > &a_AlphaVEm, vector< double ** > &a_BCForcing, double bop[5], double a_t)
 
void enforceBCfreeAtt2 (vector< Sarray > &a_Up, vector< Sarray > &a_Mu, vector< Sarray > &a_Lambda, vector< Sarray * > &a_AlphaVEp, vector< double ** > &a_BCForcing)
 
void enforceBCanisotropic (vector< Sarray > &a_U, vector< Sarray > &a_C, double t, vector< double ** > &a_BCForcing)
 
void addAttToFreeBcForcing (vector< Sarray * > &AlphaVEp, vector< double ** > &BCForcing, double bop[5])
 
void cartesian_bc_forcing (double t, vector< double ** > &a_BCForcing, vector< Source * > &a_Source)
 
void evalRHS (vector< Sarray > &a_U, vector< Sarray > &a_Mu, vector< Sarray > &a_Lambda, vector< Sarray > &a_Lu, vector< Sarray * > &a_Alpha)
 
void evalRHSanisotropic (vector< Sarray > &a_U, vector< Sarray > &a_C, vector< Sarray > &a_Uacc)
 
void evalPredictor (vector< Sarray > &a_Up, vector< Sarray > &a_U, vector< Sarray > &a_Um, vector< Sarray > &a_Rho, vector< Sarray > &a_Lu, vector< Sarray > &a_F)
 
void evalDpDmInTime (vector< Sarray > &a_Up, vector< Sarray > &a_U, vector< Sarray > &a_Um, vector< Sarray > &a_Uacc)
 
void evalCorrector (vector< Sarray > &a_Up, vector< Sarray > &a_Rho, vector< Sarray > &a_Lu, vector< Sarray > &a_F)
 
void updateMemVarPred (vector< Sarray * > &a_AlphaVEp, vector< Sarray * > &a_AlphaVEm, vector< Sarray > &a_U, double a_t)
 
void updateMemVarCorr (vector< Sarray * > &a_AlphaVEp, vector< Sarray * > &a_AlphaVEm, vector< Sarray > &a_Up, vector< Sarray > &a_U, vector< Sarray > &a_Um, double a_t)
 
void updateMemVarCorrNearInterface (Sarray &a_AlphaVEp, Sarray &a_AlphaVEm, Sarray &a_Up, Sarray &a_U, Sarray &a_Um, double a_t, int a_mech, int a_grid)
 
void evalDpDmInTimeAtt (vector< Sarray * > &a_AlphaVEp, vector< Sarray * > &a_AlphaVE, vector< Sarray * > &a_AlphaVEm)
 
void addSuperGridDamping (vector< Sarray > &a_Up, vector< Sarray > &a_U, vector< Sarray > &a_Um, vector< Sarray > &a_Rho)
 
void cycleSolutionArrays (vector< Sarray > &a_Um, vector< Sarray > &a_U, vector< Sarray > &a_Up, vector< Sarray * > &a_AlphaVEm, vector< Sarray * > &a_AlphaVE, vector< Sarray * > &a_AlphaVEp)
 
void cycleSolutionArrays (vector< Sarray > &a_Um, vector< Sarray > &a_U, vector< Sarray > &a_Up)
 
void bndryInteriorDifference (vector< Sarray > &a_Uex, vector< Sarray > &a_U, double lowZ[3], double interiorZ[3], double highZ[3])
 
void test_RhoUtt_Lu (vector< Sarray > &a_Uacc, vector< Sarray > &a_Lu, vector< Sarray > &a_F, double lowZ[3], double interiorZ[3], double highZ[3])
 
void setRestartInfo (int fromCycle, int dumpInterval, const string &filePrefix)
 
void computeDT ()
 
void computeDTanisotropic ()
 
bool proc_zero () const
 
int no_of_procs () const
 
void create_output_directory ()
 
void initialize_image_files ()
 
void update_images (int Nsteps, double time, vector< Sarray > &a_Up, vector< Sarray > &a_U, vector< Sarray > &a_Um, vector< Sarray > &a_Rho, vector< Sarray > &a_Mu, vector< Sarray > &a_Lambda, vector< Source * > &a_sources, int dminus)
 
void initialize_SAC_files ()
 
void update_SACs (int Nsteps)
 
void print_execution_times (double times[7])
 
void print_execution_time (double t1, double t2, string msg)
 
void finalizeIO ()
 
string bc_name (const boundaryConditionType bc) const
 
int mkdirs (const string &path)
 
void setOutputPath (const string &path)
 
const string & getOutputPath ()
 
const string & getObservationPath ()
 
const string & getName ()
 
void set_global_bcs (boundaryConditionType bct[6])
 
void add_mtrl_block (MaterialData *md)
 
void set_threshold_velocities (double vpmin, double vsmin)
 
void addMaterialProperty (MaterialProperty *mat)
 
void processSupergrid (char *buffer)
 
void set_sg_damping (double coeff)
 
void set_sg_thickness (int gp_thickness)
 
void set_sg_width (double sg_width)
 
bool usingSupergrid ()
 
void setup_supergrid ()
 
void assign_supergrid_damping_arrays ()
 
void assign_local_bcs ()
 
bool timeSteppingSet ()
 
bool proc_decompose_2d (int ni, int nj, int nproc, int proc_max[2])
 
void decomp1d (int nglobal, int myid, int nproc, int &s, int &e)
 
void coarsen1d (int &n, int &ifirst, int &ilast, int periodic)
 
void allocateCurvilinearArrays ()
 
void generate_grid ()
 
void setup_metric ()
 
void check_min_max_int (vector< Sarray > &a_U)
 
void setupMPICommunications ()
 
void setup2D_MPICommunications ()
 
void communicate_array (Sarray &u, int grid)
 
void communicate_arrays (vector< Sarray > &u)
 
void communicate_array_2dfinest (Sarray &u)
 
void communicate_array_2d (Sarray &u, int g, int k)
 
void communicate_array_2d_asym (Sarray &u, int g, int k)
 
void communicate_array_2d_ext (Sarray &u)
 
void set_materials ()
 
void set_anisotropic_materials ()
 
void setup_attenuation_relaxation (double minvsoh)
 
void setup_viscoelastic ()
 
void setup_viscoelastic_tw ()
 
void extrapolateInZ (int g, Sarray &field, bool lowk, bool highk)
 
void extrapolateInXY (vector< Sarray > &field)
 
void extrapolateInZvector (int g, Sarray &field, bool lowk, bool highk)
 
void extrapolateInXYvector (vector< Sarray > &field)
 
void extrapolateTopo (Sarray &field)
 
void checkTopo (Sarray &field)
 
void addImage (Image *i)
 
void addImage3D (Image3D *i)
 
void setIO_timing (bool iotiming)
 
void setParallel_IO (bool pfs, int nwriters)
 
void extractTopographyFromGridFile (string a_topoFileName)
 
void extractTopographyFromImageFile (string a_topoFileName)
 
void extractTopographyFromCartesianFile (string a_topoFileName)
 
void setEtreeFile (EtreeFile *efile)
 
void extractTopographyFromEfile (string a_topoFileName, string a_topoExtFileName, string a_QueryType, double a_EFileResolution)
 
void extractTopographyFromRfile (std::string a_topoFileName)
 
void smoothTopography (int maxIter)
 
void buildGaussianHillTopography (double amp, double Lx, double Ly, double x0, double y0)
 
void extractSurfaceFromGridFile (string a_surfaceFileName)
 
void extractSurfaceFromCartesianFile (string a_surfaceFileName)
 
void computeCartesianCoord (double &x, double &y, double lon, double lat)
 
void computeGeographicCoord (double x, double y, double &longitude, double &latitude)
 
void initializeSystemTime ()
 
void compute_epicenter (vector< Source * > &a_GlobalUniqueSources)
 
void set_epicenter (double epiLat, double epiLon, double epiDepth, double earliestTime)
 
void get_epicenter (double &epiLat, double &epiLon, double &epiDepth, double &earliestTime)
 
void computeNearestGridPoint (int &a_i, int &a_j, int &a_k, int &a_g, double a_x, double a_y, double a_z)
 
void computeNearestSurfaceGridPoint (int &a_i, int &a_j, double a_x, double a_y, double a_z)
 
void coord_from_index (int i, int j, int k, int g, double &x, double &y, double &z)
 
double distance (double a_x1, double a_y1, double a_z1, double a_x0, double a_y0, double a_z0)
 
void computeNearestLowGridPoint (int &a_i, int &a_j, int &a_k, int &a_g, double a_x, double a_y, double a_z)
 
bool interior_point_in_proc (int a_i, int a_j, int a_g)
 
bool point_in_proc (int a_i, int a_j, int a_g)
 
bool point_in_proc_ext (int a_i, int a_j, int a_g)
 
void initializePaddingCells ()
 
void check_dimensions ()
 
void convert_material_to_mulambda ()
 
void check_materials ()
 
void computeSolutionError (vector< Sarray > &U, double t, vector< Sarray * > &Alpha)
 
double localMin (vector< Sarray > &a_field)
 
double localMax (vector< Sarray > &a_field)
 
double localMinVp ()
 
double localMaxVp ()
 
double localMinVs ()
 
double localMaxVs ()
 
double localMinVpOverVs ()
 
double localMaxVpOverVs ()
 
bool topographyExists ()
 
bool usingAttenuation ()
 
bool is_onesided (int g, int side) const
 
void interpolate_between_grids (vector< Sarray > &u, vector< Sarray > &um, double t, vector< Sarray * > &AlphaVE)
 
bool interpolate_topography (double q, double r, double &Z0, bool smoothed)
 
void copy_topo_to_topogridext ()
 
bool getDepth (double x, double y, double z, double &depth)
 
bool curvilinear_grid_mapping (double q, double r, double s, double &X0, double &Y0, double &Z0)
 
bool invert_curvilinear_grid_mapping (double X0, double Y0, double Z0, double &q, double &r, double &s)
 
bool find_curvilinear_derivatives_at_point (double q, double r, double s, double qX[], double rX[], double sX[])
 
void save_errors (double max_error[3], double l2_error[3])
 
void compute_minvsoverh (double &minvsoh)
 
void set_resolution (int ppw)
 
void set_prefilter (FilterType passband, int order, int passes, double fc1, double fc2)
 
void set_scenario (const string &scenario)
 
void set_conservative_interpolation (bool onoff, double ctol, int cmaxit)
 
void set_geodyn_data (string filename, int nx, int nz, double h, double origin[3], double dt, int nsteps, int faces)
 
void impose_geodyn_ibcdata (vector< Sarray > &u, vector< Sarray > &um, double t)
 
void get_geodyn_timelevel (vector< Sarray > &geodyndata)
 
void copy_geodyn_timelevel (vector< Sarray > &geodyndata1, vector< Sarray > &geodyndata2)
 
void integrate_source ()
 
void compute_energy (double dt, bool write_file, vector< Sarray > &Um, vector< Sarray > &U, vector< Sarray > &Up, int step)
 
double scalarProduct (vector< Sarray > &U, vector< Sarray > &V)
 
void get_gridgen_info (int &order, double &zetaBreak) const
 
void cleanUpRefinementLevels ()
 
int getNumberOfCartesianGrids ()
 
int getNumberOfGrids ()
 
int getNumberOfGhostPoints ()
 
int getNumberOfParallelPaddingPoints ()
 
double getLatOrigin ()
 
double getGridAzimuth ()
 
double getMetersPerDegree ()
 
bool usingParallelFS ()
 
int getNumberOfWritersPFS ()
 
double getTimeStep () const
 
int getNumberOfTimeSteps () const
 
void get_exact_point_source (double *u, double t, int g, Source &source, int *wind=0)
 
double VerySmoothBump_x_T_Integral (double t, double R, double alpha, double beta)
 
double C6SmoothBump_x_T_Integral (double t, double R, double alpha, double beta)
 
double SmoothWave_x_T_Integral (double t, double R, double alpha, double beta)
 
double Gaussian_x_T_Integral (double t, double R, double f, double alpha, double beta)
 
double VSBTP (double Lim, double t)
 
double C6SBTP (double Lim, double t)
 
double SWTP (double Lim, double t)
 
double d_VerySmoothBump_dt (double t, double R, double c)
 
double d_C6SmoothBump_dt (double t, double R, double c)
 
double d_SmoothWave_dt (double t, double R, double c)
 
double d_Gaussian_dt (double t, double R, double c, double f)
 
double VerySmoothBump (double t, double R, double c)
 
double C6SmoothBump (double t, double R, double c)
 
double SmoothWave (double t, double R, double c)
 
double Gaussian (double t, double R, double c, double f)
 
void get_exact_lamb (vector< Sarray > &a_U, double a_t, Source &a_source)
 
void get_exact_lamb2 (vector< Sarray > &a_U, double a_t, Source &a_source)
 
double G4_Integral (double T, double t, double r, double beta)
 
double G3_Integral (double iT, double it, double ir, double ibeta)
 
double G2_Integral (double iT, double it, double ir, double ibeta)
 
void getGlobalBoundingBox (double bbox[6])
 
string getPath ()
 
void set_utcref (TimeSeries &ts)
 
void print_utc ()
 
void processCG (char *buffer)
 
void processScaleFactors (char *buffer)
 
void average_speeds (double &cp, double &cs)
 
void layered_speeds (vector< double > &cp, vector< double > &z)
 
void testsourcediff (vector< Source * > GlobalSources, double gradient[11], double hessian[121])
 
void get_scalefactors (double sf[11])
 
bool compute_sf ()
 
void compute_guess (bool &guesspos, bool &guesst0fr, bool &guessmom, bool &guessshifts, bool &output_seismograms)
 
void get_cgparameters (int &maxit, int &maxrestart, double &tolerance, bool &fletcherreeves, int &stepselection, bool &do_linesearch, int &varcase, bool &testing)
 
void parameters_to_material (int nmpar, double *xm, vector< Sarray > &rho, vector< Sarray > &mu, vector< Sarray > &lambda)
 
void material_to_parameters (int nmpar, double *xm, vector< Sarray > &rho, vector< Sarray > &mu, vector< Sarray > &lambda)
 
void get_material_parameter (int nmpar, double *xm)
 
void get_scale_factors (int nmpar, double *xm)
 
void check_anisotropic_material (vector< Sarray > &rho, vector< Sarray > &c)
 
void get_nr_of_material_parameters (int &nmvar)
 
void add_to_grad (vector< Sarray > &K, vector< Sarray > &Kacc, vector< Sarray > &Um, vector< Sarray > &U, vector< Sarray > &Up, vector< Sarray > &Uacc, vector< Sarray > &gRho, vector< Sarray > &gMu, vector< Sarray > &gLambda)
 
void get_optmethod (int &method, int &bfgs_m)
 
void get_utc (int utc[7]) const
 
void perturb_mtrl ()
 
void perturb_mtrl (int peri, int perj, int perk, double h, int grid, int var)
 
void perturb_velocities (vector< Sarray > &a_vs, vector< Sarray > &a_vp)
 
void metric_derivatives_test ()
 
void material_ic (vector< Sarray > &a_mtrl)
 
void gettopowgh (double ai, double wgh[8]) const
 
void smooth_grid (int maxIter)
 
void enforceDirichlet5 (vector< Sarray > &a_U)
 
bool check_for_nan (vector< Sarray > &a_U, int verbose, string name)
 
void define_parallel_io (vector< Parallel_IO * > &parallel_io)
 
void read_volimage (std::string &path, std::string &fname, vector< Sarray > &data)
 
void interpolate (int nx, int ny, int nz, double xmin, double ymin, double zmin, double hx, double hy, double hz, Sarray &rho, Sarray &mu, Sarray &lambda, int grid, Sarray &rhogrid, Sarray &mugrid, Sarray &lambdagrid)
 
void interpolate_to_coarse (int nx, int ny, int nz, double xmin, double ymin, double zmin, double hx, double hy, double hz, Sarray &rho, Sarray &mu, Sarray &lambda, vector< Sarray > &rhogrid, vector< Sarray > &mugrid, vector< Sarray > &lambdagrid)
 
void interpolation_gradient (int nx, int ny, int nz, double xmin, double ymin, double zmin, double hx, double hy, double hz, Sarray &gradrho, Sarray &gradmu, Sarray &gradlambda, int grid, Sarray &gradrhogrid, Sarray &gradmugrid, Sarray &gradlambdagrid)
 
void enforceIC (std::vector< Sarray > &a_Up, std::vector< Sarray > &a_U, std::vector< Sarray > &a_Um, vector< Sarray * > &a_AlphaVEp, vector< Sarray * > &a_AlphaVE, vector< Sarray * > &a_AlphaVEm, double t, bool predictor, std::vector< GridPointSource * > point_sources)
 
void enforceIC2 (std::vector< Sarray > &a_Up, std::vector< Sarray > &a_U, std::vector< Sarray > &a_Um, vector< Sarray * > &a_AlphaVEp, double t, std::vector< GridPointSource * > point_sources)
 
void dirichlet_hom_ic (Sarray &U, int g, int k, bool inner)
 
void dirichlet_twilight_ic (Sarray &U, int g, int kic, double t)
 
void dirichlet_LRic (Sarray &U, int g, int kic, double t, int adj)
 
void dirichlet_LRstress (Sarray &B, int g, int kic, double t, int adj)
 
void gridref_initial_guess (Sarray &u, int g, bool upper)
 
void compute_preliminary_corrector (Sarray &a_Up, Sarray &a_U, Sarray &a_Um, Sarray *a_AlphaVEp, Sarray *a_AlphaVE, Sarray *a_AlphaVEm, Sarray &Utt, Sarray &Unext, int g, int kic, double t, std::vector< GridPointSource * > point_sources)
 
void compute_preliminary_predictor (Sarray &a_Up, Sarray &a_U, Sarray *a_AlphaVEp, Sarray &Unext, int g, int kic, double t, vector< GridPointSource * > point_sources)
 
void compute_icstresses (Sarray &a_Up, Sarray &B, int g, int kic, double *a_str_x, double *a_str_y)
 
void add_ve_stresses (Sarray &a_Up, Sarray &B, int g, int kic, int a_a, double *a_str_x, double *a_str_y)
 
void consintp (Sarray &Uf, Sarray &Unextf, Sarray &Bf, Sarray &Muf, Sarray &Lambdaf, Sarray &Rhof, double hf, Sarray &Uc, Sarray &Unextc, Sarray &Bc, Sarray &Muc, Sarray &Lambdac, Sarray &Rhoc, double hc, double cof, int gc, int gp, int is_periodic[2])
 
void checkintp (Sarray &Uf, Sarray &Unextf, Sarray &Bf, Sarray &Muf, Sarray &Lambdaf, Sarray &Rhof, double hf, Sarray &Uc, Sarray &Unextc, Sarray &Bc, Sarray &Muc, Sarray &Lambdac, Sarray &Rhoc, double hc, double cof, int gc, int gf, int is_periodic[2], double t)
 
void check_displacement_continuity (Sarray &Uf, Sarray &Uc, int gf, int gc)
 
void check_corrector (Sarray &Uf, Sarray &Uc, Sarray &Unextf, Sarray &Unextc, int kf, int kc)
 

Public Attributes

const double NO_TOPO
 
int mNumberOfGrids
 
int mNumberOfCartesianGrids
 
vector< double > mGridSize
 
vector< int > m_iStart
 
vector< int > m_iEnd
 
vector< int > m_jStart
 
vector< int > m_jEnd
 
vector< int > m_kStart
 
vector< int > m_kEnd
 
vector< int > m_iStartAct
 
vector< int > m_iEndAct
 
vector< int > m_jStartAct
 
vector< int > m_jEndAct
 
vector< int > m_kStartAct
 
vector< int > m_kEndAct
 
vector< int > m_iStartActGlobal
 
vector< int > m_iEndActGlobal
 
vector< int > m_jStartActGlobal
 
vector< int > m_jEndActGlobal
 
vector< int > m_kStartActGlobal
 
vector< int > m_kEndActGlobal
 
vector< int > m_global_nx
 
vector< int > m_global_ny
 
vector< int > m_global_nz
 
vector< int > m_iStartInt
 
vector< int > m_iEndInt
 
vector< int > m_jStartInt
 
vector< int > m_jEndInt
 
vector< int > m_kStartInt
 
vector< int > m_kEndInt
 
int m_paddingCells [4]
 
vector< double > m_zmin
 
Sarray mX
 
Sarray mY
 
Sarray mZ
 
Sarray mJ
 
Sarray mMetric
 
bool m_prefilter_sources
 
bool m_filter_observations
 
Filterm_filter_ptr
 
Filterm_filterobs_ptr
 
int m_opttest
 
Sarray mTopo
 
Sarray mTopoGridExt
 
vector< MaterialProperty * > m_materials
 
MPI_Comm m_cartesian_communicator
 
ofstream msgStream
 
vector< SarraymMu
 
vector< SarraymLambda
 
vector< SarraymRho
 
vector< SarraymC
 
Sarray mCcurv
 

The documentation for this class was generated from the following files: