1 // Copyright (c) 2017-2026, Lawrence Livermore National Security, LLC and other CEED contributors.
2 // All Rights Reserved. See the top-level LICENSE and NOTICE files for details.
4 // SPDX-License-Identifier: BSD-2-Clause
20   Ceed         ceed = user->ceed;  in ICs_FixMultiplicity()
26 …if (user->phys->ics_time_label) PetscCallCeed(ceed, CeedOperatorSetContextDouble(ceed_data->op_ics…  in ICs_FixMultiplicity()
27   PetscCall(ApplyCeedOperatorLocalToGlobal(NULL, Q, ceed_data->op_ics_ctx));  in ICs_FixMultiplicity()
29   PetscCallCeed(ceed, CeedElemRestrictionCreateVector(ceed_data->elem_restr_q, &mult_vec, NULL));  in ICs_FixMultiplicity()
31   // -- Get multiplicity  in ICs_FixMultiplicity()
34   PetscCallCeed(ceed, CeedElemRestrictionGetMultiplicity(ceed_data->elem_restr_q, mult_vec));  in ICs_FixMultiplicity()
41   // -- Fix multiplicity  in ICs_FixMultiplicity()
61     PetscCall(VecAXPY(Qbc, -1., Q_loc));  in SetBCsFromICs()
98   *out = -13;  // appease the overzealous GCC compiler warning Gods  in BinaryReadIntoInt()
123 …token == FLUIDS_FILE_TOKEN) {  // New style format; we're reading a file with step number and time…  in LoadFluidsBinaryVec()
152 …PetscCheck(strcmp(app_ctx->test_file_path, "") != 0, comm, PETSC_ERR_FILE_READ, "File for regressi…  in RegressionTest()
153   PetscCall(PetscViewerBinaryOpen(comm, app_ctx->test_file_path, FILE_MODE_READ, &viewer));  in RegressionTest()
159   PetscCall(VecAXPY(Q, -1.0, Q_ref));  in RegressionTest()
164   if (error > app_ctx->test_tol) {  in RegressionTest()
188   // Get |exact solution - obtained solution|  in PrintError()
190   PetscCall(VecAXPY(Q, -1.0, Q_exact));  in PrintError()
200 // Post-processing
207   if (problem->compute_exact_solution_error && user->app_ctx->test_type == TESTTYPE_NONE) {  in PostProcess()
214   if (user->app_ctx->test_type == TESTTYPE_NONE) {  in PostProcess()
220   PetscCall(VecViewFromOptions(Q, NULL, "-vec_view"));  in PostProcess()
223   if (user->app_ctx->test_type == TESTTYPE_SOLVER) {  in PostProcess()
224     PetscCall(RegressionTest(user->app_ctx, Q));  in PostProcess()
238   PetscCall(PetscViewerBinaryOpen(comm, app_ctx->cont_file, FILE_MODE_READ, &viewer));  in SetupICsFromBinary()
239   PetscCall(LoadFluidsBinaryVec(comm, viewer, Q, &app_ctx->cont_time, &app_ctx->cont_steps));  in SetupICsFromBinary()
284   PetscCall(DMDestroy(&context->dm));  in NodalProjectionDataDestroy()
285   PetscCall(KSPDestroy(&context->ksp));  in NodalProjectionDataDestroy()
287   PetscCall(OperatorApplyContextDestroy(context->l2_rhs_ctx));  in NodalProjectionDataDestroy()
373   Ceed     ceed = user->ceed;  in PrintRunInfo()
384                         "\n-- Navier-Stokes solver - libCEED + PETSc --\n"  in PrintRunInfo()
391   PetscCall(problem->print_info(user, problem, user->app_ctx));  in PrintRunInfo()
396   PetscCallCeed(ceed, CeedGetResource(user->ceed, &used_resource));  in PrintRunInfo()
397   PetscCallCeed(ceed, CeedGetPreferredMemType(user->ceed, &mem_type_backend));  in PrintRunInfo()
406   if (problem->dim == 2) box_faces_str[3] = '\0';  in PrintRunInfo()
407 …PetscCall(PetscOptionsGetString(NULL, NULL, "-dm_plex_box_faces", box_faces_str, sizeof(box_faces_…  in PrintRunInfo()
408   PetscCall(DMGetVecType(user->dm, &vec_type));  in PrintRunInfo()
414                         box_faces_str, vec_type, phys_ctx->implicit ? "implicit" : "explicit"));  in PrintRunInfo()
441   if (user->app_ctx->cont_steps) {  in PrintRunInfo()
447 …                   user->app_ctx->cont_file, user->app_ctx->cont_steps, user->app_ctx->cont_time));  in PrintRunInfo()
452   const CeedInt  num_P = user->app_ctx->degree + 1, num_Q = num_P + user->app_ctx->q_extra;  in PrintRunInfo()
453   PetscCall(DMGetGlobalVectorInfo(user->dm, &owned_dofs, &glob_dofs, NULL));  in PrintRunInfo()
454   PetscCall(DMGetLocalVectorInfo(user->dm, &local_dofs, NULL, NULL));  in PrintRunInfo()
461 …                    "    Global %" PetscInt_FMT "-DoF nodes                 : %" PetscInt_FMT "\n",  in PrintRunInfo()
463   // -- Get Partition Statistics  in PrintRunInfo()
468     PetscInt  median_index = comm_size % 2 ? comm_size / 2 : comm_size / 2 - 1;  in PrintRunInfo()
475       part_owned_dofs[1]             = gather_buffer[comm_size - 1];  // max  in PrintRunInfo()
479 …"    Global Vector %" PetscInt_FMT "-DoF nodes          : %" PetscInt_FMT ", %" PetscInt_FMT ", %"…  in PrintRunInfo()
488       part_local_dofs[1]             = gather_buffer[comm_size - 1];  // max  in PrintRunInfo()
492 …"    Local Vector %" PetscInt_FMT "-DoF nodes           : %" PetscInt_FMT ", %" PetscInt_FMT ", %"…  in PrintRunInfo()
505         PetscCall(DMGetSectionSF(user->dm, &sf));  in PrintRunInfo()
509           if (rranks[i] == rank) continue;  // Ignore same-part global->local transfers  in PrintRunInfo()
510           num_remote_roots_total += roffset[i + 1] - roffset[i];  in PrintRunInfo()
515           num_remote_leaves_total += ioffset[i + 1] - ioffset[i];  in PrintRunInfo()
523         part_boundary_dofs[1]           = gather_buffer[comm_size - 1];  // max  in PrintRunInfo()
527 …"    Ghost Interface %" PetscInt_FMT "-DoF nodes        : %" PetscInt_FMT ", %" PetscInt_FMT ", %"…  in PrintRunInfo()
537         part_neighbors[1]              = gather_buffer[comm_size - 1];  // max  in PrintRunInfo()
548         part_boundary_dofs[1]           = gather_buffer[comm_size - 1];  // max  in PrintRunInfo()
552 …"    Owned Interface %" PetscInt_FMT "-DoF nodes        : %" PetscInt_FMT ", %" PetscInt_FMT ", %"…  in PrintRunInfo()
562         part_neighbors[1]              = gather_buffer[comm_size - 1];  // max  in PrintRunInfo()