debug higher order

inputFiles/wavePropagation/benchmarks/acouselas3D_flu_komatitsch.xml

@@ -31,7 +31,7 @@
     <FiniteElements>
       <FiniteElementSpace
         name="FE1"
-        order="3"
+        order="1"
         formulation="SEM"/>
     </FiniteElements>
   </NumericalMethods>

src/coreComponents/physicsSolvers/wavePropagation/AcousticElasticWaveEquationSEM.cpp

@@ -39,12 +39,23 @@ void AcousticElasticWaveEquationSEM::initializePostInitialConditionsPreSubGroups
   auto acousNodesSet = acousticSolver()->getSolverNodesSet();
   auto elasNodesSet = elasticSolver()->getSolverNodesSet();
 
+  if (false)  // debug
+  {
+    std::cout << "acoustic nodes" << std::endl;
+    for (auto val : acousNodesSet) std::cout << val << ' ';
+    std::cout << std::endl;
+    std::cout << "elastic nodes" << std::endl;
+    for (auto val : elasNodesSet) std::cout << val << ' ';
+    std::cout << std::endl;
+  }
+
   for (auto val : acousNodesSet)
   {
     if (elasNodesSet.contains(val)) m_interfaceNodesSet.insert(val);
   }
 
-  // std::cout << "\t[AcousticElasticWaveEquationSEM::initializePostInitialConditionsPreSubGroups] size=" << m_solverTargetNodesSet.size() << std::endl;
+  std::cout << "\t[AcousticElasticWaveEquationSEM::initializePostInitialConditionsPreSubGroups] "
+  << "m_interfaceNodesSet.size()=" << m_interfaceNodesSet.size() << std::endl;
   GEOS_THROW_IF( m_interfaceNodesSet.size() == 0, "Failed to compute interface: check xml input (solver order)", std::runtime_error );
 }
 

src/coreComponents/physicsSolvers/wavePropagation/AcousticWaveEquationSEM.cpp

@@ -320,7 +320,7 @@ void AcousticWaveEquationSEM::initializePostInitialConditionsPreSubGroups()
       } );
     } );
   } );
-  // std::cout << "\t[AcousticWaveEquationSEM::initializePostInitialConditionsPreSubGroups] size=" << m_solverTargetNodesSet.size() << std::endl;
+  std::cout << "\t[AcousticWaveEquationSEM::initializePostInitialConditionsPreSubGroups] m_solverTargetNodesSet.size()=" << m_solverTargetNodesSet.size() << std::endl;
 
 }
 

src/coreComponents/physicsSolvers/wavePropagation/ElasticWaveEquationSEM.cpp

@@ -599,7 +599,7 @@ void ElasticWaveEquationSEM::initializePostInitialConditionsPreSubGroups()
     } );
   } );
 
-  // std::cout << "\t[ElasticWaveEquationSEM::initializePostInitialConditionsPreSubGroups] size=" << m_solverTargetNodesSet.size() << std::endl;
+  std::cout << "\t[ElasticWaveEquationSEM::initializePostInitialConditionsPreSubGroups] m_solverTargetNodesSet.size()=" << m_solverTargetNodesSet.size() << std::endl;
 }