Home
last modified time | relevance | path

Searched refs:sum (Results 1 – 25 of 53) sorted by relevance

123

/libCEED/tests/
H A Dt305-basis.c43 CeedScalar sum = 0; in main() local
44 for (int k = 0; k < p; k++) sum += M[p * i + k] * X[p * k + j]; in main()
45 work_array[p * i + j] = sum; in main()
50 CeedScalar sum = 0; in main() local
51 for (int k = 0; k < p; k++) sum += X[p * k + i] * work_array[p * k + j]; in main()
52 M[p * i + j] = sum; in main()
68 CeedScalar sum = 0; in main() local
69 for (int k = 0; k < p; k++) sum += K[p * i + k] * X[p * k + j]; in main()
70 work_array[p * i + j] = sum; in main()
75 CeedScalar sum = 0; in main() local
[all …]
H A Dt302-basis.c14 CeedScalar sum = 0.0; in main() local
44 sum = 0.0; in main()
45 for (CeedInt j = 0; j < p; j++) sum += collocated_gradient_1d[j + p * i] * x_2[j]; in main()
46 …if (fabs(sum - 2 * q_ref[i]) > 100 * CEED_EPSILON) printf("Error in collocated gradient %f != %f\n… in main()
59 sum = 0.0; in main()
60 for (CeedInt j = 0; j < p + 2; j++) sum += collocated_gradient_1d[j + (p + 2) * i] * x_2[j]; in main()
61 …if (fabs(sum - 2 * q_ref[i]) > 100 * CEED_EPSILON) printf("Error in collocated gradient %f != %f\n… in main()
H A Dt304-basis.c26 CeedScalar sum = 0; in main() local
27 …for (int k = 0; k < p; k++) sum += interpolation[p * k + i] * quadrature_weights[k] * interpolatio… in main()
28 M[p * i + j] = sum; in main()
29 Q[p * i + j] = sum; in main()
38 CeedScalar sum = 0; in main() local
39 for (int k = 0; k < p; k++) sum += Q[p * i + k] * lambda[k] * Q[p * j + k]; in main()
40 Q_lambda_Qt[p * i + j] = sum; in main()
H A Dt552-operator.c109 CeedScalar sum = 0.; in main() local
113 sum += v_array[i]; in main()
116 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Area: 2.0… in main()
128 CeedScalar sum = 0.; in main() local
131 sum += v_array[i]; in main()
134 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Fine Grid: %f != True Area: 2.0\n… in main()
143 CeedScalar sum = 0.; in main() local
145 sum = 0.; in main()
147 sum += v_array[i]; in main()
150 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Area: 2.0… in main()
H A Dt553-operator.c100 CeedScalar sum = 0.; in main() local
104 sum += v_array[i]; in main()
107 …if (fabs(sum - 1.) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Area: 1.0… in main()
119 CeedScalar sum = 0.; in main() local
123 sum += v_array[i]; in main()
126 …if (fabs(sum - 1.) > 1000. * CEED_EPSILON) printf("Computed Area Fine Grid: %f != True Area: 1.0\n… in main()
135 CeedScalar sum = 0.; in main() local
139 sum += v_array[i]; in main()
142 …if (fabs(sum - 1.) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Area: 1.0… in main()
H A Dt550-operator.c104 CeedScalar sum = 0.; in main() local
108 sum += v_array[i]; in main()
111 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Area: 2.0… in main()
123 CeedScalar sum = 0.; in main() local
127 sum += v_array[i]; in main()
131 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Fine Grid: %f != True Area: 2.0\n… in main()
139 CeedScalar sum = 0.; in main() local
143 sum += v_array[i]; in main()
146 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Area: 2.0… in main()
H A Dt551-operator.c121 CeedScalar sum = 0.; in main() local
125 sum += v_array[i]; in main()
128 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Area: 2.0… in main()
140 CeedScalar sum = 0.; in main() local
144 sum += v_array[i]; in main()
147 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Fine Grid: %f != True Area: 2.0\n… in main()
156 CeedScalar sum = 0.; in main() local
160 sum += v_array[i]; in main()
163 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Area: 2.0… in main()
H A Dt554-operator.c167 CeedScalar sum = 0.; in main() local
171 sum += v_array[i]; in main()
174 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Area: 2.0… in main()
186 CeedScalar sum = 0.; in main() local
189 sum += v_array[i]; in main()
193 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Fine Grid: %f != True Area: 2.0\n… in main()
201 CeedScalar sum = 0.; in main() local
205 sum += v_array[i]; in main()
208 …if (fabs(sum - 2.) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Area: 2.0… in main()
H A Dt505-operator.c87 CeedScalar sum = 0.; in main() local
90 for (CeedInt i = 0; i < num_nodes_u; i++) sum += v_array[i]; in main()
92 if (fabs(sum - 1.) > 1000. * CEED_EPSILON) printf("Computed Area: %f != True Area: 1.0\n", sum); in main()
102 CeedScalar sum = -num_nodes_u; in main() local
105 for (CeedInt i = 0; i < num_nodes_u; i++) sum += v_array[i]; in main()
107 if (fabs(sum - 1.) > 1000. * CEED_EPSILON) printf("Computed Area: %f != True Area: 1.0\n", sum); in main()
H A Dt598-operator.c209 CeedScalar sum = 0.; in main() local
213 sum += v_array[i]; in main()
216 …if (fabs(sum - num_elem) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Are… in main()
228 CeedScalar sum = 0.; in main() local
232 sum += v_array[i]; in main()
236 …if (fabs(sum - num_elem) > 1000. * CEED_EPSILON) printf("Computed Area Fine Grid: %f != True Area:… in main()
244 CeedScalar sum = 0.; in main() local
248 sum += v_array[i]; in main()
251 …if (fabs(sum - num_elem) > 1000. * CEED_EPSILON) printf("Computed Area Coarse Grid: %f != True Are… in main()
H A Dt331-basis.c59 CeedScalar sum = 0.; in main() local
61 sum += u_array[i]; in main()
63 if (fabs(sum) > 100. * CEED_EPSILON) printf("sum of array %f != %f\n", sum, 0.0); in main()
H A Dt322-basis.c47 CeedScalar sum = 0; in main() local
51 for (int i = 0; i < q; i++) sum += v_array[i] * w_array[i]; in main()
52 if (fabs(sum - 17. / 24.) > 100. * CEED_EPSILON) printf("%f != %f\n", sum, 17. / 24.); in main()
H A Dt312-basis.c22 CeedScalar sum, error, pint[ALEN(p) + 1]; in main() local
67 sum = 0; in main()
68 for (CeedInt i = 0; i < q; i++) sum += w_array[i] * u_q_array[i]; in main()
75 error = sum - Eval(1, ALEN(pint), pint) + Eval(-1, ALEN(pint), pint); in main()
78 …printf("Error %e sum %g exact %g\n", error, sum, Eval(1, ALEN(pint), pint) - Eval(-1, ALEN(pint)… in main()
H A Dt414-qfunction.c76 CeedScalar sum = 0; in main() local
77 for (CeedInt i = 0; i < num_qpts; i++) sum += v_array[i + c * num_qpts]; in main()
78 if (fabs(sum - (c + 1)) > 10 * CEED_EPSILON) { in main()
80 …%" CeedInt_FMT "D volume error in component %" CeedInt_FMT ": %f != %f\n", dim, c, sum, (c + 1.0)); in main()
H A Dt524-operator.c177 CeedScalar sum = 0.; in main() local
180 for (CeedInt i = 0; i < num_dofs; i++) sum += v_array[i]; in main()
182 if (fabs(sum - 1.) > 1000. * CEED_EPSILON) printf("Computed Area: %f != True Area: 1.0\n", sum); in main()
192 CeedScalar sum = -num_dofs; in main() local
195 for (CeedInt i = 0; i < num_dofs; i++) sum += v_array[i]; in main()
197 if (fabs(sum - 1.) > 1000. * CEED_EPSILON) printf("Computed Area: %f != True Area: 1.0\n", sum); in main()
H A Dt415-qfunction.c80 CeedScalar sum = 0; in main() local
82 for (CeedInt g = 0; g < dim; g++) sum += v_array[i + (g * num_comp + c) * num_qpts]; in main()
84 if (fabs(sum - dim * (c + 1)) > 10 * CEED_EPSILON) { in main()
86 …dInt_FMT "D volume error in component %" CeedInt_FMT ": %f != %f\n", dim, c, sum, dim * (c + 1.0)); in main()
H A Dt503-operator.c88 CeedScalar sum = 0.; in main() local
91 for (CeedInt i = 0; i < num_nodes_u; i++) sum += v_array[i]; in main()
93 if (fabs(sum - 1.) > 1000. * CEED_EPSILON) printf("Computed Area: %f != True Area: 1.0\n", sum); in main()
H A Dt501-operator.c86 CeedScalar sum = 0.; in main() local
89 for (CeedInt i = 0; i < num_nodes_u; i++) sum += v_array[i]; in main()
91 if (fabs(sum - 1.) > 1000. * CEED_EPSILON) printf("Computed Area: %f != True Area: 1.0\n", sum); in main()
H A Dt599-operator.c119 CeedScalar sum = 0.0; in main() local
123 for (CeedInt i = 0; i < num_nodes; i++) sum += v_array[i]; in main()
126 if (fabs(sum - 4.0 * num_elem) > CEED_EPSILON * 5e3) { in main()
127 …printf("Incorrect area computed, %g != %g (abs error %g)\n", sum, 4.0 * num_elem, fabs(sum - 4.0 *… in main()
H A Dt590-operator.c97 CeedScalar sum = 0.0; in main() local
101 for (CeedInt i = 0; i < num_nodes; i++) sum += v_array[i]; in main()
104 …if (fabs(sum - 4.0 * num_elem) > CEED_EPSILON * 5e3) printf("Incorrect area computed, %f != %f\n",… in main()
H A Dt511-operator.c106 CeedScalar sum = 0.; in main() local
109 for (CeedInt i = 0; i < num_dofs; i++) sum += v_array[i]; in main()
111 if (fabs(sum - 1.) > 1000. * CEED_EPSILON) printf("Computed Area: %f != True Area: 1.0\n", sum); in main()
H A Dt591-operator.c167 CeedScalar sum = 0.0; in main() local
171 for (CeedInt i = 0; i < num_nodes; i++) sum += v_array[i]; in main()
174 …if (fabs(sum - 1.0 * num_elem) > CEED_EPSILON * 5e3) printf("Incorrect area computed, %f != %f\n",… in main()
H A Dt521-operator.c177 CeedScalar sum = 0.; in main() local
180 for (CeedInt i = 0; i < num_dofs; i++) sum += v_array[i]; in main()
182 if (fabs(sum - 1.) > 1000. * CEED_EPSILON) printf("Computed Area: %f != True Area: 1.0\n", sum); in main()
/libCEED/julia/LibCEED.jl/docs/src/
H A Dindex.md75 // Compute and print the sum of the entries of 'v' giving the mesh volume.
86 # Compute and print the sum of the entries of 'v' giving the mesh volume.
87 vol = witharray_read(sum, v, MEM_HOST)
92 // Initialize 'u' with sum of coordinates, x+y+z.
109 sum!(u_host, x_host)))
/libCEED/examples/solids/qfunctions/
H A Dfinite-strain-neo-hookean.h37 CeedScalar sum = 0; in log1p_series_shifted() local
40 sum -= log(2.) / 2; in log1p_series_shifted()
43 sum += log(2.) / 2; in log1p_series_shifted()
49 sum += y; in log1p_series_shifted()
51 sum += y / 3; in log1p_series_shifted()
53 sum += y / 5; in log1p_series_shifted()
55 sum += y / 7; in log1p_series_shifted()
56 return 2 * sum; in log1p_series_shifted()

123