[aspect-devel] Error from violated condition on parameters
FELIPE ORELLANA ROVIROSA
f_orellana at berkeley.edu
Mon Nov 21 22:56:55 PST 2016
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>
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] *On
> Behalf Of *FELIPE ORELLANA ROVIROSA
> *Sent:* 21 November 2016 04:35
> *To:* 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
> 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/20161121/9d8f7053/attachment-0001.html>
More information about the Aspect-devel
mailing list