[deal.II] Computing DFT

Dear all I want to compute the DFT of a PDE solution. I have seen some Fourier Series in deal.II https://www.dealii.org/developer/doxygen/deal.II/classFESeries_1_1Fourier.html and it says it is used to compute the DFT on the reference cell. I have seen examples of its use to estimate local s

Re: [deal.II] system_to_component_index() usage in face_worker for DG method

Thank you for the answer. I have been looking at a way of doing what you suggested and I found the function* FEInterfaceValues ::interface_dof_to_dof_indices()* which seems to be doing what I want. I provide it with the

Re: [deal.II] Computing DFT

On 7/23/21 3:00 AM, Praveen C wrote: I want to compute the DFT of a PDE solution. I have seen some Fourier Series in deal.II https://nam10.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwww.dealii.org%2Fdeveloper%2Fdoxygen%2Fdeal.II%2FclassFESeries_1_1Fourier.html&data=04%7C01%7CWolfgang.B

Re: [deal.II] system_to_component_index() usage in face_worker for DG method

On 7/23/21 8:59 AM, Sylvain Mathonnière wrote: Thank you for the answer. I have been looking at a way of doing what you suggested and I found the function*FEInterfaceValues

Re: [deal.II] system_to_component_index() usage in face_worker for DG method

Oh, I think I know what those values appears and it is my bad ... It is probably not 22116 but 22 and 116 and it is probably not 89130 but rather 89 and 130. I used something like *std::cout << dof_index[0] << std::endl* to get those values and sometimes they just get attached to each other lik

Re: [deal.II] system_to_component_index() usage in face_worker for DG method

Sylvain, Oh, I think I know what those values appears and it is my bad ... > It is probably not 22116 but 22 and 116 and it is probably not 89130 but > rather 89 and 130. I used something like *std::cout << dof_index[0] << > std::endl* to get those values and sometimes they just get attached to

Re: [deal.II] system_to_component_index() usage in face_worker for DG method

On 7/23/21 9:41 AM, Sylvain Mathonnière wrote: As for the 4294967295 occurence, is it because I am dealing with an interface which is a boundary as well ? Isn't this case handled by the *boundary_worker* and not the *face_worker* (using step-12 notation) ? No, it's because you are using DG e

Re: [deal.II] Computing DFT

Thank you. Related to this question. Can I save a parallel distributed solution using serialize and read it back as a serial triangulation/solution ? This may help me to sample the solution on a Cartesian grid and do FFT. Thanks praveen > On 23-Jul-2021, at 8:30 PM, Wolfgang Bangerth wrote: