In stiffness_matrix file, in stiffness_matrix func qsh moved to one outer loop.

stiffness_matrix.py

@@ -8,7 +8,7 @@ ex = np.array((1, 0, 0))
 ey = np.array((0, 1, 0))
 ez = np.array((0, 0, 1))
 gauss_pw = [[[-0.577350269189626, 1],\
-             [-0.577350269189626, 1]],\
+             [0.577350269189626, 1]],\
                     [[-0.774596669241483, 0.555555555555555],\
                       [0, 0.888888888888888], \
                       [-0.774596669241483, 0.555555555555555]],\
@@ -366,7 +366,7 @@ def stiffness_matrix(surface, lobatto_pw, number_gauss_point=2, elastic_modulus=
             lag_xi2 = lagfunc(lobatto_pw, xi2)
             der_lag_dxi1 = der_lagfunc_dxi(lobatto_pw, xi1)
             der_lag_dxi2 = der_lagfunc_dxi(lobatto_pw, xi2)
-        
+            qsh = mat_transform_matrix(surface, xi1, xi2)
             for k in range(len_gauss):
                 t = gauss_pw_nparray[k, 0]
                 w3 =gauss_pw_nparray[k, 1]
@@ -374,7 +374,6 @@ def stiffness_matrix(surface, lobatto_pw, number_gauss_point=2, elastic_modulus=
                     der_lag_dxi1, der_lag_dxi2, t)
                 der_disp_var = der_disp (surface, lobatto_pw, der_disp_dxi_var, xi1,\
                     xi2, t)                
-                qsh = mat_transform_matrix(surface, xi1, xi2)
                 elastic_matrix = np.transpose(qsh).\
                     dot(elastic_matrix_planestress).dot(qsh)
                 b = b_matrix(lobatto_pw, der_disp_var)