Am 16.02.2012 07:59, schrieb Zhicheng Wang:
Basically, what I have did is add two functions to tutorial
step-35 : refine_mesh(), which refines the mesh and transfer the
solution; setup_dofs(), which reset the matrices. However, when I
try to run the program, it always give me errors as follows:
--------------------------------------------------------
An error occurred in line <346> of file
</home/zwang/Academy/Codes/deal.II/include/deal.II/lac/constraint_matrix.templates.h>
in function
void
dealii::ConstraintMatrix::condense(dealii::SparseMatrix<number>&,
VectorType&) const [with number = double, VectorType =
dealii::Vector<double>]
The violated condition was:
sorted == true
The name and call sequence of the exception was:
ExcMatrixNotClosed()
Additional Information:
(none)
Dear Zhicheng,
after you have added all your constraints into the ConstraintMatrix, you
have to 'close' the matrix, see the documentation:
http://www.dealii.org/developer/doxygen/deal.II/classConstraintMatrix.html
I did not look at your code, but his should solve the problem.
Best regards,
Christian
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
