Is this 
what you are looking for?



On Monday, June 22, 2020 at 1:32:43 AM UTC-4, Jimmy Ho wrote:
> Dear All,
> I am trying to write a code to give initial condition for a vector of 
> FullMatrix<double> objects. The initial value of each entry of each matrix 
> is uniquely defined by the corresponding nodal location. Hence, I am 
> looking for an easy way to iterate all nodes (including mid-side nodes) in 
> the mesh and access their coordinates. Is there a way to do so?
> I am aware of the Interpolate() function in VectorTools. However, the 
> documentation for this function show that the input function to interpolate 
> must have the same number of components as the number of components in the 
> vector-valued finite element. This is not the case for me, since the 
> initial values of the matrix are not directly linked to the DoFs of the 
> finite element. I also considered making the entries of the matrices as 
> extra DoFs of the finite element, then I will be able to use Interpolate() 
> and some component masks to assign initial values. But this idea is 
> sub-optimal, since in subsequent calculations, the values of the entries 
> can be directly computed using other DoFs of the finite element as well as 
> their previous values. So there's really no point posing them as additional 
> DoFs of the finite element. Is there any suggestions on how to achieve this?
> Thanks a lot for any discussions and suggestions! Your help is much 
> appreciated!
> Best,
> Jimmy

The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see 
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 dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 

Reply via email to