Dear Chucui, Thanks for giving this a go. This is better than before, but I think that you still misinterpret what a minimal working example is. This is a 1000 line piece of code, but I estimate only 100 lines of code are required to reproduce the problem. If the problem arises in the first call to setup_dofs(), which to the best of my understanding it does, then you really don’t need anything other than - the set of calls used to create the Triangulation, and - the call to setup_dofs() itself.
So could you please remove everything below the call to setup_dofs() in run(), and everything related to the assembly, solving the linear system, exact solution, etc. Forget about solving the original problem -- really try to make this compact while reproducing the error message hit when creating the sparse matrix. Then we will have a focussed piece of code that we can use to debug the issue, and if there is something to be changed in the library we will add this code to the test suite and would naturally give you attribution for it. Best, Jean-Paul P.S. A huge plus would be to strip all of the remaining functional code out of the Problem class and have it all in the run() function. This would be nice, but is not entirely necessary. > On 08 Aug 2018, at 15:40, [email protected] wrote: > > Dear Daniel, > > I'm sorry to make this mistake, (and thank Jean-Paul very much). The newest > code focus on build the sparsity pattern in "setup_dof()", while other parts > (assemble_system(), assemble_rhs(), run() ) need to be left to make the code > work and maintain the structure of 3 components of the blockvector solution. > > Thank you much! > > Best, > Chucui > > -- > The deal.II project is located at http://www.dealii.org/ > <http://www.dealii.org/> > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > <https://groups.google.com/d/forum/dealii?hl=en> > --- > You received this message because you are subscribed to the Google Groups > "deal.II User Group" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to [email protected] > <mailto:[email protected]>. > For more options, visit https://groups.google.com/d/optout > <https://groups.google.com/d/optout>. > <LittleCase0-3.cc> -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. For more options, visit https://groups.google.com/d/optout.
