[aspect-devel] Error from violated condition on parameters
Juliane Dannberg
dannberg at gfz-potsdam.de
Tue Nov 22 07:41:08 PST 2016
Hi Felipe,
if you run your input file in optimized mode you will see that the
density is indeed negative in a part of your model domain (probably
because your temperature goes up to almost 5000 K), and that is exactly
what the error message is telling you.
I would also like to add that normally it is not a good idea to just run
your models in optimized mode if you get an error message in debug mode,
because usually the error message tells you that something in your code
or your model setup is not working as intended, and in these cases it
helps to look at the setup again and fix these things.
Best,
Juliane
On 11/22/2016 01:56 AM, FELIPE ORELLANA ROVIROSA wrote:
>
> Thanks Mohamed,
>
> I will check this thing.
>
> tank you,
>
> Felipe
>
>
>
> On Mon, Nov 21, 2016 at 1:17 AM, Mohamed Gouiza <M.Gouiza at leeds.ac.uk
> <mailto:M.Gouiza at leeds.ac.uk>> wrote:
>
> Hi Felipe,
>
> Your prm file does run on my machine, see below.
>
> Either your ASPECT didn’t compile correctly (Try to run an example
> from the cookbook), or try in optimized mode rather than debug mode.
>
> Mohamed
>
> [earmgo at login1.polaris test]$ aspect RR3Dbic.prm
>
> -----------------------------------------------------------------------------
>
> -- This is ASPECT, the Advanced Solver for Problems in Earth's
> ConvecTion.
>
> -- . version 1.5.0-pre
>
> -- . running in OPTIMIZED mode
>
> -- . running with 1 MPI process
>
> -- . using Trilinos
>
> -----------------------------------------------------------------------------
>
> -----------------------------------------------------------------------------
>
> The output directory <output/> provided in the input file appears
> not to exist.
>
> ASPECT will create it for you.
>
> -----------------------------------------------------------------------------
>
> Number of active cells: 512 (on 4 levels)
>
> Number of degrees of freedom: 20,381 (14,739+729+4,913)
>
> *** Timestep 0: t=0 years
>
> Solving temperature system... 0 iterations.
>
> Rebuilding Stokes preconditioner...
>
> Solving Stokes system... 30+2 iterations.
>
> Number of active cells: 1,296 (on 5 levels)
>
> Number of degrees of freedom: 51,950 (37,617+1,794+12,539)
>
> *** Timestep 0: t=0 years
>
> Solving temperature system... 0 iterations.
>
> Rebuilding Stokes preconditioner...
>
> Solving Stokes system... 30+2 iterations.
>
> Number of active cells: 4,544 (on 6 levels)
>
> Number of degrees of freedom: 179,987 (130,407+6,111+43,469)
>
> *** Timestep 0: t=0 years
>
> Solving temperature system... 0 iterations.
>
> Rebuilding Stokes preconditioner...
>
> Solving Stokes system... 30+3 iterations.
>
> Number of active cells: 10,368 (on 6 levels)
>
> Number of degrees of freedom: 377,927 (274,239+12,275+91,413)
>
> *** Timestep 0: t=0 years
>
> Solving temperature system... 0 iterations.
>
> Rebuilding Stokes preconditioner...
>
> Solving Stokes system... 30+4 iterations.
>
> Postprocessing:
>
> RMS, max velocity: 0.00454 m/year, 0.0155 m/year
>
> Temperature min/avg/max: 300 K, 2.474e+05 K, 4.922e+05 K
>
> Heat fluxes through boundary parts: 0.001263 W, 0.0006153 W,
> 0.005619 W, 0.002371 W, 1.895e+14 W, 1.971e+12 W
>
> Writing graphical output: output/solution/solution-00000
>
> *** Timestep 1: t=1 years
>
> Solving temperature system... 3 iterations.
>
> Solving Stokes system... 8+0 iterations.
>
> Postprocessing:
>
> RMS, max velocity: 0.00454 m/year, 0.0155 m/year
>
> Temperature min/avg/max: 300 K, 2.474e+05 K, 4.922e+05 K
>
> Heat fluxes through boundary parts: 91.23 W, 113.4 W, 363.3 W,
> 454.7 W, 1.895e+14 W, 1.971e+12 W
>
> Termination requested by criterion: end time
>
> +---------------------------------------------+------------+------------+
>
> | Total wallclock time elapsed since start | 107s | |
>
> | | | |
>
> | Section | no. calls | wall time | % of total |
>
> +---------------------------------+-----------+------------+------------+
>
> | Assemble Stokes system | 5 | 33.6s | 31% |
>
> | Assemble temperature system | 5 | 7.01s | 6.5% |
>
> | Build Stokes preconditioner | 4 | 17.4s | 16% |
>
> | Build temperature preconditioner| 5 | 2.36s | 2.2% |
>
> | Solve Stokes system | 5 | 30.1s | 28% |
>
> | Solve temperature system | 5 | 0.158s | 0.15% |
>
> | Initialization | 1 | 0.123s | 0.11% |
>
> | Postprocessing | 2 | 1.6s | 1.5% |
>
> | Refine mesh structure, part 1 | 3 | 2.17s | 2% |
>
> | Refine mesh structure, part 2 | 3 | 0.451s | 0.42% |
>
> | Setup dof systems | 4 | 11.6s | 11% |
>
> | Setup initial conditions | 4 | 0.416s | 0.39% |
>
> +---------------------------------+-----------+------------+------------+
>
> *From:*Aspect-devel [mailto:aspect-devel-bounces at geodynamics.org
> <mailto:aspect-devel-bounces at geodynamics.org>] *On Behalf Of
> *FELIPE ORELLANA ROVIROSA
> *Sent:* 21 November 2016 04:35
> *To:* aspect-devel at geodynamics.org
> <mailto:aspect-devel at geodynamics.org>
> *Subject:* [aspect-devel] Error from violated condition on parameters
>
> Hi all,
>
> I am running aspect 1.5.0-pre
>
> I have tried to run my job several times, and I get this
> awkward error
>
> -- This is ASPECT, the Advanced Solver for Problems in Earth's
> ConvecTion.
> -- . version 1.5.0-pre
> -- . running in DEBUG mode
> -- . running with 32 MPI processes
> -- . using Trilinos
> -----------------------------------------------------------------------------
>
>
> -----------------------------------------------------------------------------
> The output directory </home1/04020/unfelipe/brickbic/> provided in
> the input file appears not to exist.
> ASPECT will create it for you.
> -----------------------------------------------------------------------------
>
>
> Number of active cells: 512 (on 4 levels)
> Number of degrees of freedom: 20,381 (14,739+729+4,913)
>
> *** Timestep 0: t=0 years
>
> --------------------------------------------------------
> An error occurred in line <1152> of file
> </home1/04020/unfelipe/packages/aspect/source/simulator/assembly.cc>
> in function
> void
> aspect::Assemblers::CompleteEquations<dim>::local_assemble_advection_system(const
> typename aspect::Simulator<dim>::Adv
> ectionField&, double,
> aspect::internal::Assembly::Scratch::AdvectionSystem<dim>&,
> aspect::internal::Assembly::CopyData::Advecti
> onSystem<dim>&) const [with int dim = 3; typename
> aspect::Simulator<dim>::AdvectionField =
> aspect::Simulator<3>::AdvectionField
> ]
> The violated condition was:
> density_c_P >= 0
> The name and call sequence of the exception was:
> ExcMessage ("The product of density and c_P needs to be a "
> "non-negative quantity.")
> Additional Information:
> The product of density and c_P needs to be a non-negative quantity.
>
> Stacktrace:
> -----------
> #0 /home1/04020/unfelipe/packages/aspect/build/aspect:
> aspect::Assemblers::CompleteEquations<3>::local_assemble_advection_syst
> em(aspect::Simulator<3>::AdvectionField const&, double,
> aspect::internal::Assembly::Scratch::AdvectionSystem<3>&,
> aspect::inter
> nal::Assembly::CopyData::AdvectionSystem<3>&) const
> #1 /home1/04020/unfelipe/packages/aspect/build/aspect: s...
> ..
>
> of course, my parameters are positive numbers, so I don't know
> where the error comes from.
>
> I am attaching my parameter file.
>
> any help would be appreciated,
>
> thanks,
>
> Felipe
>
>
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org <mailto:Aspect-devel at geodynamics.org>
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
> <http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel>
>
>
>
>
> _______________________________________________
> Aspect-devel mailing list
> Aspect-devel at geodynamics.org
> http://lists.geodynamics.org/cgi-bin/mailman/listinfo/aspect-devel
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161122/f978e403/attachment-0001.html>
More information about the Aspect-devel
mailing list