I work on FEM in Mathematica and want to speed up the computation when dealing with very fined mesh (i.e., a huge number of DOF). I found that a very much longer time is consumed when the Coefficients of PDE contains interpolation functions, compared to the case that the Coefficients are constants. To be moreRead more