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
gra
>
> > 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 usin
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 no
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 watc
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.
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 t
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 wrote
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 “intersectin
