Hi Ibrahim,
My guess is that the mesh has been exported without the connectivity between
the parts being enforced. So even though the physical vertices of the mesh
overlap where you have different parts, there are two or more distinct vertices
that lie on top of one another and the
Yeah! .block(0) is essential. Without this, the data obtained is garbage! Thanks for this great hint! From: Wolfgang BangerthSent: Wednesday, August 18, 2021 8:37 PMTo: dealii@googlegroups.comSubject: Re: [deal.II] How to watch the variable of BlockVector type? On 8/18/21 7:18 PM, Michael Lee
On 8/18/21 9:20 AM, ibrahim zawra wrote:
Thanks a lot that helped, I applied symmetric<2,dim> tensor because it is a
heat problem. I also implemented the convection boundary condition. I have a
problem in the results, it seems the contact between surfaces is not detected
automatically, I want
On 8/18/21 7:18 PM, Michael Lee wrote:
BlockVector & system_rhs;
system_rhs.reinit(dofs_per_block);
For normal vectors, I can use *((double(*)[10]) v.begin()), but it does not
work for BlockVector. I'm using VS code.
Try *( (double(*)[10]) v.block(0).begin())
Best
W.
Sorry, it seems to work using *((double(*)[20])residual.begin()). From: Michael LeeSent: Wednesday, August 18, 2021 7:18 PMTo: deal.II User GroupSubject: [deal.II] How to watch the variable of BlockVector type? I want to check if I have system_rhs assembled correctly. Is there an easy way to
I want to check if I have *system_rhs* assembled correctly. Is there an
easy way to watch the contents of this BlockVector during debugging?
BlockVector & system_rhs;
system_rhs.reinit(dofs_per_block);
For normal vectors, I can use *((double(*)[10]) v.begin()), but it does
>
> > 2. The ExceptionBase::print_stack_trace() function does not print
> > anything if I catch an exception thrown using AssetThrow(cond,
> > dealii::ExcMessage()) or even my own MyException class derived from
> > ExceptionBase. How do I ensure that the stack trace is "populated". I
> > tried
Thanks a lot that helped, I applied symmetric<2,dim> tensor because it is a
heat problem. I also implemented the convection boundary condition. I have
a problem in the results, it seems the contact between surfaces is not
detected automatically, I want to model such contact.Please I will be