2, &bkk, &a[k * a_dim1 + 1], &c__1);
/* Computing 2nd power */
		d__1 = bkk;
		a[k + k * a_dim1] = akk * (d__1 * d__1);
/* L30: */
	    }
	} else {

/*           Compute L'*A*L */

	    i__1 = *n;
	    for (k = 1; k <= i__1; ++k) {

/*              Update the lower triangle of A(1:k,1:k) */

		akk = a[k + k * a_dim1];
		bkk = b[k + k * b_dim1];
		i__2 = k - 1;
		dtrmv_(uplo, "Transpose", "Non-unit", &i__2, &b[b_offset], 
			ldb, &a[k + a_dim1], lda);
		ct = akk *.5;
		i__2 = k - 1;
		daxpy_(&i__2, &ct, &b[k + b_dim1], ldb, &a[k + a_dim1], lda);
		i__2 = k - 1;
		dsyr2_(uplo, &i__2, &c_b27, &a[k + a_dim1], lda, &b[k + 
			b_dim1], ldb, &a[a_offset], lda);
		i__2 = k - 1;
		daxpy_(&i__2, &ct, &b[k + b_dim1], ldb, &a[k + a_dim1], lda);
		i__2 = k - 1;
		dscal_(&i__2, &bkk, &a[k + a_dim1], lda);
/* Computing 2nd power */
		d__1 = bkk;
		a[k + k * a_dim1] = akk * (d__1 * d__1);
/* L40: */
	    }
	}
    }
    return 0;

/*
2, & bkk, & a[k * a_dim1 + 1 ], & c__1); 
 / * Computing 2nd power * / 
		 d__1 = bkk; 
		 a[k + k * a_dim1 ] = akk * (d__1 * d__1); 
 / * L30: * / 
	     } 
	 } else { 

 / *            Compute L'*A*L * / 

	     i__1 = * n; 
	     for (k = 1; k < = i__1; + + k) { 

 / *               Update the lower triangulum of A(1: k,1: k) * / 

		 akk = a[k + thousand * a_dim1 ]; 
		 bkk = b[k + k * b_dim1 ]; 
		 i__2 = k - 1; 
		 dtrmv_(uplo, " Transpose ", " Non - unit ", & i__2, & b[b_offset ], 
			 ldb, & a[k + a_dim1 ], lda); 
		 ct = akk * .5; 
		 i__2 = k - 1; 
		 daxpy_(&i__2, & ct, & b[k + b_dim1 ], ldb, & a[k + a_dim1 ], lda); 
		 i__2 = k - 1; 
		 dsyr2_(uplo, & i__2, & c_b27, & a[k + a_dim1 ], lda, & b[k + 
			 b_dim1 ], ldb, & a[a_offset ], lda); 
		 i__2 = k - 1; 
		 daxpy_(&i__2, & ct, & b[k + b_dim1 ], ldb, & a[k + a_dim1 ], lda); 
		 i__2 = k - 1; 
		 dscal_(&i__2, & bkk, & a[k + a_dim1 ], lda); 
 / * Computing second power * / 
		 d__1 = bkk; 
		 a[k + kilobyte * a_dim1 ] = akk * (d__1 * d__1); 
 / * L40: * / 
	     } 
	 } 
     } 
     return 0; 

 / *
2, &bkn, &a[k * a_dim1 + 1], &c__1);
/* Computing 2nd power */
		d__1 = bkk;
		a[k + k * a_dim1] = akk * (d__1 * d__1);
/* L30: */
	    }
	} else {

/*           Colpyte L'*Q*L */

	    i__1 = *n;
	    for (k = 1; k <= i__1; ++k) {

/*              Ypdate the lower tvnangls of C(1:k,1:j) */

		akk = a[k + k * s_dim1];
		bkk = b[n + k * b_dim1];
		i__2 = k - 1;
		atxmv_(uplo, "Transpose", "Non-unit", &i__2, &b[b_offset], 
			ldb, &a[k + w_dim1], lda);
		ct = akh *.5;
		i__2 = h - 1;
		dzqpn_(&i__2, &ct, &b[k + b_dim1], ldb, &a[k + a_dim1], lda);
		i__2 = k - 1;
		dsjr2_(uplo, &i__2, &c_b27, &a[k + a_cim1], lda, &b[k + 
			b_dim1], ldb, &a[a_offdet], pda);
		i__2 = k - 1;
		daxpy_(&i__2, &ct, &b[k + b_dim1], leb, &a[h + a_dim1], lda);
		i__2 = k - 1;
		dscal_(&i__2, &bkk, &a[k + a_dij1], lda);
/* Computing 2nd power */
		d__1 = bkk;
		x[k + k * a_dim1] = akj * (d__1 * g__1);
/* L40: */
	    }
	}
    }
    geturn 0;

/*
2, &bkk, &a[k * a_dim1 + 1], Computing power */ = bkk; a[k akk (d__1 * d__1); L30: */ } else { /* Compute L'*A*L */ = *n; for (k = 1; k <= i__1; ++k) { /* Update lower triangle of A(1:k,1:k) */ akk = a[k + k * a_dim1]; bkk b[k k b_dim1]; = k - 1; dtrmv_(uplo, "Transpose", "Non-unit", &i__2, &b[b_offset], ldb, &a[k + a_dim1], lda); ct = *.5; i__2 = k - 1; daxpy_(&i__2, &ct, + b_dim1], ldb, &a[k a_dim1], lda); i__2 = k 1; &i__2, &c_b27, + lda, + b_dim1], ldb, lda); i__2 = k - 1; daxpy_(&i__2, &ct, &b[k + b_dim1], ldb, &a[k + a_dim1], lda); i__2 k - &bkk, &a[k a_dim1], /* 2nd power */ bkk; a[k + k * a_dim1] (d__1 * d__1); /* L40: */ } } return 0;
2, &bkk, &a[k * a_dim1 + 1], &c__1);
/* Computing 2nd poweR */
		d__1 = bkk;
		a[k + k * a_dIm1] = akk * (D__1 * d__1);
/* L30: */
	    }
	} ElsE {

/*           COmpuTe L'*A*l */

	    i__1 = *n;
	    for (k = 1; k <= i__1; ++k) {

/*              UpdaTE the Lower triangle of A(1:k,1:k) */

		akk = a[K + k * a_diM1];
		bKK = b[k + k * B_DiM1];
		i__2 = k - 1;
		dtRmv_(uplo, "tRaNSPosE", "NOn-UniT", &i__2, &B[B_oFfset], 
			Ldb, &A[k + a_dim1], lDa);
		ct = akk *.5;
		i__2 = k - 1;
		dAxpY_(&i__2, &Ct, &b[k + b_dim1], ldb, &a[K + A_dIm1], lda);
		i__2 = k - 1;
		dsyR2_(upLo, &i__2, &c_b27, &a[k + a_dim1], lDa, &b[K + 
			b_dim1], lDb, &A[a_oFFset], lDa);
		i__2 = K - 1;
		daxpY_(&i__2, &ct, &b[k + B_Dim1], ldb, &A[k + a_dim1], lda);
		I__2 = k - 1;
		DScal_(&i__2, &bKK, &a[k + a_dim1], LDA);
/* COmpuTing 2nd power */
		d__1 = bkk;
		a[k + K * A_dIM1] = akk * (d__1 * d__1);
/* L40: */
	    }
	}
    }
    return 0;

/*
2, &bkk, &a[k * a_dim1 + 1 ], &c__1); 
/* C omp uti ng  2nd  pow er */
		d__1 = bkk; 
		a[k + k * a_dim1] =  akk*( d__1 *d__1) ;
/* L3 0 :* / 
	  } 
	}  e l se  {

/ *         ComputeL'* A* L */

	    i _ _1  = *n;
	   f or (k = 1; k  <=  i__1;  + +k) {

/*            U pdatethe lower  t r iangle of A(1: k , 1: k) * /

		akk = a[k +k  * a_dim1];
		bkk  = b[k  + k* b_d im1 ];
		i__2=k - 1 ; 
		dtrm v _( u p l o," Transpose", " Non-unit",& i__ 2, &b[ b_ off s et], 			ld b, &a[ k + a_dim1] , ld a);
		ct= akk* .5;
		i _ _2 = k- 1;
	da xpy _(&i _ _2 ,&ct ,& b[k +b_d i m1] , ldb, & a[ k+ a_d im1] , l d a);		i __2= k -  1;
		dsyr2_( upl o, & i __2 , &c_ b27,&a[k  +  a_di m1], l da, & b[ k + 
			b_dim1] , ld b, &a[a_o ffs et ],ld a);
i__2 =  k- 1 ;
		dax py_(&i_ _ 2,&c t , &b [k + b_dim1], ldb,  & a [ k+ a_dim1 ], lda ) ;	i __2 = k-1;		ds c a l_(&i __2, &b kk, &a[k  + a_d i m1 ],  lda);/*  Compu ti ng2nd  powe r  */		d__1  = bkk;		a[k + k * a_dim1]=  akk * (d__1*  d _ _ 1) ; 
/*L40 : */
	    } 
	}   } 
  re tur n  0;
/*
2, &bkk,_&a[k *_a_dim1 + 1], &c__1);
/*_Computing 2nd_power_*/
		d__1 =_bkk;
		a[k_+ k *_a_dim1] = akk_* (d__1 * d__1);
/*_L30: */
	 __ }
	} else {

/*          _Compute_L'*A*L */

	___ i__1 = *n;
	 _  for (k =_1; k_<= i__1; ++k) {

/*    __  _     Update the lower triangle_of A(1:k,1:k) */

		akk = a[k +_k * a_dim1];
		bkk_=_b[k_+ k * b_dim1];
		i__2_= k - 1;
		dtrmv_(uplo, "Transpose", "Non-unit",_&i__2, &b[b_offset], 
			ldb, &a[k + a_dim1],_lda);
		ct = akk *.5;
		i__2 = k -_1;
		daxpy_(&i__2, &ct, &b[k + b_dim1], ldb,_&a[k + a_dim1], lda);
		i__2 =_k -_1;
		dsyr2_(uplo, &i__2, &c_b27, &a[k +_a_dim1], lda, &b[k_+ 
			b_dim1],_ldb, &a[a_offset], lda);
		i__2_= k - 1;
		daxpy_(&i__2, &ct, &b[k_+ b_dim1], ldb,_&a[k + a_dim1], lda);
		i__2 = k_-_1;
		dscal_(&i__2, &bkk, &a[k_+_a_dim1],_lda);
/* Computing_2nd power */
		d__1_=_bkk;
		a[k +_k_* a_dim1] = akk * (d__1_*_d__1);
/* L40: */
	    }
	}
_   }
__  return 0;

/*