As source term for calculations I am using data which is written in an external matrix (because it has been calculated by an external function), together with the x- and y-values for each matrix entry. My current approach for using this data is to use an interpolation function from GSL, and call it for each cell while looping over them. This gives me the value at that position, interpolated. Now, though, I would like to improve that approach. I was thinking about using the project-function and interpolating the matrix onto a component in my FESystem, but as far as I can see the project-function does not allow me to only change one component in my solution while leaving the other components untouched. Of course, I also could loop once over all cells before starting the main loop and solve the equation u = f, with f the value at each position, but I am not sure how efficient that is, especially when I can use the matrix-free vectorized approach for solving for the other components. Therefore, are there other (maybe more efficient) approaches which I can use for transferring and interpolating the values from the original matrix to something which I then can use in my deal.II-functions? Thanks!
-- 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 [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/290765bf-ae24-404d-9e36-64cb887a1e22%40googlegroups.com.
