[aspect-devel] Error from violated condition on parameters

Mon Nov 21 01:17:49 PST 2016

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20161121/fa5dd304/attachment-0001.html>