794   CeedScalar        *vals, *BTD_mat = NULL, *elem_mat = NULL, *elem_mat_b = NULL;  in CeedOperatorAssembleSingle()  local
798   CeedCall(CeedCalloc(elem_size_out * elem_size_in, &elem_mat));  in CeedOperatorAssembleSingle()
827                                            false, B_mat_in, elem_mat));  in CeedOperatorAssembleSingle()
832 …CeedCall(CeedMatrixMatrixMultiply(ceed, BTD_mat, B_mat_in, elem_mat, elem_size_out, elem_size_in, …  in CeedOperatorAssembleSingle()
844               elem_mat[i * elem_size_in + j] *= orient;  in CeedOperatorAssembleSingle()
851           memcpy(elem_mat_b, elem_mat, elem_size_out * elem_size_in * sizeof(CeedScalar));  in CeedOperatorAssembleSingle()
854 … elem_mat[i * elem_size_in + j] = elem_mat_b[i * elem_size_in + j] * elem_curl_orients[3 * i + 1] +  in CeedOperatorAssembleSingle()
865               elem_mat[i * elem_size_in + j] *= elem_orients[j] ? -1.0 : 1.0;  in CeedOperatorAssembleSingle()
872           memcpy(elem_mat_b, elem_mat, elem_size_out * elem_size_in * sizeof(CeedScalar));  in CeedOperatorAssembleSingle()
875 … elem_mat[i * elem_size_in + j] = elem_mat_b[i * elem_size_in + j] * elem_curl_orients[3 * j + 1] +  in CeedOperatorAssembleSingle()
885             vals[offset + count] = elem_mat[i * elem_size_in + j];  in CeedOperatorAssembleSingle()
897   CeedCall(CeedFree(&elem_mat));  in CeedOperatorAssembleSingle()