GEOS-DEV / GEOS

GEOS Simulation Framework
GNU Lesser General Public License v2.1
210 stars 84 forks source link

[Bug] Integrated tests failing on Lassen #1201

Closed corbett5 closed 3 years ago

corbett5 commented 3 years ago

The following integrated tests are failing on Lassen with clang

 FAIL RUN : 6
 (  KGD_ZeroViscosity_06 Sneddon_01 PhaseFieldFracture_DamageOnly_01 PhaseFieldFracture_DamageOnly_08 PhaseFieldFracture_DamageAndLoad_01 PhaseFieldFracture_DamageAndLoad_08  ) 

 FAIL CHECK : 1
 (  Hydrofracture_KGD_NodeBased_C3D6_09  ) 

(https://github.com/GEOSX/nightlyTests/blob/4d868a5a92cb38d48725ede20d8c9cbb1f9b463d/lassen/lassen_integrated_test.log)

Of these it appears that the PhaseFieldFracture problems use mathpresso which isn't available for PowerPC. If they could be rewritten without mathpresso that would most likely fix the issue.

Sneddon_01 has the following run time error

***** ERROR
***** LOCATION: /usr/WS2/corbett5/geosx/src/coreComponents/linearAlgebra/interfaces/trilinos/TrilinosSolver.cpp:100
***** Controlling expression (should be false): ierr != 0
***** Rank 0: Error in call to solver->NumericFactorization()
Expected ierr == 0
  ierr = -22
  0 = 0

** StackTrace of 9 frames **
Frame 0: geosx::TrilinosSolver::solve(geosx::EpetraMatrix&, geosx::EpetraVector&, geosx::EpetraVector&, geosx::DofManager const*) 
Frame 1: geosx::SolverBase::SolveSystem(geosx::DofManager const&, geosx::EpetraMatrix&, geosx::EpetraVector&, geosx::EpetraVector&) 
Frame 2: geosx::SolverBase::NonlinearImplicitStep(double const&, double const&, int, geosx::DomainPartition&) 
Frame 3: geosx::SolidMechanicsEmbeddedFractures::SolverStep(double const&, double const&, int, geosx::DomainPartition&) 
Frame 4: geosx::SolverBase::Execute(double, double, int, int, double, geosx::dataRepository::Group*) 
Frame 5: geosx::EventBase::Execute(double, double, int, int, double, geosx::dataRepository::Group*) 
Frame 6: geosx::EventManager::Run(geosx::dataRepository::Group*) 
Frame 7: geosx::ProblemManager::RunSimulation() 
Frame 8: main 
Frame 9:  
=====

And Hydrofracture_KGD_NodeBased_C3D6_09 could be fixed by upping the absolute tolerance by 2.5x to 0.0005.

corbett5 commented 3 years ago

The following fail with gcc

FAIL RUN:
Hydrofracture_KGD_NodeBased_C3D6_01 KGD_ZeroToughness_01 KGD_ZeroToughness_06 KGD_ZeroViscosity_01 
KGD_ZeroViscosity_06 Sneddon_conforming01 Sneddon_conforming04 SurfaceGenerator_01 SurfaceGenerator_08 
DryFrac_StaticPenny_PrismElem_01 DryFrac_StaticPenny_PrismElem_08 DryFrac_ThreeNodesPinched_HorizontalFrac_01 
DryFrac_ThreeNodesPinched_SlantFrac_01 co2_flux_3d_01 Sneddon_01 PhaseFieldFracture_DamageOnly_01 
PhaseFieldFracture_DamageOnly_08 PhaseFieldFracture_DamageAndLoad_01 PhaseFieldFracture_DamageAndLoad_08 
PoroElastic_Terzaghi_01 FlowProppantTransport_2d_01 staircase_3d_01 beamBending_01 beamBending_08 beamBending_27 

 FAIL CHECK : 1
 Hydrofracture_KGD_NodeBased_C3D6_09
corbett5 commented 3 years ago

Don't need two issues for this. https://github.com/GEOSX/GEOSX/issues/1372