74     integer i__, j, k;  in invert_()  local
95     for (j = 1; j <= i__1; ++j) {  in invert_()
97 	for (i__ = j; i__ <= i__2; ++i__) {  in invert_()
98 	    v[i__ + j * 10 - 11] = phi[i__ + j * phi_dim1];  in invert_()
100 	i__2 = j - 1;  in invert_()
102 	    save = v[j + k * 10 - 11] * v[k + k * 10 - 11];  in invert_()
104 	    for (i__ = j; i__ <= i__3; ++i__) {  in invert_()
105 		v[i__ + j * 10 - 11] -= v[i__ + k * 10 - 11] * save;  in invert_()
109 	if ((r__1 = v[j + j * 10 - 11], abs(r__1)) < 1e-10f) {  in invert_()
112 	rc[j] = psi[j];  in invert_()
113 	i__2 = j - 1;  in invert_()
115 	    rc[j] -= rc[k] * v[j + k * 10 - 11];  in invert_()
117 	v[j + j * 10 - 11] = 1.f / v[j + j * 10 - 11];  in invert_()
118 	rc[j] *= v[j + j * 10 - 11];  in invert_()
121 	r__2 = rc[j];  in invert_()
123 	rc[j] = max(r__1,-.999f);  in invert_()
129     for (i__ = j; i__ <= i__1; ++i__) {  in invert_()