Hello, everyone,
I'm wondering if you have a incompressible N-S solver with adaptive
mesh refinement to share with me, it seems hard for me
to work out this issue, since I have struggled with it for several days.
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)
*******************************************************************************************
BTW, i attached the *.cc file to this email.
Any help is appreciated !
Regards,
zhicheng
--
Zhicheng Wang
Ph.D. Candidate
RM 402, No.11 Beisihuan West Road
Institute of Engineering Thermophysics, Chinese Academy of Sciences
Beijing, 100190, P.R. China
Tel: (86) 10-82543092
E-mail: [email protected]; [email protected]
AdaptiveNS.cc.tar.gz
Description: GNU Zip compressed data
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
