This problem lies in the user-defined array "dphi_fwd[][]", which is set as std::vector<std::vector<RealGradient> >. Because the program is able to run when I drop off all "dphi_fwd" arrays. This dphi_fwd[][] is updated in different quadrature points and varied with time. So my question is setting up this sort of RealGradient array to take part in the computation.
I suppose the system given dphi[][] comes from the reference of fe->get_dphi(), this is should be unchanged. Is any ideal way to build up this array, which could be updated and memory-saving? Thanks in advance. Bests, Liang ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
