48 CEED_QFUNCTION_HELPER void MaxEntry(const CeedScalar *A, CeedInt N, CeedInt *max_idx_row, CeedInt *…  in MaxEntry()  argument
50   *j_max               = max_idx_row[*i_max];  in MaxEntry()
53     CeedInt j = max_idx_row[i];  in MaxEntry()
159 …void ApplyRot(CeedScalar *A, CeedInt N, CeedInt i, CeedInt j, CeedInt *max_idx_row, CeedScalar *ro…  in ApplyRot()  argument
175     if (i == max_idx_row[w]) max_idx_row[w] = MaxEntryRow(A, N, w);  in ApplyRot()
176     else if (fabs(A[w * N + i]) > fabs(A[w * N + max_idx_row[w]])) max_idx_row[w] = i;  in ApplyRot()
188   max_idx_row[i] = MaxEntryRow(A, N, i);  in ApplyRot()
193     if (j == max_idx_row[w]) max_idx_row[w] = MaxEntryRow(A, N, w);  in ApplyRot()
194     else if (fabs(A[w * N + j]) > fabs(A[w * N + max_idx_row[w]])) max_idx_row[w] = j;  in ApplyRot()
198     if (j == max_idx_row[w]) max_idx_row[w] = MaxEntryRow(A, N, w);  in ApplyRot()
199     else if (fabs(A[w * N + j]) > fabs(A[w * N + max_idx_row[w]])) max_idx_row[w] = j;  in ApplyRot()
205   max_idx_row[j] = MaxEntryRow(A, N, j);  in ApplyRot()
281 …Int Diagonalize(CeedScalar *A, CeedInt N, CeedScalar *eval, CeedScalar *evec, CeedInt *max_idx_row,  in Diagonalize()  argument
289   for (CeedInt i = 0; i < N - 1; i++) max_idx_row[i] = MaxEntryRow(A, N, i);  in Diagonalize()
296     MaxEntry(A, N, max_idx_row, &i, &j);  in Diagonalize()
301       max_idx_row[i] = MaxEntryRow(A, N, i);  in Diagonalize()
307     ApplyRot(A, N, i, j, max_idx_row, rotmat_cst);  // Apply this rotation to the A matrix.  in Diagonalize()
323 …Diagonalize3(CeedScalar A[3][3], CeedScalar eval[3], CeedScalar evec[3][3], CeedInt max_idx_row[3],  in Diagonalize3()
325 …(CeedScalar *)A, 3, (CeedScalar *)eval, (CeedScalar *)evec, (CeedInt *)max_idx_row, sort_criteria,…  in Diagonalize3()