It is difficult to write it as a single integral. The operation is similar to the split-step fourier method, i.e. transforming the column vector f(r) once using g(rho)=2\pi\int_0^\infty rf(r)J_0(2\pi\rho r)dr, multiplying it with a vector, and transforming it back using f(r) = 2\pi\int_0^\infty\rho g(\rho)J_0(2\pi\rho r)d\rho The operation is for radially symmetric systems, i.e. with z along the x-axis, and r along the y-axis. When starting on the left border with f_0, i.e. at position z = 0, doing the operation mentioned above gives the values for the nodes at z = 1, when enumerating the nodes from 0 to n along the z axis, and having equidistant nodes along the z-axis. Those integrals can be replaced by a matrix-vector-multiplication, thus making it easier to implement numerically. Does that make sense?
Am Freitag, 12. Juli 2019 18:24:22 UTC+2 schrieb Daniel Arndt: > > Maxi, > > can you clarify the operator evaluation you want to perform in > mathematical terms (maybe as an integral)? > > Best, > Daniel > -- 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/da921bbb-aad9-4ee9-af37-24971c981016%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.