Searched refs:fdm_diagonal (Results 1 – 1 of 1) sorted by relevance
3037 CeedScalar *q_data_array, *fdm_diagonal; in CeedOperatorCreateFDMElementInverse() local3039 CeedCall(CeedCalloc(num_comp * num_nodes, &fdm_diagonal)); in CeedOperatorCreateFDMElementInverse()3043 if (interp) fdm_diagonal[c * num_nodes + n] = 1.0; in CeedOperatorCreateFDMElementInverse()3047 fdm_diagonal[c * num_nodes + n] += lambda[i]; in CeedOperatorCreateFDMElementInverse()3050 …if (fabs(fdm_diagonal[c * num_nodes + n]) < fdm_diagonal_bound) fdm_diagonal[c * num_nodes + n] = … in CeedOperatorCreateFDMElementInverse()3059 …q_data_array[(e * num_comp + c) * num_nodes + n] = 1. / (elem_avg[e] * fdm_diagonal[c * num_nodes … in CeedOperatorCreateFDMElementInverse()3064 CeedCall(CeedFree(&fdm_diagonal)); in CeedOperatorCreateFDMElementInverse()