VecLoadIntoVector() takes as INPUT a vector like dmmg[DMMGGetLevels(dmmg)-1]->x That is a vector that comes from the finest level of the DMMG. Here you are passing a null into VecLoadIntoVector() in other words you are not passing a vector at all into the routine. This cannot possibly work. Your FormInitialGuess_physical(DMMG dmmg, Vec X) takes a DMMG as input, and a vector X. Likely this vector X is what you pass into VecLoadIntoVector().
Barry On Mar 10, 2010, at 9:50 AM, (Rebecca) Xuefei YUAN wrote: > 5 0x085e985a in VecLoadIntoVector (viewer=0x8a08d10, vec=0x0) > at vector.c:1016