30   user->spanstats.num_comp_stats = TURB_NUM_COMPONENTS;  in CreateStatsDM()
43   user->spanstats.span_width = domain_max[2] - domain_min[2];  in CreateStatsDM()
74 …  PetscCall(DMPlexFilter(user->dm, label, 1, PETSC_FALSE, PETSC_FALSE, NULL, &user->spanstats.dm));  in CreateStatsDM()
75 …PetscCall(DMSetCoordinateDisc(user->spanstats.dm, NULL, PETSC_TRUE));  // Ensure that a coordinate…  in CreateStatsDM()
77     PetscCall(DMPlexDistribute(user->spanstats.dm, 0, NULL, &parent_distributed_dm));  in CreateStatsDM()
80       PetscCall(DMDestroy(&user->spanstats.dm));  in CreateStatsDM()
81       user->spanstats.dm = parent_distributed_dm;  in CreateStatsDM()
88     PetscCall(DMPlexIsSimplex(user->spanstats.dm, &is_simplex));  in CreateStatsDM()
92   PetscCall(PetscObjectSetName((PetscObject)user->spanstats.dm, "Spanwise_Stats"));  in CreateStatsDM()
93   PetscCall(DMSetOptionsPrefix(user->spanstats.dm, "turbulence_spanstats_"));  in CreateStatsDM()
94   PetscCall(DMSetFromOptions(user->spanstats.dm));  in CreateStatsDM()
95   PetscCall(DMViewFromOptions(user->spanstats.dm, NULL, "-dm_view"));  in CreateStatsDM()
98 …, PETSC_TRUE, user->app_ctx->degree, 1, user->app_ctx->q_extra, 1, &user->spanstats.num_comp_stats,  in CreateStatsDM()
99                                user->spanstats.dm));  in CreateStatsDM()
102   PetscCall(DMGetLocalSection(user->spanstats.dm, &section));  in CreateStatsDM()
186   DM       dm = user->spanstats.dm;  in SpanStatsSetupDataCreate()
187   CeedInt  num_comp_x, num_comp_stats = user->spanstats.num_comp_stats;  in SpanStatsSetupDataCreate()
299   CeedInt             q_data_size, num_comp_stats = user->spanstats.num_comp_stats;  in SetupL2ProjectionStats()
315 …(OperatorApplyContextCreate(NULL, user->spanstats.dm, ceed, op_proj_rhs, NULL, NULL, NULL, NULL, &…  in SetupL2ProjectionStats()
316 …rCreateLocalVecs(op_proj_rhs, DMReturnVecType(user->spanstats.dm), PETSC_COMM_SELF, &user->spansta…  in SetupL2ProjectionStats()
317 …PetscCall(QDataGet(ceed, user->spanstats.dm, domain_label, label_value, stats_data->elem_restr_par…  in SetupL2ProjectionStats()
331     PetscCall(MatCreateCeed(user->spanstats.dm, user->spanstats.dm, op_mass, NULL, &mat_mass));  in SetupL2ProjectionStats()
333     PetscCall(KSPCreate(PetscObjectComm((PetscObject)user->spanstats.dm), &ksp));  in SetupL2ProjectionStats()
345     user->spanstats.ksp = ksp;  in SetupL2ProjectionStats()
361 …CeedInt                     num_comp_stats = user->spanstats.num_comp_stats, num_comp_x = problem-…  in CreateStatisticCollectionOperator()
384   if (user->spanstats.do_mms_test) {  in CreateStatisticCollectionOperator()
420 …dOperatorGetContextFieldLabel(op_stats_collect, "solution time", &user->spanstats.solution_time_la…  in CreateStatisticCollectionOperator()
421 …dOperatorGetContextFieldLabel(op_stats_collect, "previous time", &user->spanstats.previous_time_la…  in CreateStatisticCollectionOperator()
423 …PetscCall(OperatorApplyContextCreate(user->dm, user->spanstats.dm, user->ceed, op_stats_collect, u…  in CreateStatisticCollectionOperator()
424                                        &user->spanstats.op_stats_collect_ctx));  in CreateStatisticCollectionOperator()
426 …PetscCall(CeedOperatorCreateLocalVecs(op_stats_collect, DMReturnVecType(user->spanstats.dm), PETSC…  in CreateStatisticCollectionOperator()
427                                         &user->spanstats.Child_Stats_loc));  in CreateStatisticCollectionOperator()
428   PetscCall(VecZeroEntries(user->spanstats.Child_Stats_loc));  in CreateStatisticCollectionOperator()
437   CeedInt             num_comp_stats = user->spanstats.num_comp_stats, num_comp_x, q_data_size;  in SetupMMSErrorChecking()
447 …PetscCall(QDataGet(ceed, user->spanstats.dm, domain_label, label_value, stats_data->elem_restr_par…  in SetupMMSErrorChecking()
465 …PetscCall(OperatorApplyContextCreate(user->spanstats.dm, user->spanstats.dm, user->ceed, op_error,…  in SetupMMSErrorChecking()
466                                        &user->spanstats.mms_error_ctx));  in SetupMMSErrorChecking()
494 …all(CreateStatsSF(ceed, ceed_data, stats_data, user->dm, user->spanstats.dm, &user->spanstats.sf));  in TurbulenceStatisticsSetup()
502 …tionsGetBool(NULL, NULL, "-ts_monitor_turbulence_spanstats_mms", &user->spanstats.do_mms_test, NUL…  in TurbulenceStatisticsSetup()
503   if (user->spanstats.do_mms_test) {  in TurbulenceStatisticsSetup()
523   SpanStatsData user_stats = user->spanstats;  in CollectStatistics()
544   SpanStatsData      user_stats = user->spanstats;  in ProcessStatistics()
608       PetscCall(DMSetOutputSequenceNumber(user->spanstats.dm, steps, solution_time));  in TSMonitor_TurbulenceStatistics()
609       PetscCall(DMGetGlobalVector(user->spanstats.dm, &stats));  in TSMonitor_TurbulenceStatistics()
619       if (user->spanstats.do_mms_test && reason != TS_CONVERGED_ITERATING) {  in TSMonitor_TurbulenceStatistics()
622         PetscCall(ApplyCeedOperatorGlobalToGlobal(stats, error, user->spanstats.mms_error_ctx));  in TSMonitor_TurbulenceStatistics()
628       PetscCall(DMRestoreGlobalVector(user->spanstats.dm, &stats));  in TSMonitor_TurbulenceStatistics()
636   PetscCall(VecDestroy(&user->spanstats.Child_Stats_loc));  in TurbulenceStatisticsDestroy()
637   PetscCall(VecDestroy(&user->spanstats.Parent_Stats_loc));  in TurbulenceStatisticsDestroy()
639   PetscCall(OperatorApplyContextDestroy(user->spanstats.op_stats_collect_ctx));  in TurbulenceStatisticsDestroy()
640   PetscCall(OperatorApplyContextDestroy(user->spanstats.op_proj_rhs_ctx));  in TurbulenceStatisticsDestroy()
641   PetscCall(OperatorApplyContextDestroy(user->spanstats.mms_error_ctx));  in TurbulenceStatisticsDestroy()
643   PetscCall(KSPDestroy(&user->spanstats.ksp));  in TurbulenceStatisticsDestroy()
644   PetscCall(PetscSFDestroy(&user->spanstats.sf));  in TurbulenceStatisticsDestroy()
645   PetscCall(DMDestroy(&user->spanstats.dm));  in TurbulenceStatisticsDestroy()