Skip to content
Open
Show file tree
Hide file tree
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Original file line number Diff line number Diff line change
Expand Up @@ -34,6 +34,7 @@ using namespace geos;
using namespace geos::dataRepository;
using namespace geos::constitutive;
using namespace geos::testing;
using namespace geos::wells::kernels::compositional;

CommandLineOptions g_commandLineOptions;
void writeTableToFile( string const & filename, char const * str )
Expand Down Expand Up @@ -492,7 +493,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DPRES,
wellElemDofNumber[iwelem] + ColOffset::DPRES,
dP,
jacobianFD.toViewConstSizes() );
}
Expand All @@ -513,7 +514,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + jc,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + jc,
dRho,
jacobianFD.toViewConstSizes() );
}
Expand All @@ -538,7 +539,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + NC,
dRate,
jacobianFD.toViewConstSizes() );
}
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -34,6 +34,7 @@ using namespace geos;
using namespace geos::dataRepository;
using namespace geos::constitutive;
using namespace geos::testing;
using namespace geos::wells::kernels::compositional;

CommandLineOptions g_commandLineOptions;
void writeTableToFile( string const & filename, char const * str )
Expand Down Expand Up @@ -541,7 +542,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DPRES,
wellElemDofNumber[iwelem] + ColOffset::DPRES,
dP,
jacobianFD.toViewConstSizes() );
}
Expand All @@ -562,7 +563,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + jc,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + jc,
dRho,
jacobianFD.toViewConstSizes() );
}
Expand All @@ -584,7 +585,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC+1,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + NC+1,
dT,
jacobianFD.toViewConstSizes() );
}
Expand All @@ -610,7 +611,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + NC,
dRate,
jacobianFD.toViewConstSizes() );
}
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -34,6 +34,7 @@ using namespace geos;
using namespace geos::dataRepository;
using namespace geos::constitutive;
using namespace geos::testing;
using namespace geos::wells::kernels::compositional;

CommandLineOptions g_commandLineOptions;

Expand Down Expand Up @@ -395,7 +396,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells<> & solver,

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DPRES,
wellElemDofNumber[iwelem] + ColOffset::DPRES,
dP,
jacobianFD.toViewConstSizes() );
}
Expand All @@ -416,7 +417,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells<> & solver,

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + jc,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + jc,
dRho,
jacobianFD.toViewConstSizes() );
}
Expand All @@ -441,7 +442,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells<> & solver,

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + NC,
dRate,
jacobianFD.toViewConstSizes() );
}
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -37,6 +37,7 @@ using namespace geos;
using namespace geos::dataRepository;
using namespace geos::constitutive;
using namespace geos::testing;
using namespace geos::wells::kernels::singlePhase;

CommandLineOptions g_commandLineOptions;

Expand Down Expand Up @@ -287,7 +288,7 @@ void testNumericalJacobian( SinglePhaseReservoirAndWells<> & solver,
// this is computing J_RW and J_WW
fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + singlePhaseWellKernels::ColOffset::DPRES,
wellElemDofNumber[iwelem] + ColOffset::DPRES,
dP,
jacobianFD.toViewConstSizes() );
}
Expand Down Expand Up @@ -315,7 +316,7 @@ void testNumericalJacobian( SinglePhaseReservoirAndWells<> & solver,
// this is computing J_RW and J_WW
fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + singlePhaseWellKernels::ColOffset::DRATE,
wellElemDofNumber[iwelem] + ColOffset::DRATE,
dRate,
jacobianFD.toViewConstSizes() );
}
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -37,6 +37,7 @@ using namespace geos;
using namespace geos::dataRepository;
using namespace geos::constitutive;
using namespace geos::testing;
using namespace geos::wells::kernels::singlePhase;

CommandLineOptions g_commandLineOptions;

Expand Down Expand Up @@ -412,7 +413,7 @@ void testNumericalJacobian( SinglePhaseReservoirAndWells<> & solver,
// this is computing J_RW and J_WW
fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + singlePhaseWellKernels::ColOffset_WellJac< 1 >::dP,
wellElemDofNumber[iwelem] + ColOffset_WellJac< 1 >::dP,
dP,
jacobianFD.toViewConstSizes() );
}
Expand Down Expand Up @@ -440,13 +441,13 @@ void testNumericalJacobian( SinglePhaseReservoirAndWells<> & solver,
// this is computing J_RW and J_WW
fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + singlePhaseWellKernels::ColOffset_WellJac< 1 >::dT,
wellElemDofNumber[iwelem] + ColOffset_WellJac< 1 >::dT,
dT,
jacobianFD.toViewConstSizes() );

if( diag_check && !wellControls.isProducer() && iwelem == 0 )
{
localIndex rowIndex = wellElemDofNumber[iwelem] + singlePhaseWellKernels::ColOffset_WellJac< 1 >::dT;
localIndex rowIndex = wellElemDofNumber[iwelem] + ColOffset_WellJac< 1 >::dT;
setNumericalJacobianValue( rowIndex, rowIndex, 1.0, jacobianFD.toViewConstSizes() );
}
}
Expand All @@ -472,7 +473,7 @@ void testNumericalJacobian( SinglePhaseReservoirAndWells<> & solver,
// this is computing J_RW and J_WW
fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + singlePhaseWellKernels::ColOffset_WellJac< 1 >::dQ,
wellElemDofNumber[iwelem] + ColOffset_WellJac< 1 >::dQ,
dRate,
jacobianFD.toViewConstSizes() );
}
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -37,6 +37,7 @@ using namespace geos;
using namespace geos::dataRepository;
using namespace geos::constitutive;
using namespace geos::testing;
using namespace geos::wells::kernels::singlePhase;

CommandLineOptions g_commandLineOptions;

Expand Down Expand Up @@ -495,7 +496,7 @@ void testNumericalJacobian( SinglePhaseReservoirAndWells<> & solver,
// this is computing J_RW and J_WW
fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + singlePhaseWellKernels::ColOffset_WellJac< 1 >::dP,
wellElemDofNumber[iwelem] + ColOffset_WellJac< 1 >::dP,
dP,
jacobianFD.toViewConstSizes() );
}
Expand Down Expand Up @@ -523,13 +524,13 @@ void testNumericalJacobian( SinglePhaseReservoirAndWells<> & solver,
// this is computing J_RW and J_WW
fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + singlePhaseWellKernels::ColOffset_WellJac< 1 >::dT,
wellElemDofNumber[iwelem] + ColOffset_WellJac< 1 >::dT,
dT,
jacobianFD.toViewConstSizes() );

if( diag_check && !wellControls.isProducer() && iwelem == 0 )
{
localIndex rowIndex = wellElemDofNumber[iwelem] + singlePhaseWellKernels::ColOffset_WellJac< 1 >::dT;
localIndex rowIndex = wellElemDofNumber[iwelem] + ColOffset_WellJac< 1 >::dT;
setNumericalJacobianValue( rowIndex, rowIndex, 1.0, jacobianFD.toViewConstSizes() );
}
}
Expand All @@ -555,7 +556,7 @@ void testNumericalJacobian( SinglePhaseReservoirAndWells<> & solver,
// this is computing J_RW and J_WW
fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + singlePhaseWellKernels::ColOffset_WellJac< 1 >::dQ,
wellElemDofNumber[iwelem] + ColOffset_WellJac< 1 >::dQ,
dRate,
jacobianFD.toViewConstSizes() );
}
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -34,6 +34,7 @@ using namespace geos;
using namespace geos::dataRepository;
using namespace geos::constitutive;
using namespace geos::testing;
using namespace geos::wells::kernels::compositional;

CommandLineOptions g_commandLineOptions;

Expand Down Expand Up @@ -567,7 +568,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DPRES,
wellElemDofNumber[iwelem] + ColOffset::DPRES,
dP,
jacobianFD.toViewConstSizes() );
}
Expand All @@ -588,7 +589,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + jc,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + jc,
dRho,
jacobianFD.toViewConstSizes() );
}
Expand All @@ -610,13 +611,13 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio
assembleFunction( jacobian.toViewConstSizes(), residual.toView() );
fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC+1,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + NC+1,
dT,
jacobianFD.toViewConstSizes() );
if( iwelem == 1 )
{
real64 dRdX = 0.0;
localIndex rowIndex = wellElemDofNumber[0] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC+1;;
localIndex rowIndex = wellElemDofNumber[0] + ColOffset::DCOMP + NC+1;;
for( integer ider=0; ider< 3; ider++ )
{
globalIndex colIndex = wellElemDofNumber[0]+ ider;
Expand All @@ -628,8 +629,8 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio
}
else
{
localIndex rowIndex = wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC+1;;
globalIndex colIndex = wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC+1;;
localIndex rowIndex = wellElemDofNumber[iwelem] + ColOffset::DCOMP + NC+1;;
globalIndex colIndex = wellElemDofNumber[iwelem] + ColOffset::DCOMP + NC+1;;
setNumericalJacobianValue( rowIndex, colIndex, 1.0, jacobianFD.toViewConstSizes() );
}

Expand Down Expand Up @@ -657,7 +658,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + NC,
dRate,
jacobianFD.toViewConstSizes() );
}
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -34,6 +34,7 @@ using namespace geos;
using namespace geos::dataRepository;
using namespace geos::constitutive;
using namespace geos::testing;
using namespace geos::wells::kernels::compositional;

CommandLineOptions g_commandLineOptions;

Expand Down Expand Up @@ -569,7 +570,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DPRES,
wellElemDofNumber[iwelem] + ColOffset::DPRES,
dP,
jacobianFD.toViewConstSizes() );
}
Expand All @@ -590,7 +591,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + jc,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + jc,
dRho,
jacobianFD.toViewConstSizes() );
}
Expand All @@ -612,14 +613,14 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio
assembleFunction( jacobian.toViewConstSizes(), residual.toView() );
fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC+1,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + NC+1,
dT,
jacobianFD.toViewConstSizes() );
}
else
{
localIndex rowIndex = wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC+1;;
globalIndex colIndex = wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC+1;;
localIndex rowIndex = wellElemDofNumber[iwelem] + ColOffset::DCOMP + NC+1;;
globalIndex colIndex = wellElemDofNumber[iwelem] + ColOffset::DCOMP + NC+1;;
setNumericalJacobianValue( rowIndex, colIndex, 1.0, jacobianFD.toViewConstSizes() );
}
}
Expand All @@ -646,7 +647,7 @@ void testNumericalJacobian( CompositionalMultiphaseReservoirAndWells< Compositio

fillNumericalJacobian( residual.toViewConst(),
residualOrig.toViewConst(),
wellElemDofNumber[iwelem] + compositionalMultiphaseWellKernels::ColOffset::DCOMP + NC,
wellElemDofNumber[iwelem] + ColOffset::DCOMP + NC,
dRate,
jacobianFD.toViewConstSizes() );
}
Expand Down
10 changes: 10 additions & 0 deletions src/coreComponents/physicsSolvers/fluidFlow/CMakeLists.txt
Original file line number Diff line number Diff line change
Expand Up @@ -109,6 +109,16 @@ set( fluidFlowSolvers_headers
kernels/compositional/ThermalResidualNormKernel.hpp
kernels/compositional/ThermalSolutionCheckKernel.hpp
kernels/compositional/ThermalSolutionScalingKernel.hpp
kernels/compositional/zFormulation/AccumulationKernel.hpp
kernels/compositional/zFormulation/FluxComputeKernel.hpp
kernels/compositional/zFormulation/PhaseComponentFlux.hpp
kernels/compositional/zFormulation/PhaseVolumeFractionKernel.hpp
kernels/compositional/zFormulation/SolutionScalingAndCheckingKernelBase.hpp
kernels/compositional/zFormulation/DirichletFluxComputeKernel.hpp
kernels/compositional/zFormulation/PPUPhaseFlux.hpp
kernels/compositional/zFormulation/PhaseMobilityKernel.hpp
kernels/compositional/zFormulation/PotGrad.hpp
kernels/compositional/zFormulation/SolutionScalingKernel.hpp
kernels/immiscibleMultiphase/ImmiscibleMultiphaseKernels.hpp
wells/CompositionalMultiphaseWell.hpp
wells/CompositionalMultiphaseWellFields.hpp
Expand Down
Loading
Loading