132   double fdx,fdy,fdz,fOx,fOy,fOz,fHx,fHy,fHz;  in eval()  local
337             fHz = 0.5*alpha * fdz;  in eval()
345             f[iH1].z += fHz;  in eval()
349             f[iH2].z += fHz;  in eval()
356               v[2] = x[i].z*fOz + xH1.z*fHz + xH2.z*fHz;  in eval()
358               v[4] = x[i].x*fOz + xH1.x*fHz + xH2.x*fHz;  in eval()
359               v[5] = x[i].y*fOz + xH1.y*fHz + xH2.y*fHz;  in eval()
396             fHz = 0.5*alpha * fdz;  in eval()
404             f[jH1].z += fHz;  in eval()
408             f[jH2].z += fHz;  in eval()
415               v[2] += x[j].z*fOz + xH1.z*fHz + xH2.z*fHz;  in eval()
417               v[4] += x[j].x*fOz + xH1.x*fHz + xH2.x*fHz;  in eval()
418               v[5] += x[j].y*fOz + xH1.y*fHz + xH2.y*fHz;  in eval()