Method commitSensitivity

SRC/element/brick/BbarBrickWithSensitivity.cpp:1993–2216

int
BbarBrickWithSensitivity::commitSensitivity(int gradNumber, int numGrads)
{
/*	const Vector &disp1 = theNodes[0]->getDispSensitivity(gradNumber);
	const Vector &disp2 = theNodes[1]->getDispSensitivity(gradNumber);
	const Vector &disp3 = theNodes[2]->getDispSensitivity(gradNumber);
	const Vector &disp4 = theNodes[3]->getDispSensitivity(gradNumber);

  static double u[2][4];

	u[0][0] = disp1(0);
	u[1][0] = disp1(1);
	u[0][1] = disp2(0);
	u[1][1] = disp2(1);
	u[0][2] = disp3(0);
	u[1][2] = disp3(1);
	u[0][3] = disp4(0);
	u[1][3] = disp4(1);


	refer
		sens(0) = theNodes[0]->getDispSensitivity(1,gradNumber);
		sens(1) = theNodes[0]->getDispSensitivity(2,gradNumber);
		sens(2) = theNodes[1]->getDispSensitivity(1,gradNumber);
		sens(3) = theNodes[1]->getDispSensitivity(2,gradNumber);
		sens(4) = theNodes[2]->getDispSensitivity(1,gradNumber);
		sens(5) = theNodes[2]->getDispSensitivity(2,gradNumber);
		sens(6) = theNodes[3]->getDispSensitivity(1,gradNumber);
		sens(7) = theNodes[3]->getDispSensitivity(2,gradNumber);

*/

		//strains ordered : eps11, eps22, eps33, 2*eps12, 2*eps23, 2*eps31


	  static int markGUU=0;

	   const int ndm = 3 ;

	   const int ndf = 3 ;

	   const int nstress = 6 ;

	   const int numberNodes = 8 ;

	   const int numberGauss = 8 ;

	   const int nShape = 4 ;

	  int i, j, k, p, q ;
//	  int jj, kk ;

	  int success ;

	   double volume ;

	   double xsj ;  // determinant jacaobian matrix