@@ -845,7 +845,8 @@ int CeedBasisCreateHdiv(Ceed ceed, CeedElemTopology topo, CeedInt num_comp, Ceed
845
845
@param[in] num_nodes Total number of nodes (dofs per element)
846
846
@param[in] num_qpts Total number of quadrature points
847
847
@param[in] interp Row-major (dim * num_qpts * num_nodes) matrix expressing the values of basis functions at quadrature points
848
- @param[in] curl Row-major (cdim * num_qpts * num_nodes, cdim = 1 if dim < 3 else dim) matrix expressing curl of basis functions at quadrature points
848
+ @param[in] curl Row-major (cdim * num_qpts * num_nodes, cdim = 1 if dim < 3 else dim) matrix expressing curl of basis functions at quadrature
849
+ points
849
850
@param[in] q_ref Array of length num_qpts * dim holding the locations of quadrature points on the reference element
850
851
@param[in] q_weight Array of length num_qpts holding the quadrature weights on the reference element
851
852
@param[out] basis Address of the variable where the newly created CeedBasis will be stored.
@@ -1094,7 +1095,7 @@ int CeedBasisApply(CeedBasis basis, CeedInt num_elem, CeedTransposeMode t_mode,
1094
1095
break ;
1095
1096
case CEED_EVAL_NONE :
1096
1097
case CEED_EVAL_INTERP :
1097
- qdim = (fe_space == CEED_FE_SPACE_H1 ) ? 1 : dim ;
1098
+ qdim = (fe_space == CEED_FE_SPACE_H1 ) ? 1 : dim ;
1098
1099
bad_dims = ((t_mode == CEED_TRANSPOSE && (u_length < num_elem * num_comp * num_qpts * qdim || v_length < num_elem * num_comp * num_nodes )) ||
1099
1100
(t_mode == CEED_NOTRANSPOSE && (v_length < num_elem * num_qpts * num_comp * qdim || u_length < num_elem * num_comp * num_nodes )));
1100
1101
break ;
@@ -1107,7 +1108,7 @@ int CeedBasisApply(CeedBasis basis, CeedInt num_elem, CeedTransposeMode t_mode,
1107
1108
(t_mode == CEED_NOTRANSPOSE && (v_length < num_elem * num_qpts * num_comp || u_length < num_elem * num_comp * num_nodes )));
1108
1109
break ;
1109
1110
case CEED_EVAL_CURL :
1110
- cdim = (dim < 3 ) ? 1 : dim ;
1111
+ cdim = (dim < 3 ) ? 1 : dim ;
1111
1112
bad_dims = ((t_mode == CEED_TRANSPOSE && (u_length < num_elem * num_comp * num_qpts * cdim || v_length < num_elem * num_comp * num_nodes )) ||
1112
1113
(t_mode == CEED_NOTRANSPOSE && (v_length < num_elem * num_qpts * num_comp * cdim || u_length < num_elem * num_comp * num_nodes )));
1113
1114
break ;
@@ -1871,26 +1872,27 @@ CeedPragmaOptimizeOn
1871
1872
}
1872
1873
CeedPragmaOptimizeOn
1873
1874
1874
- /**
1875
- @brief Apply Householder Q matrix
1875
+ /**
1876
+ @brief Apply Householder Q matrix
1876
1877
1877
- Compute A = Q A, where Q is mxm and A is mxn.
1878
+ Compute A = Q A, where Q is mxm and A is mxn.
1878
1879
1879
- @param[in,out] A Matrix to apply Householder Q to, in place
1880
- @param[in] Q Householder Q matrix
1881
- @param[in] tau Householder scaling factors
1882
- @param[in] t_mode Transpose mode for application
1883
- @param[in] m Number of rows in A
1884
- @param[in] n Number of columns in A
1885
- @param[in] k Number of elementary reflectors in Q, k<m
1886
- @param[in] row Row stride in A
1887
- @param[in] col Col stride in A
1880
+ @param[in,out] A Matrix to apply Householder Q to, in place
1881
+ @param[in] Q Householder Q matrix
1882
+ @param[in] tau Householder scaling factors
1883
+ @param[in] t_mode Transpose mode for application
1884
+ @param[in] m Number of rows in A
1885
+ @param[in] n Number of columns in A
1886
+ @param[in] k Number of elementary reflectors in Q, k<m
1887
+ @param[in] row Row stride in A
1888
+ @param[in] col Col stride in A
1888
1889
1889
- @return An error code: 0 - success, otherwise - failure
1890
+ @return An error code: 0 - success, otherwise - failure
1890
1891
1891
- @ref Developer
1892
- **/
1893
- int CeedHouseholderApplyQ (CeedScalar * A , const CeedScalar * Q , const CeedScalar * tau , CeedTransposeMode t_mode , CeedInt m , CeedInt n , CeedInt k ,
1892
+ @ref Developer
1893
+ **/
1894
+ int
1895
+ CeedHouseholderApplyQ (CeedScalar * A , const CeedScalar * Q , const CeedScalar * tau , CeedTransposeMode t_mode , CeedInt m , CeedInt n , CeedInt k ,
1894
1896
CeedInt row , CeedInt col ) {
1895
1897
CeedScalar * v ;
1896
1898
CeedCall (CeedMalloc (m , & v ));
@@ -1931,4 +1933,4 @@ int CeedMatrixMatrixMultiply(Ceed ceed, const CeedScalar *mat_A, const CeedScala
1931
1933
return CEED_ERROR_SUCCESS ;
1932
1934
}
1933
1935
1934
- /// @}
1936
+ /// @}
0 commit comments