Hi, The problem has been figured out. I should refine the mesh before starting running steps.
On Thu, Feb 16, 2012 at 3:54 PM, Zhicheng Wang <[email protected]>wrote: > Hi, Christian, > > Thank you so much for your reply. > actually, I did close the 'close' the ConstraintMatrix in stepup_dofs(): > > > ****************************************************************************** > velocity_constraints.clear (); > DoFTools::make_hanging_node_constraints (dof_handler_velocity, > velocity_constraints); > velocity_constraints.close (); > > ************************************************************************************ > > Regards, > > Zhicheng > > On Thu, Feb 16, 2012 at 3:24 PM, Christian Goll < > [email protected]> wrote: > >> 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<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<http://poisson.dealii.org/mailman/listinfo/dealii> >> > > > > -- > 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] > -- 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]
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
