@@ -227,31 +227,37 @@ averageToVertex(Omega_h::Mesh &mesh,
227227 return sizeField;
228228}
229229
230-
231230/* *
232- * @brief Computes the SPR (Superconvergent Patch Recovery) size field for mesh adaptation
231+ * @brief Computes the SPR (Superconvergent Patch Recovery) size field for mesh
232+ * adaptation
233233 *
234- * This function performs error estimation using the SPR method and computes desired element
235- * sizes based on the error between the original field and the recovered field. The resulting
236- * size field can be used to guide adaptive mesh refinement.
234+ * This function performs error estimation using the SPR method and computes
235+ * desired element sizes based on the error between the original field and the
236+ * recovered field. The resulting size field can be used to guide adaptive mesh
237+ * refinement.
237238 *
238- * @tparam EstimationT The estimation type containing mesh, fields, tolerance, and other parameters
239+ * @tparam EstimationT The estimation type containing mesh, fields, tolerance,
240+ * and other parameters
239241 * @tparam OmegahMeshField The Omega_h mesh field type for field evaluation
240242 * @tparam FieldElement The field element type for coordinate field operations
241243 *
242- * @param e (In) The estimation object containing the mesh, input field (eps), recovered field (eps_star),
243- * tolerance, and storage for intermediate results
244- * @param omf (In) The Omega_h mesh field object used for evaluating fields at integration points
245- * @param coordFe (In) The coordinate field element used for integration over mesh elements
244+ * @param e (In) The estimation object containing the mesh, input field (eps),
245+ * recovered field (eps_star), tolerance, and storage for
246+ * intermediate results
247+ * @param omf (In) The Omega_h mesh field object used for evaluating fields at
248+ * integration points
249+ * @param coordFe (In) The coordinate field element used for integration over
250+ * mesh elements
246251 *
247252 * @return A tuple containing:
248- * - Kokkos::View<MeshField::Real *>: The vertex-averaged size field indicating desired
249- * element sizes at mesh vertices
250- * - MeshField::Real: The sum of error norms raised to the power (2d/(2p+d)) over all elements,
251- * where d is the mesh dimension and p is the recovered field polynomial order
253+ * - Kokkos::View<MeshField::Real *>: The vertex-averaged size field
254+ * indicating desired element sizes at mesh vertices
255+ * - MeshField::Real: The sum of error norms raised to the power
256+ * (2d/(2p+d)) over all elements, where d is the mesh dimension and p is the
257+ * recovered field polynomial order
252258 */
253259template <typename EstimationT, typename OmegahMeshField, typename FieldElement>
254- std::tuple< Kokkos::View<MeshField::Real *>, MeshField::Real >
260+ std::tuple<Kokkos::View<MeshField::Real *>, MeshField::Real>
255261getSprSizeField (EstimationT &e, OmegahMeshField &omf, FieldElement &coordFe) {
256262 Error errorIntegrator (e, omf);
257263 errorIntegrator.process (coordFe);
0 commit comments