Re: [deal.II] evaluating basis functions at random points

2021-12-13 Thread Daniel Arndt
Andrea,

have a look at GridTools::compute_point_locations(
https://www.dealii.org/developer/doxygen/deal.II/namespaceGridTools.html#a8e8bb9211264d2106758ac4d7184117e
).
Step-60 demonstrates its use.

Best,
Daniel

Am Mo., 13. Dez. 2021 um 16:14 Uhr schrieb Andrea Bonito :

> Howdy:
>
> this might be a simple question but I do not seem to find a way to
> efficiently evaluate basis functions at random points. More precisely, I am
> given random points in my computational domain and I need to assemble the
> terms
>
> phi_i(p) phi_j(p)  for all random points p and all i,j =1,..., n_dofs.
>
> As of now, I am going cell by cell and check whether p is inside the cell
> (including the boundary). If it is, I am creating a quadrature with one
> point in the reference cell (using the mapping).
>
> I am sure one of you has a better way... any help would be appreciated.
>
> Andrea
>
> --
> 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 dealii+unsubscr...@googlegroups.com.
> To view this discussion on the web visit
> https://groups.google.com/d/msgid/dealii/18f9edd4-c848-4983-9724-e87b31bd219cn%40googlegroups.com
> 
> .
>

-- 
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 dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/CAOYDWbL0sLUAKK2wOzu8uV2oKOhu-QY7OXSyQPUSaX1NVMxmsQ%40mail.gmail.com.


[deal.II] evaluating basis functions at random points

2021-12-13 Thread Andrea Bonito
Howdy:

this might be a simple question but I do not seem to find a way to 
efficiently evaluate basis functions at random points. More precisely, I am 
given random points in my computational domain and I need to assemble the 
terms 

phi_i(p) phi_j(p)  for all random points p and all i,j =1,..., n_dofs. 

As of now, I am going cell by cell and check whether p is inside the cell 
(including the boundary). If it is, I am creating a quadrature with one 
point in the reference cell (using the mapping).

I am sure one of you has a better way... any help would be appreciated.

Andrea

-- 
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 dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/18f9edd4-c848-4983-9724-e87b31bd219cn%40googlegroups.com.