Hello,
I can run successfuly the testcases proposed in implicit/2p, but those from
multidomain directory do not run. The error given by the multidomain run is in
the attached file. I can't trace this error message, but I guess that it comes
from Boost::FUSION library, which is required explicitly by multidomain module.
Regarding the Boost::Fusion module, I notices that the Boost flags are
different in the files
dune-common/cmake/modules/FindBoostFusion.cmake uses the flag BOOST_FOUND
dune-common/cmake/modules/DuneBoost.cmake uses the flag Boost_FOUND
Fixing this difference, the compiler finds the Fusion library
-- Boost version: 1.57.0
-- Checking whether the Boost::FUSION library is available.
-- Performing Test HAVE_BOOST_FUSION
-- Performing Test HAVE_BOOST_FUSION - Success
-- Boost::FUSION is available
but BoostFusion is still in the list of "The following OPTIONAL packages have
not been found:". Can be this the reason for the failure of multidomain
testcase?
Thank you in advance,
Calin
-----Original Message-----
From: Dan, Calin-Marius
Sent: Freitag, 17. April 2015 16:53
To: DuMuX User Forum
Subject: RE: [DuMuX] dunecontrol options file - example
Hello Christoph,
Thank you for explanations, it compiles now! I have another problem with
Boost::Fusion needed by dune-multidomaingrid. I have installed boost_1.5[7,8],
both contain the Fusion libraries in boost_1_5?_0/include/boost/fusion During
the compilation, it finds boost, but it does not check for Fusion:
-- Boost version: 1.58.0
-- Skipping check for Boost::FUSION as Boost is not available.
I attached the compilation outputs. Do you have any idea what could be wrong?
regards,
calin
-----Original Message-----
From: Dumux [mailto:[email protected]] On Behalf Of
Christoph Grüninger
Sent: Freitag, 17. April 2015 11:27
To: [email protected]
Subject: Re: [DuMuX] dunecontrol options file - example
Hi Calin,
I assume you use DuMuX 2.7 where we provide opts files for CMake.
Unfortunately I forgot to add "USE_CMAKE=1" in them.
1. You can use the CMAKE_PREFIX_PATH as described in [1].
2. Alternatively you can set -DSUPERLU_ROOT and -DUG_ROOT.
Bye
Christoph
[1]
http://users.dune-project.org/projects/main-wiki/wiki/Using_cmake_to_build_your_DUNE_project
--
People who are more than casually interested in computers should have at least
some idea of what the underlying hardware is like. Otherwise
the programs they write will be pretty weird. -- Donald Knuth
No parameter file given.
Defaulting to './test_2cstokes2p2c.input' for input file.
Starting with timestep size = 0.5s, simulation end = 7200s
YaspGridParameterBlock: Parameter 'overlap' not specified, defaulting to '0'.
rglobal = 0 0.25, rlocal = 0 0.5
lengthLeft = 0, lengthRight = 1, hLeft = 0.0333333, hRight = 0.0333333, nLeft =
0, nRight = 30
lengthLeft = 0.5, lengthRight = 0.5, hLeft = 0.0015046, hRight = 0.0015046,
nLeft = 31, nRight = 31
ParameterNewtonRelTol= , 1e-05
=>
Dune::mdgrid::subdomain::SubDomainGrid<Dune::mdgrid::MultiDomainGrid<Dune::UGGrid<2>,
Dune::mdgrid::FewSubDomainsTraits<2, 4ul, Dune::mdgrid::AllCodims> > > (dim=2,
dimworld=2)
level 0 codim[0]=930 codim[1]=1921 codim[2]=992
leaf codim[0]=930 codim[1]=1921 codim[2]=992
leaf dim=2 geomTypes=((cube, 2)[0]=930,(simplex, 1)[1]=1921,(simplex, 0)[2]=992)
=>
Dune::mdgrid::subdomain::SubDomainGrid<Dune::mdgrid::MultiDomainGrid<Dune::UGGrid<2>,
Dune::mdgrid::FewSubDomainsTraits<2, 4ul, Dune::mdgrid::AllCodims> > > (dim=2,
dimworld=2)
level 0 codim[0]=930 codim[1]=1921 codim[2]=992
leaf codim[0]=930 codim[1]=1921 codim[2]=992
leaf dim=2 geomTypes=((cube, 2)[0]=930,(simplex, 1)[1]=1921,(simplex, 0)[2]=992)
Using fast H2O-Air fluid system
# Run-time specified parameters:
[ FreeFlow ]
RefMassfrac = "0.008"
RefPressure = "1e5"
RefTemperature = "298.15"
RefVelocity = "3.5"
SinusConcentrationAmplitude = "0.0"
SinusConcentrationPeriod = "3600"
SinusPressureAmplitude = "0.0"
SinusPressurePeriod = "3600"
SinusTemperatureAmplitude = "0.0"
SinusTemperaturePeriod = "3600"
SinusVelAmplitude = "0.0"
SinusVelPeriod = "3600"
[ Grid ]
CellsX = "30"
CellsY = "62"
File = "grids/interfacedomain.dgf"
GradingFactorY = "1.13"
InterfacePosY = "0.25"
NoDarcyX = "0.0"
RunUpDistanceX = "0.0"
UseInterfaceMeshCreator = "true"
XMax = "0.25"
XMin = "0.0"
YMax = "0.5"
YMin = "0.0"
[ Newton ]
MaxRelativeShift = "1e-5"
MaxSteps = "12"
TargetSteps = "8"
[ Output ]
FreqFluxOutput = "1000"
FreqMassOutput = "2"
FreqOutput = "50"
FreqRestart = "1000"
FreqVaporFluxOutput = "2"
[ PorousMedium ]
InitialSw = "0.98"
RefPressure = "1e5"
RefTemperature = "298.15"
[ SpatialParams ]
AlphaBJ = "1.0"
LambdaSolid = "5.3"
Permeability = "2.65e-10"
Porosity = "0.41"
Snr = "0.01"
Swr = "0.005"
VgAlpha = "6.371e-4"
VgN = "6.9"
[ Stokes ]
StabilizationAlpha = "-1.0"
[ TimeManager ]
DtInitial = "5e-1"
EpisodeLength = "3600"
InitTime = "0"
TEnd = "7200"
# Compile-time specified parameters:
[ Implicit ]
EnableHints = "0"
EnableJacobianRecycling = "0"
EnablePartialReassemble = "0"
MassUpwindWeight = "1"
NumericDifferenceMethod = "1"
[ Newton ]
EnableResidualCriterion = "0"
EnableShiftCriterion = "1"
ResidualReduction = "1e-05"
SatisfyResidualAndShiftCriterion = "0"
UseLineSearch = "0"
[ Problem ]
EnableGravity = "1"
[ Stokes ]
StabilizationBeta = "0"
# UNUSED parameters:
TimeManager.MaxTimeStepSize = "360"
Grid.Refinement = "0"
Output.NameFF = "stokes2c"
Output.NamePM = "2p2c"
FreeFlow.ExponentMTC = "0.0"
FreeFlow.BoundaryLayerOffset = "0.0"
FreeFlow.ConstThickness = "0.0016"
PorousMedium.CharPoreDiameter = "1e-4"
Newton.WriteConvergence = "false"
Newton.MaxTimeStepDivisions = "20"
LinearSolver.ResidualReduction = "1e-10"
LinearSolver.Verbosity = "0"
LinearSolver.MaxIterations = "200"
Initializing problem 'simCoupled'
Writing flux data at interface
Writing result file for "stokes2c"
Writing result file for "2p2c"
Initialization took 1.38004 seconds on 1 processes.
The cumulative CPU time was 1.38004 seconds.
Properties for TwoCStokesTwoPTwoCProblem:
opaque FluidSystem defined at
/home/cdan/dev/dune/dumux/test/multidomain/2cstokes2p2c/2cstokes2p2cproblem.hh:95
type Grid = 'Dune::UGGrid<2>' defined at
/home/cdan/dev/dune/dumux/test/multidomain/2cstokes2p2c/2cstokes2p2cproblem.hh:64
type LinearSolver = 'SuperLUBackend<TypeTag>' defined at
/home/cdan/dev/dune/dumux/test/multidomain/2cstokes2p2c/2cstokes2p2cproblem.hh:107
type MultiDomainCouplingLocalOperator =
'Dumux::TwoCStokesTwoPTwoCLocalOperator<TypeTag>' defined at
/home/cdan/dev/dune/dumux/test/multidomain/2cstokes2p2c/2cstokes2p2cproblem.hh:76
type Problem = 'Dumux::TwoCStokesTwoPTwoCProblem<TypeTag>' defined at
/home/cdan/dev/dune/dumux/test/multidomain/2cstokes2p2c/2cstokes2p2cproblem.hh:72
type SubDomain1TypeTag = '::Dumux::Properties::TTag::Stokes2cSubProblem'
defined at
/home/cdan/dev/dune/dumux/test/multidomain/2cstokes2p2c/2cstokes2p2cproblem.hh:79
type SubDomain2TypeTag = '::Dumux::Properties::TTag::TwoPTwoCSubProblem'
defined at
/home/cdan/dev/dune/dumux/test/multidomain/2cstokes2p2c/2cstokes2p2cproblem.hh:80
Inherited from TwoCStokesTwoPTwoC:
int FreeFlowBoundaryLayerModel = '0' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/2cstokes2p2c/2cstokes2p2cpropertydefaults.hh:72
int FreeFlowMassTransferModel = '0' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/2cstokes2p2c/2cstokes2p2cpropertydefaults.hh:75
opaque GridView defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/2cstokes2p2c/2cstokes2p2cpropertydefaults.hh:43
type JacobianAssembler = 'Dumux::MultiDomainAssembler<TypeTag>' defined
at
/home/cdan/dev/dune/dumux/dumux/multidomain/2cstokes2p2c/2cstokes2p2cpropertydefaults.hh:60
type NewtonController =
'Dumux::TwoCStokesTwoPTwoCNewtonController<TypeTag>' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/2cstokes2p2c/2cstokes2p2cpropertydefaults.hh:63
bool NewtonWriteConvergence = 'false' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/2cstokes2p2c/2cstokes2p2cpropertydefaults.hh:69
int NumEq = '1' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/2cstokes2p2c/2cstokes2p2cpropertydefaults.hh:66
opaque SolutionVector defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/2cstokes2p2c/2cstokes2p2cpropertydefaults.hh:51
Inherited from MultiDomain:
opaque JacobianMatrix defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:170
int LinearSolverBlockSize = '::Dumux::Properties::GetProperty<TypeTag,
::Dumux::Properties::PTag::NumEq>::p::value' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:178
int LinearSolverMaxIterations = '250' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:238
scalar LinearSolverResidualReduction = '1e-6' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:235
type Model = 'MultiDomainModel<TypeTag>' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:183
opaque MultiDomainCondition defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:89
opaque MultiDomainConstraintsTrafo defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:143
opaque MultiDomainCoupling defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:129
opaque MultiDomainGrid defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:62
opaque MultiDomainGridFunctionSpace defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:71
opaque MultiDomainGridOperator defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:152
opaque MultiDomainSubProblem1 defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:97
opaque MultiDomainSubProblem2 defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:113
scalar NewtonMaxRelativeShift = '1e-5' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainnewtoncontroller.hh:43
int NewtonMaxSteps = '15' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainnewtoncontroller.hh:42
int NewtonMaxTimeStepDivisions = '10' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:241
type NewtonMethod = 'NewtonMethod<TypeTag>' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:194
int NewtonTargetSteps = '8' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainnewtoncontroller.hh:41
opaque NumEq1 defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:215
opaque NumEq2 defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:223
type SplitAndMerge = 'SplitAndMerge<TypeTag>' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:244
type TimeManager = 'TimeManager<TypeTag>' defined at
/home/cdan/dev/dune/dumux/dumux/multidomain/common/multidomainpropertydefaults.hh:200
Inherited from ImplicitBase:
type BaseModel = 'Dumux::ImplicitModel<TypeTag>' defined at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:80
type BoundaryTypes =
'Dumux::BoundaryTypes<::Dumux::Properties::GetProperty<TypeTag,
::Dumux::Properties::PTag::NumEq>::p::value>' defined at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:107
type ElementMapper =
'Dune::MultipleCodimMultipleGeomTypeMapper<typename
::Dumux::Properties::GetProperty<TypeTag,
::Dumux::Properties::PTag::GridView>::p::type, Dune::MCMGElementLayout>'
defined at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:77
type ElementSolutionVector = 'Dune::BlockVector<typename
::Dumux::Properties::GetProperty<TypeTag,
::Dumux::Properties::PTag::PrimaryVariables>::p::type>' defined at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:96
bool ImplicitEnableHints = 'false' defined at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:113
bool ImplicitEnableJacobianRecycling = 'false' defined at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:116
bool ImplicitEnablePartialReassemble = 'false' defined at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:119
int ImplicitMaxTimeStepDivisions = '10' defined at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:149
int ImplicitNumericDifferenceMethod = '+1' defined at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:110
type LocalJacobian = 'Dumux::ImplicitLocalJacobian<TypeTag>' defined
at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:86
type PrimaryVariables = 'Dune::FieldVector<typename
::Dumux::Properties::GetProperty<TypeTag,
::Dumux::Properties::PTag::Scalar>::p::type,
::Dumux::Properties::GetProperty<TypeTag,
::Dumux::Properties::PTag::NumEq>::p::value>' defined at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:102
type VertexMapper =
'Dune::MultipleCodimMultipleGeomTypeMapper<typename
::Dumux::Properties::GetProperty<TypeTag,
::Dumux::Properties::PTag::GridView>::p::type, Dune::MCMGVertexLayout>' defined
at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:71
type VolumeVariables = 'Dumux::ImplicitVolumeVariables<TypeTag>'
defined at
/home/cdan/dev/dune/dumux/dumux/implicit/common/implicitpropertydefaults.hh:83
Inherited from ImplicitModel: (none)
Inherited from NumericModel:
opaque Components defined at
/home/cdan/dev/dune/dumux/dumux/material/fluidsystems/defaultcomponents.hh:85
opaque DefaultComponents defined at
/home/cdan/dev/dune/dumux/dumux/material/fluidsystems/defaultcomponents.hh:58
opaque EnableComplicatedFluidSystem defined at
/home/cdan/dev/dune/dumux/dumux/material/fluidsystems/defaultcomponents.hh:100
type GridCreator = 'Dumux::DgfGridCreator<TypeTag>' defined at
/home/cdan/dev/dune/dumux/dumux/common/basicproperties.hh:155
string ModelParameterGroup = '""' defined at
/home/cdan/dev/dune/dumux/dumux/common/basicproperties.hh:152
opaque ParameterTree defined at
/home/cdan/dev/dune/dumux/dumux/common/basicproperties.hh:116
type Scalar = 'double' defined at
/home/cdan/dev/dune/dumux/dumux/common/basicproperties.hh:110
scalar TimeManagerMaxTimeStepSize = 'std::numeric_limits<typename
::Dumux::Properties::GetProperty<TypeTag,
::Dumux::Properties::PTag::Scalar>::p::type>::max()' defined at
/home/cdan/dev/dune/dumux/dumux/common/basicproperties.hh:113
opaque VtkMultiWriter defined at
/home/cdan/dev/dune/dumux/dumux/common/basicproperties.hh:161
int VtkOutputLevel = '0' defined at
/home/cdan/dev/dune/dumux/dumux/common/basicproperties.hh:158
Inherited from LinearSolverTypeTag:
int LinearSolverGMResRestart = '10' defined at
/home/cdan/dev/dune/dumux/dumux/linear/linearsolverproperties.hh:84
int LinearSolverPreconditionerIterations = '1' defined at
/home/cdan/dev/dune/dumux/dumux/linear/linearsolverproperties.hh:81
scalar LinearSolverPreconditionerRelaxation = '1.0' defined at
/home/cdan/dev/dune/dumux/dumux/linear/linearsolverproperties.hh:78
int LinearSolverVerbosity = '0' defined at
/home/cdan/dev/dune/dumux/dumux/linear/linearsolverproperties.hh:75
Inherited from NewtonMethod:
scalar NewtonAbsTolerance = '1e-5' defined at
/home/cdan/dev/dune/dumux/dumux/nonlinear/newtoncontroller.hh:146
bool NewtonEnableAbsoluteCriterion = 'false' defined at
/home/cdan/dev/dune/dumux/dumux/nonlinear/newtoncontroller.hh:131
bool NewtonEnableRelativeCriterion = 'true' defined at
/home/cdan/dev/dune/dumux/dumux/nonlinear/newtoncontroller.hh:126
bool NewtonEnableResidualCriterion =
'::Dumux::Properties::GetProperty<TypeTag,
::Dumux::Properties::PTag::NewtonEnableAbsoluteCriterion>::p::value' defined at
/home/cdan/dev/dune/dumux/dumux/nonlinear/newtoncontroller.hh:130
bool NewtonEnableShiftCriterion =
'::Dumux::Properties::GetProperty<TypeTag,
::Dumux::Properties::PTag::NewtonEnableRelativeCriterion>::p::value' defined at
/home/cdan/dev/dune/dumux/dumux/nonlinear/newtoncontroller.hh:125
scalar NewtonRelTolerance = '1e-8' defined at
/home/cdan/dev/dune/dumux/dumux/nonlinear/newtoncontroller.hh:141
scalar NewtonResidualReduction =
'::Dumux::Properties::GetProperty<TypeTag,
::Dumux::Properties::PTag::NewtonAbsTolerance>::p::value' defined at
/home/cdan/dev/dune/dumux/dumux/nonlinear/newtoncontroller.hh:145
bool NewtonSatisfyAbsAndRel = 'false' defined at
/home/cdan/dev/dune/dumux/dumux/nonlinear/newtoncontroller.hh:136
bool NewtonSatisfyResidualAndShiftCriterion =
'::Dumux::Properties::GetProperty<TypeTag,
::Dumux::Properties::PTag::NewtonSatisfyAbsAndRel>::p::value' defined at
/home/cdan/dev/dune/dumux/dumux/nonlinear/newtoncontroller.hh:135
bool NewtonUseLineSearch = 'false' defined at
/home/cdan/dev/dune/dumux/dumux/nonlinear/newtoncontroller.hh:121
Inherited from LinearSolverTypeTag: (none)
Inherited from NewtonMethod: (none)
Assemble: r(x^k) = dS/dt + div F - q; M = grad rtest_2cstokes2p2c:
/home/cdan/dev/dune/dune-multidomaingrid/dune/grid/multidomaingrid/indexsets.hh:502:
Dune::mdgrid::IndexSetWrapper<GridImp, HostGridViewType>::IndexType
Dune::mdgrid::IndexSetWrapper<GridImp,
HostGridViewType>::getSubIndexForSubDomain::invoke() const [with int codim = 2;
GridImp = const Dune::mdgrid::MultiDomainGrid<Dune::UGGrid<2>,
Dune::mdgrid::FewSubDomainsTraits<2, 4ul, Dune::mdgrid::AllCodims> >;
HostGridViewType = Dune::GridView<Dune::DefaultLeafGridViewTraits<const
Dune::UGGrid<2>, (Dune::PartitionIteratorType)4u> >;
Dune::mdgrid::IndexSetWrapper<GridImp, HostGridViewType>::IndexType = unsigned
int]: Assertion `me.domains.contains(_subDomain)' failed.
_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
