Lines Matching full:0
47 const CeedScalar *X = in[0], *J = in[1], *w = in[2]; in SetupMassGeoSphere()
49 CeedScalar *q_data = out[0]; in SetupMassGeoSphere()
52 CeedPragmaSIMD for (CeedInt i = 0; i < Q; i++) { in SetupMassGeoSphere()
54 const CeedScalar xx[3][1] = {{X[i + 0 * Q]}, {X[i + 1 * Q]}, {X[i + 2 * Q]}}; in SetupMassGeoSphere()
57 // 0 3 in SetupMassGeoSphere()
61 {J[i + Q * 0], J[i + Q * 3]}, in SetupMassGeoSphere()
67 const CeedScalar mod_xx_sq = xx[0][0] * xx[0][0] + xx[1][0] * xx[1][0] + xx[2][0] * xx[2][0]; in SetupMassGeoSphere()
69 for (int j = 0; j < 3; j++) { in SetupMassGeoSphere()
70 for (int k = 0; k < 3; k++) { in SetupMassGeoSphere()
71 xx_sq[j][k] = 0; in SetupMassGeoSphere()
72 … for (int l = 0; l < 1; l++) xx_sq[j][k] += xx[j][l] * xx[k][l] / (sqrt(mod_xx_sq) * mod_xx_sq); in SetupMassGeoSphere()
77 …{1. / sqrt(mod_xx_sq) - xx_sq[0][0], -xx_sq[0][1], -xx_sq[0][2] … in SetupMassGeoSphere()
78 …{-xx_sq[1][0], 1. / sqrt(mod_xx_sq) - xx_sq[1][1], -xx_sq[1][2] … in SetupMassGeoSphere()
79 …{-xx_sq[2][0], -xx_sq[2][1], 1. / sqrt(mod_xx_sq) - xx… in SetupMassGeoSphere()
83 for (int j = 0; j < 3; j++) { in SetupMassGeoSphere()
84 for (int k = 0; k < 2; k++) { in SetupMassGeoSphere()
85 dxdX[j][k] = 0; in SetupMassGeoSphere()
86 for (int l = 0; l < 3; l++) dxdX[j][k] += dxdxx[j][l] * dxxdX[l][k]; in SetupMassGeoSphere()
91 const CeedScalar J[3][1] = {{dxdX[1][0] * dxdX[2][1] - dxdX[2][0] * dxdX[1][1]}, in SetupMassGeoSphere()
92 {dxdX[2][0] * dxdX[0][1] - dxdX[0][0] * dxdX[2][1]}, in SetupMassGeoSphere()
93 {dxdX[0][0] * dxdX[1][1] - dxdX[1][0] * dxdX[0][1]}}; in SetupMassGeoSphere()
95 const CeedScalar mod_J = sqrt(J[0][0] * J[0][0] + J[1][0] * J[1][0] + J[2][0] * J[2][0]); in SetupMassGeoSphere()
96 q_data[i + Q * 0] = mod_J * w[i]; in SetupMassGeoSphere()
98 return 0; in SetupMassGeoSphere()