[aspect-devel] Compositional vector output / ASPECT 2.0

Diogo Louro Lourenco dlourenco at ucdavis.edu
Fri Jun 1 14:19:31 PDT 2018


Hi all,

When updating to ASPECT 2.0, and using compositional vector as an output
variable in section Postprocess/Visualization I get the following error,
which seems to be related with deal.II:

An error occurred in line <77> of file
</Users/ldiogo/Documents/aspect/source/postprocess/visualization/compositional_vector.cc>
in function
    void
aspect::Postprocess::VisualizationPostprocessors::CompositionalVector<dim>::evaluate_vector_field(const
dealii::DataPostprocessorInputs::Vector<dim>&,
std::vector<dealii::Vector<double> >&) const [with int dim = 3]
The violated condition was:
    computed_quantities.size() == n_quadrature_points
Additional information:
    This exception -- which is used in many places in the library --
usually indicates that some condition which the author of the code thought
must be satisfied at a certain point in an algorithm, is not fulfilled. An
example would be that the first part of an algorithm sorts elements of an
array in ascending order, and a second part of the algorithm later
encounters an element that is not larger than the previous one.

There is usually not very much you can do if you encounter such an
exception since it indicates an error in deal.II, not in your own program.
Try to come up with the smallest possible program that still demonstrates
the error and contact the deal.II mailing lists with it to obtain help.

When not using compositional vector as an output, ASPECT runs normally. I
am attaching a prm file, but it does not seem like the error is related
with the model input.

Cheers,
Diogo
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20180601/17b5e19f/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: par.prm
Type: application/octet-stream
Size: 2081 bytes
Desc: not available
URL: <http://lists.geodynamics.org/pipermail/aspect-devel/attachments/20180601/17b5e19f/attachment.obj>


More information about the Aspect-devel mailing list