Commit 295b7ad8 authored by Giovanni La Mura's avatar Giovanni La Mura
Browse files

Fix access to VJ variable in GHIT (untested)

parent 50b7ecc5

src/libnptm/clu_subs.cpp

+2 −2
Original line number Diff line number Diff line
@@ -980,7 +980,7 @@ dcomplex ghit( 
	      int ny = l3 * l3 + lt34;

	      double aors = 1.0 * (l3 + lt34);

	      double f3j = (c1->v3j0[i3j0 - 1] * c1->v3j0[i3j0 - 1] * sqrt(aors)) * jsn;

	      cfun = (c1->vh[nbhj + lt34 - 1] * c1->vyhj[nby + ny - 1]) * f3j;

	      cfun = (c1->vj * c1->vyhj[nby + ny - 1]) * f3j;

	      csum += cfun;

	      jsn *= -1;

	      lt34 += 2;
@@ -998,7 +998,7 @@ dcomplex ghit( 
		int ny = l3 * l3  + lt40 + m1mm2;

		double aors = 1.0 * (l3 + lt40);

		double f3j = (c6->rac3j[il - 1] * c1->v3j0[i3j0 - 1] * sqrt(aors)) * jsn;

		cfun = (c1->vh[nbhj + lt40 - 1] * c1->vyhj[nby + ny - 1]) * f3j;

		cfun = (c1->vj * c1->vyhj[nby + ny - 1]) * f3j;

		csum += cfun;

	      }

	      // label 40