diff --git a/src/coreComponents/physicsSolvers/wavePropagation/sem/acoustic/secondOrderEqn/isotropic/AcousticWaveEquationSEM.cpp b/src/coreComponents/physicsSolvers/wavePropagation/sem/acoustic/secondOrderEqn/isotropic/AcousticWaveEquationSEM.cpp index 1fc032d98cf..655438384dc 100644 --- a/src/coreComponents/physicsSolvers/wavePropagation/sem/acoustic/secondOrderEqn/isotropic/AcousticWaveEquationSEM.cpp +++ b/src/coreComponents/physicsSolvers/wavePropagation/sem/acoustic/secondOrderEqn/isotropic/AcousticWaveEquationSEM.cpp @@ -246,8 +246,6 @@ void AcousticWaveEquationSEM::addSourceToRightHandSide( integer const & cycleNum arrayView2d< real64 const > const sourceConstants = m_sourceConstants.toViewConst(); arrayView1d< localIndex const > const sourceIsAccessible = m_sourceIsAccessible.toViewConst(); arrayView2d< real32 const > const sourceValue = m_sourceValue.toViewConst(); - //bool useSourceWaveletTables = m_useSourceWaveletTables; - //arrayView1d< TableFunction::KernelWrapper const > const sourceWaveletTableWrappers = m_sourceWaveletTableWrappers.toViewConst(); GEOS_THROW_IF( cycleNumber > sourceValue.size( 0 ), getDataContext() << ": Too many steps compared to array size", @@ -256,14 +254,10 @@ void AcousticWaveEquationSEM::addSourceToRightHandSide( integer const & cycleNum { if( sourceIsAccessible[isrc] == 1 ) { - //real64 const srcValue = - // useSourceWaveletTables ? sourceWaveletTableWrappers[ isrc ].compute( &time_n ) : WaveSolverUtils::evaluateRicker( time_n, - // timeSourceFrequency, timeSourceDelay, rickerOrder ); for( localIndex inode = 0; inode < sourceConstants.size( 1 ); ++inode ) { real32 const localIncrement = sourceConstants[isrc][inode] * sourceValue[cycleNumber][isrc]; - //real32 const localIncrement = sourceConstants[isrc][inode] * srcValue; RAJA::atomicAdd< ATOMIC_POLICY >( &rhs[sourceNodeIds[isrc][inode]], localIncrement ); } } @@ -288,7 +282,7 @@ void AcousticWaveEquationSEM::initializePostInitialConditionsPreSubGroups() - forDiscretizationOnMeshTargets( domain.getMeshBodies(), [&] ( string const &, + forDiscretizationOnMeshTargets( domain.getMeshBodies(), [&] ( string const & meshBodyName, MeshLevel & mesh, arrayView1d< string const > const & regionNames ) { @@ -376,7 +370,7 @@ void AcousticWaveEquationSEM::initializePostInitialConditionsPreSubGroups() MeshLevel & mesh, arrayView1d< string const > const & regionNames ) { - + GEOS_UNUSED_VAR( meshBodyName ); MeshLevel & baseMesh = domain.getMeshBodies().getGroup< MeshBody >( meshBodyName ).getBaseDiscretization(); precomputeSourceAndReceiverTerm( baseMesh, mesh, regionNames ); diff --git a/src/coreComponents/physicsSolvers/wavePropagation/shared/WaveSolverBase.cpp b/src/coreComponents/physicsSolvers/wavePropagation/shared/WaveSolverBase.cpp index db80f5fb2d5..342033418f0 100644 --- a/src/coreComponents/physicsSolvers/wavePropagation/shared/WaveSolverBase.cpp +++ b/src/coreComponents/physicsSolvers/wavePropagation/shared/WaveSolverBase.cpp @@ -58,7 +58,7 @@ WaveSolverBase::WaveSolverBase( const std::string & name, setInputFlag( InputFlags::FALSE ). setRestartFlags( RestartFlags::NO_WRITE ). setSizedFromParent( 0 ). - setDescription( "Source Value of the sources" ); + setDescription( "Array which contains the value of the Ricker wavelets at each time-steps" ); registerWrapper( viewKeyStruct::timeSourceDelayString(), &m_timeSourceDelay ).