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]
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
