Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

fix: an issue of established f-m connections being duplicated in the run time. #3466

Merged
merged 12 commits into from
Dec 19, 2024
Original file line number Diff line number Diff line change
Expand Up @@ -176,7 +176,8 @@
int const cycleNumber,
DomainPartition & domain )
{
if( cycleNumber == 0 && time_n <= 0 )
// for initial fracture initialization
if( time_n <= 0 )

Check warning on line 180 in src/coreComponents/physicsSolvers/multiphysics/HydrofractureSolver.cpp

View check run for this annotation

Codecov / codecov/patch

src/coreComponents/physicsSolvers/multiphysics/HydrofractureSolver.cpp#L180

Added line #L180 was not covered by tests
{
initializeNewFractureFields( domain );
}
Expand Down Expand Up @@ -366,11 +367,11 @@
maxHydraulicAperture = MpiWrapper::max( maxHydraulicAperture );

GEOS_LOG_LEVEL_INFO_RANK_0( logInfo::Solution, GEOS_FMT( " {}: Max aperture change: {} m, max hydraulic aperture change: {} m",
this->getName(), fmt::format( "{:.{}f}", maxApertureChange, 6 ), fmt::format( "{:.{}f}", maxHydraulicApertureChange, 6 ) ) );
this->getName(), fmt::format( "{:.{}e}", maxApertureChange, 6 ), fmt::format( "{:.{}e}", maxHydraulicApertureChange, 6 ) ) );
GEOS_LOG_LEVEL_INFO_RANK_0( logInfo::Solution, GEOS_FMT( " {}: Min aperture: {} m, max aperture: {} m",
this->getName(), fmt::format( "{:.{}f}", minAperture, 6 ), fmt::format( "{:.{}f}", maxAperture, 6 ) ) );
this->getName(), fmt::format( "{:.{}e}", minAperture, 6 ), fmt::format( "{:.{}e}", maxAperture, 6 ) ) );
GEOS_LOG_LEVEL_INFO_RANK_0( logInfo::Solution, GEOS_FMT( " {}: Min hydraulic aperture: {} m, max hydraulic aperture: {} m",
this->getName(), fmt::format( "{:.{}f}", minHydraulicAperture, 6 ), fmt::format( "{:.{}f}", maxHydraulicAperture, 6 ) ) );
this->getName(), fmt::format( "{:.{}e}", minHydraulicAperture, 6 ), fmt::format( "{:.{}e}", maxHydraulicAperture, 6 ) ) );
}
template< typename POROMECHANICS_SOLVER >
void HydrofractureSolver< POROMECHANICS_SOLVER >::setupCoupling( DomainPartition const & domain,
Expand Down Expand Up @@ -737,7 +738,6 @@
for( localIndex kf=0; kf<2; ++kf )
{
localIndex const faceIndex = elemsToFaces[kfe][kf];

for( localIndex a=0; a<numNodesPerFace; ++a )
{

Expand Down Expand Up @@ -1160,7 +1160,6 @@
}
} );
}

subRegion.m_recalculateConnectionsFor2dFaces.clear();
subRegion.m_newFaceElements.clear();
} );
Expand Down
Loading