Hi Kai, are you sure this is within the try & catch block of the Newton? It be very surprised if the try/catch does not catch the exception. You can use a debugger to find out where this exactly happens. The time step reduction mechanism is only activated if you call newton.solve(x, timeLoop) with the time loop. Are you doing that?
Timo > On 26 Aug 2022, at 15:53, Kai Wendel <kaiwende...@googlemail.com> wrote: > > Hello, > > by trying to implement some experimental features, I did now for the second > time encounter the following problem: > A porous medium flow simulations runs until an error occurs, that should > cause the newton solver to retry e.g. with a smaller timestep. The initial > solution therefore has already successfully been built. > I don't understand, why the exception is not catched correclty and instead of > doing the call in the newtonsolver, that catches this exception, the programm > is terminated in a very harsh way. > The error message that I get from the compiled programs reads like the > following example: > > terminate called after throwing an instance of 'Dumux::NumericalProblem' > what(): NumericalProblem > [oilSaturationPressure:/home/kai/PaperWork/Codes/SPE-tests/spe-tests/modulerelated/material/fluidsystems/lauserstandardblackoilfluidsystem.hh:659]: > Could not find the oil saturation pressure for X_o^g = -0.0815226 > [computer:12345] *** Process received signal *** > [computer:12345] Signal: Aborted (6) > [computer:12345] Signal code: (-6) > [computer:12345] [ 0] > /lib/x86_64-linux-gnu/libpthread.so.0(+0x14420)[0x7f583c456420] > [computer:12345] [ 1] > /lib/x86_64-linux-gnu/libc.so.6(gsignal+0xcb)[0x7f583bd7300b] > [computer:12345] [ 2] > /lib/x86_64-linux-gnu/libc.so.6(abort+0x12b)[0x7f583bd52859] > [computer:12345] [ 3] > /lib/x86_64-linux-gnu/libstdc++.so.6(+0x9e911)[0x7f583c12c911] > [computer:12345] [ 4] > /lib/x86_64-linux-gnu/libstdc++.so.6(+0xaa38c)[0x7f583c13838c] > [computer:12345] [ 5] > /lib/x86_64-linux-gnu/libstdc++.so.6(+0xaa3f7)[0x7f583c1383f7] > [computer:12345] [ 6] > /lib/x86_64-linux-gnu/libstdc++.so.6(+0xaa6a9)[0x7f583c1386a9] > [computer:12345] [ 7] > ./SPE2_box(_ZN5Dumux12FluidSystems8BlackOilIdE21oilSaturationPressureEd+0x39b)[0x55ad9a96261b] > [computer:12345] [ 8] ./SPE2_box(+0x1d3fb8)[0x55ad9a91afb8] > [computer:12345] [ 9] > ./SPE2_box(_ZN5Dumux23BlackOilVolumeVariablesINS_10Properties15VolumeVariablesINS1_4TTag18OilyCakeBoxTypeTagENS3_8BlackOilEE8NCTraitsINS_27ThreePVolumeVariablesTraitsIN4Dune11FieldVectorIdLi3EEENS_12FluidSystems8BlackOilIdEENS_23CompositionalFluidStateIdSE_EENS_12SolidSystems9OneCSolidIdNS_10Components8ConstantILi1EdEELb1EEENS_15InertSolidStateIdSM_EENS9_11FieldMatrixIdLi3ELi3EEENS_19BlackOilModelTraitsILb0ELb0EEEEENS_22FicksLawImplementationIS4_NS_21DiscretizationMethods4CVFEINSV_11CVFEMethods3PQ1EEELNS_26ReferenceSystemFormulationE0EEENS_26DiffusivityMillingtonQuirkIdEEEEE6updateINS_18BoxElementSolutionINS_20BoxFVElementGeometryINS_17BoxFVGridGeometryIdNS9_8GridViewINS9_23ALU3dLeafGridViewTraitsIKNS9_7ALUGridILi3ELi3ELNS9_18ALUGridElementTypeE1ELNS9_21ALUGridRefinementTypeE1ENS9_14ALUGridMPICommEEELNS9_21PartitionIteratorTypeE4EEEEELb0ENS_28BoxDefaultGridGeometryTraitsIS1K_NS_19DefaultMapperTraitsIS1K_NS9_35MultipleCodimMultipleGeomTypeMapperIS1K_EES1O_EEEEE ELb0EEESB_EENS_15OilyCakeProblemIS4_EENS9_6EntityILi0ELi3EKNS9_9ALU3dGridILi3ELi3ELNS9_20ALU3dGridElementTypeE7ES1F_EENS9_15ALU3dGridEntityEEENS_19BoxSubControlVolumeIS1K_NS_27BoxDefaultScvGeometryTraitsIS1K_EEEEEEvRKT_RKT0_RKT1_RKT2_+0x50d)[0x55ad9a9bd81d] > [computer:12345] [10] > ./SPE2_box(_ZNR5Dumux26CVFEElementVolumeVariablesINS_23CVFEGridVolumeVariablesINS_36CVFEDefaultGridVolumeVariablesTraitsINS_15OilyCakeProblemINS_10Properties4TTag18OilyCakeBoxTypeTagEEENS_23BlackOilVolumeVariablesINS4_15VolumeVariablesIS6_NS5_8BlackOilEE8NCTraitsINS_27ThreePVolumeVariablesTraitsIN4Dune11FieldVectorIdLi3EEENS_12FluidSystems8BlackOilIdEENS_23CompositionalFluidStateIdSJ_EENS_12SolidSystems9OneCSolidIdNS_10Components8ConstantILi1EdEELb1EEENS_15InertSolidStateIdSR_EENSE_11FieldMatrixIdLi3ELi3EEENS_19BlackOilModelTraitsILb0ELb0EEEEENS_22FicksLawImplementationIS6_NS_21DiscretizationMethods4CVFEINS10_11CVFEMethods3PQ1EEELNS_26ReferenceSystemFormulationE0EEENS_26DiffusivityMillingtonQuirkIdEEEEEEEELb0EEELb0EE11bindElementINS_20BoxFVElementGeometryINS_17BoxFVGridGeometryIdNSE_8GridViewINSE_23ALU3dLeafGridViewTraitsIKNSE_7ALUGridILi3ELi3ELNSE_18ALUGridElementTypeE1ELNSE_21ALUGridRefinementTypeE1ENSE_14ALUGridMPICommEEELNSE_21PartitionIteratorTypeE4EEEEEL b0ENS_28BoxDefaultGridGeometryTraitsIS1R_NS_19DefaultMapperTraitsIS1R_NSE_35MultipleCodimMultipleGeomTypeMapperIS1R_EES1V_EEEEEELb0EEENSE_11BlockVectorISG_SaISG_EEEEEvRKNT_12GridGeometry8GridView5CodimILi0EE6EntityERKS23_RKT0_+0x353)[0x55ad9a9bf213] > [computer:12345] [11] > ./SPE2_box(_ZN5Dumux21BoxLocalAssemblerBaseINS_10Properties4TTag18OilyCakeBoxTypeTagENS_11FVAssemblerIS3_LNS_10DiffMethodE0ELb1EEENS_17BoxLocalAssemblerIS3_S6_LS5_0ELb1EEELb1EE27assembleJacobianAndResidualINS_18PartialReassemblerIS6_EEEEvRN4Dune10BCRSMatrixINSD_11FieldMatrixIdLi3ELi3EEESaISG_EEERNSD_11BlockVectorINSD_11FieldVectorIdLi3EEESaISM_EEERNS_15FVGridVariablesINS_17BoxFVGridGeometryIdNSD_8GridViewINSD_23ALU3dLeafGridViewTraitsIKNSD_7ALUGridILi3ELi3ELNSD_18ALUGridElementTypeE1ELNSD_21ALUGridRefinementTypeE1ENSD_14ALUGridMPICommEEELNSD_21PartitionIteratorTypeE4EEEEELb0ENS_28BoxDefaultGridGeometryTraitsIS12_NS_19DefaultMapperTraitsIS12_NSD_35MultipleCodimMultipleGeomTypeMapperIS12_EES16_EEEEEENS_23CVFEGridVolumeVariablesINS_36CVFEDefaultGridVolumeVariablesTraitsINS_15OilyCakeProblemIS3_EENS_23BlackOilVolumeVariablesINS1_15VolumeVariablesIS3_NS2_8BlackOilEE8NCTraitsINS_27ThreePVolumeVariablesTraitsISM_NS_12FluidSystems8BlackOilIdEENS_23CompositionalFluidSt ateIdS1M_EENS_12SolidSystems9OneCSolidIdNS_10Components8ConstantILi1EdEELb1EEENS_15InertSolidStateIdS1U_EESG_NS_19BlackOilModelTraitsILb0ELb0EEEEENS_22FicksLawImplementationIS3_NS_21DiscretizationMethods4CVFEINS21_11CVFEMethods3PQ1EEELNS_26ReferenceSystemFormulationE0EEENS_26DiffusivityMillingtonQuirkIdEEEEEEEELb0EEENS_26CVFEGridFluxVariablesCacheIS1D_NS_44PorousMediumFluxVariablesCacheImplementationIS3_S25_EELb0ENS_24CVFEDefaultGridFVCTraitsIS1D_S2G_EEEEEEPKT_+0x96)[0x55ad9a9e4b86] > [computer:12345] [12] > ./SPE2_box(_ZZN5Dumux11FVAssemblerINS_10Properties4TTag18OilyCakeBoxTypeTagELNS_10DiffMethodE0ELb1EE27assembleJacobianAndResidualINS_18PartialReassemblerIS5_EEEEvRKN4Dune11BlockVectorINS9_11FieldVectorIdLi3EEESaISC_EEEPKT_ENKUlRKNS9_6EntityILi0ELi3EKNS9_9ALU3dGridILi3ELi3ELNS9_20ALU3dGridElementTypeE7ENS9_14ALUGridMPICommEEENS9_15ALU3dGridEntityEEEE_clEST_+0x59)[0x55ad9a9e6279] > [computer:12345] [13] ./SPE2_box(+0x29f4b7)[0x55ad9a9e64b7] > [computer:12345] [14] > /lib/x86_64-linux-gnu/libgomp.so.1(+0x1a78e)[0x7f583c47f78e] > [computer:12345] [15] > /lib/x86_64-linux-gnu/libpthread.so.0(+0x8609)[0x7f583c44a609] > [computer:12345] [16] > /lib/x86_64-linux-gnu/libc.so.6(clone+0x43)[0x7f583be4f133] > [computer:12345] *** End of error message *** > Abgebrochen (Speicherabzug geschrieben) > > Thanks for your help! > > Best > Kai Wendel > _______________________________________________ > DuMux mailing list > DuMux@listserv.uni-stuttgart.de > https://listserv.uni-stuttgart.de/mailman/listinfo/dumux _______________________________________________ DuMux mailing list DuMux@listserv.uni-stuttgart.de https://listserv.uni-stuttgart.de/mailman/listinfo/dumux