Lines Matching refs:user
29 PetscErrorCode CreateStatsDM(User user, ProblemData problem, PetscInt degree) { in CreateStatsDM() argument
30 user->spanstats.num_comp_stats = TURB_NUM_COMPONENTS; in CreateStatsDM()
34 MPI_Comm comm = PetscObjectComm((PetscObject)user->dm); in CreateStatsDM()
42 PetscCall(DMGetBoundingBox(user->dm, domain_min, domain_max)); in CreateStatsDM()
43 user->spanstats.span_width = domain_max[2] - domain_min[2]; in CreateStatsDM()
52 PetscCall(DMPlexGetIsoperiodicFaceSF(user->dm, &num_isoperiodicface, &isoperiodicface)); in CreateStatsDM()
63 PetscCall(DMCreateLabel(user->dm, "Periodic Face")); in CreateStatsDM()
64 PetscCall(DMGetLabel(user->dm, "Periodic Face", &label)); in CreateStatsDM()
70 PetscCall(DMGetLabel(user->dm, "Face Sets", &label)); in CreateStatsDM()
73 PetscCall(DMPlexLabelComplete(user->dm, label)); 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 …all(DMSetupByOrder_FEM(PETSC_TRUE, PETSC_TRUE, user->app_ctx->degree, 1, user->app_ctx->q_extra, 1… in CreateStatsDM()
99 user->spanstats.dm)); in CreateStatsDM()
102 PetscCall(DMGetLocalSection(user->spanstats.dm, §ion)); in CreateStatsDM()
185 PetscErrorCode SpanStatsSetupDataCreate(Ceed ceed, User user, CeedData ceed_data, ProblemData probl… in SpanStatsSetupDataCreate() argument
186 DM dm = user->spanstats.dm; in SpanStatsSetupDataCreate()
187 CeedInt num_comp_x, num_comp_stats = user->spanstats.num_comp_stats; in SpanStatsSetupDataCreate()
296 PetscErrorCode SetupL2ProjectionStats(Ceed ceed, User user, CeedData ceed_data, SpanStatsSetupData … in SetupL2ProjectionStats() argument
299 CeedInt q_data_size, num_comp_stats = user->spanstats.num_comp_stats; in SetupL2ProjectionStats()
315 …etscCall(OperatorApplyContextCreate(NULL, user->spanstats.dm, ceed, op_proj_rhs, NULL, NULL, NULL,… in SetupL2ProjectionStats()
316 …edOperatorCreateLocalVecs(op_proj_rhs, DMReturnVecType(user->spanstats.dm), PETSC_COMM_SELF, &user… 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()
360 PetscErrorCode CreateStatisticCollectionOperator(Ceed ceed, User user, CeedData ceed_data, SpanStat… in CreateStatisticCollectionOperator() argument
361 …CeedInt num_comp_stats = user->spanstats.num_comp_stats, num_comp_x = problem-… in CreateStatisticCollectionOperator()
372 switch (user->phys->state_var) { in CreateStatisticCollectionOperator()
384 if (user->spanstats.do_mms_test) { in CreateStatisticCollectionOperator()
394 PetscCallCeed(ceed, CeedQFunctionContextCreate(user->ceed, &collect_context)); in CreateStatisticCollectionOperator()
420 …PetscCallCeed(ceed, CeedOperatorGetContextFieldLabel(op_stats_collect, "solution time", &user->spa… in CreateStatisticCollectionOperator()
421 …PetscCallCeed(ceed, CeedOperatorGetContextFieldLabel(op_stats_collect, "previous time", &user->spa… 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()
436 PetscErrorCode SetupMMSErrorChecking(Ceed ceed, User user, CeedData ceed_data, SpanStatsSetupData s… in SetupMMSErrorChecking() argument
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()
478 PetscErrorCode TurbulenceStatisticsSetup(Ceed ceed, User user, CeedData ceed_data, ProblemData prob… in TurbulenceStatisticsSetup() argument
488 PetscCall(CreateStatsDM(user, problem, user->app_ctx->degree)); in TurbulenceStatisticsSetup()
491 PetscCall(SpanStatsSetupDataCreate(ceed, user, ceed_data, problem, &stats_data)); in TurbulenceStatisticsSetup()
494 …PetscCall(CreateStatsSF(ceed, ceed_data, stats_data, user->dm, user->spanstats.dm, &user->spanstat… in TurbulenceStatisticsSetup()
497 PetscCall(CreateStatisticCollectionOperator(ceed, user, ceed_data, stats_data, problem)); in TurbulenceStatisticsSetup()
500 PetscCall(SetupL2ProjectionStats(ceed, user, ceed_data, stats_data)); in TurbulenceStatisticsSetup()
502 …PetscCall(PetscOptionsGetBool(NULL, NULL, "-ts_monitor_turbulence_spanstats_mms", &user->spanstats… in TurbulenceStatisticsSetup()
503 if (user->spanstats.do_mms_test) { in TurbulenceStatisticsSetup()
504 PetscCall(SetupMMSErrorChecking(ceed, user, ceed_data, stats_data)); in TurbulenceStatisticsSetup()
509 PetscCall(PetscViewerGetType(user->app_ctx->turb_spanstats_viewer, &viewer_type)); in TurbulenceStatisticsSetup()
512 …PetscCall(PetscViewerSetOptionsPrefix(user->app_ctx->turb_spanstats_viewer, "ts_monitor_turbulence… in TurbulenceStatisticsSetup()
513 PetscCall(PetscViewerSetFromOptions(user->app_ctx->turb_spanstats_viewer)); in TurbulenceStatisticsSetup()
522 PetscErrorCode CollectStatistics(User user, PetscScalar solution_time, Vec Q) { in CollectStatistics() argument
523 SpanStatsData user_stats = user->spanstats; in CollectStatistics()
531 PetscCall(UpdateBoundaryValues(user, user->Q_loc, solution_time)); in CollectStatistics()
532 …PetscCallCeed(user->ceed, CeedOperatorSetContextDouble(user_stats.op_stats_collect_ctx->op, user_s… in CollectStatistics()
533 PetscCall(DMGlobalToLocal(user->dm, Q, INSERT_VALUES, user->Q_loc)); in CollectStatistics()
534 …PetscCall(ApplyAddCeedOperatorLocalToLocal(user->Q_loc, user_stats.Child_Stats_loc, user_stats.op_… in CollectStatistics()
536 …PetscCallCeed(user->ceed, CeedOperatorSetContextDouble(user_stats.op_stats_collect_ctx->op, user_s… in CollectStatistics()
543 PetscErrorCode ProcessStatistics(User user, Vec stats) { in ProcessStatistics() argument
544 SpanStatsData user_stats = user->spanstats; in ProcessStatistics()
576 PetscReal summing_duration = solution_time - user->app_ctx->cont_time; in ProcessStatistics()
592 User user = (User)ctx; in TSMonitor_TurbulenceStatistics() local
595 …PetscInt collect_interval = user->app_ctx->turb_spanstats_collect_interval, viewer_interval = user… in TSMonitor_TurbulenceStatistics()
600 if (steps == user->app_ctx->cont_steps) PetscFunctionReturn(PETSC_SUCCESS); in TSMonitor_TurbulenceStatistics()
605 PetscCall(CollectStatistics(user, solution_time, Q)); in TSMonitor_TurbulenceStatistics()
608 PetscCall(DMSetOutputSequenceNumber(user->spanstats.dm, steps, solution_time)); in TSMonitor_TurbulenceStatistics()
609 PetscCall(DMGetGlobalVector(user->spanstats.dm, &stats)); in TSMonitor_TurbulenceStatistics()
610 PetscCall(ProcessStatistics(user, stats)); in TSMonitor_TurbulenceStatistics()
611 if (user->app_ctx->test_type == TESTTYPE_NONE) { in TSMonitor_TurbulenceStatistics()
612 …PetscCall(PetscViewerPushFormat(user->app_ctx->turb_spanstats_viewer, user->app_ctx->turb_spanstat… in TSMonitor_TurbulenceStatistics()
613 PetscCall(VecView(stats, user->app_ctx->turb_spanstats_viewer)); in TSMonitor_TurbulenceStatistics()
614 PetscCall(PetscViewerPopFormat(user->app_ctx->turb_spanstats_viewer)); in TSMonitor_TurbulenceStatistics()
616 if (user->app_ctx->test_type == TESTTYPE_TURB_SPANSTATS && reason != TS_CONVERGED_ITERATING) { in TSMonitor_TurbulenceStatistics()
617 PetscCall(RegressionTest(user->app_ctx, 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()
634 PetscErrorCode TurbulenceStatisticsDestroy(User user, CeedData ceed_data) { in TurbulenceStatisticsDestroy() argument
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()