Closed bendudson closed 14 years ago
Seems to be a bug in mesh.load(). Inserting the following in bout++.cpp after the load call gives a segfault on 4 or 8 processors:
Field3D a; a = 0.0; Field3D b = a;
Caused by mistyping "mesh->xend" instead of "mesh->yend" in Field3D::checkData(). Fixed by 67cfe462c2156a303cb0
Running with 4 or more processors, the lapd_drift example case segfaults. Running valgrind on a single processor doesn't give anything obvious. Adding a line like
Field3D tmp = Ni0;
after reading in Ni0 also gives a segfault. Am investigating...