| /honee/src/ |
| H A D | honee-file.c | 19 PetscErrorCode HoneeCheckFilenameExtension(MPI_Comm comm, const char filename[], const char extensi… in HoneeCheckFilenameExtension() argument 25 PetscCheck(ext_len, comm, PETSC_ERR_ARG_WRONG, "Zero-size extension: %s", extension); in HoneeCheckFilenameExtension() 62 MPI_Comm comm; in HoneeLoadInitialCondition() local 67 PetscCall(PetscObjectGetComm((PetscObject)Q, &comm)); in HoneeLoadInitialCondition() 68 PetscCall(HoneeCheckFilenameExtension(comm, filename, ".bin", &isBin)); in HoneeLoadInitialCondition() 69 PetscCall(HoneeCheckFilenameExtension(comm, filename, ".cgns", &isCGNS)); in HoneeLoadInitialCondition() 72 PetscCall(PetscViewerBinaryOpen(comm, filename, FILE_MODE_READ, &viewer)); in HoneeLoadInitialCondition() 82 PetscCall(PetscViewerCGNSOpen(comm, filename, FILE_MODE_READ, &viewer)); in HoneeLoadInitialCondition() 93 if (!set) PetscCall(PetscPrintf(comm, "WARNING: Couldn't find solution time in file\n")); in HoneeLoadInitialCondition() 103 PetscCheck(isFlowSolution, comm, PETSC_ERR_FILE_UNEXPECTED, in HoneeLoadInitialCondition() [all …]
|
| H A D | honee-meshtransform.c | 11 static PetscErrorCode GetYNodeLocs(const MPI_Comm comm, const char path[PETSC_MAX_PATH_LEN], PetscR… in GetYNodeLocs() argument 19 PetscCall(PetscFOpen(comm, path, "r", &fp)); in GetYNodeLocs() 20 PetscCall(PetscSynchronizedFGets(comm, fp, char_array_len, line)); in GetYNodeLocs() 35 PetscCall(PetscSynchronizedFGets(comm, fp, char_array_len, line)); in GetYNodeLocs() 37 PetscCheck(ndims == dims[1], comm, PETSC_ERR_FILE_UNEXPECTED, in GetYNodeLocs() 43 PetscCall(PetscFClose(comm, fp)); in GetYNodeLocs() 59 static PetscErrorCode HoneeMeshTransform_PlateMesh(MPI_Comm comm, DM dm, PetscReal growth, PetscInt… in HoneeMeshTransform_PlateMesh() argument 112 PetscCheck(*num_node_locs >= faces[1] + 1, comm, PETSC_ERR_FILE_UNEXPECTED, in HoneeMeshTransform_PlateMesh() 116 PetscCall(PetscPrintf(comm, in HoneeMeshTransform_PlateMesh() 143 MPI_Comm comm = PetscObjectComm((PetscObject)dm); in HoneeMeshTransformFromOptions() local [all …]
|
| H A D | misc.c | 98 MPI_Comm comm = PetscObjectComm((PetscObject)Q); in RegressionTest() local 103 …PetscCheck(strcmp(app_ctx->test_file_path, "") != 0, comm, PETSC_ERR_FILE_READ, "File for regressi… in RegressionTest() 104 PetscCall(PetscViewerBinaryOpen(comm, app_ctx->test_file_path, FILE_MODE_READ, &viewer)); in RegressionTest() 266 MPI_Comm comm = PetscObjectComm((PetscObject)ts); in PrintRunInfo() local 273 PetscCallMPI(MPI_Comm_rank(comm, &rank)); in PrintRunInfo() 274 PetscCallMPI(MPI_Comm_size(comm, &comm_size)); in PrintRunInfo() 275 PetscCall(PetscPrintf(comm, in PrintRunInfo() 290 PetscCall(PetscPrintf(comm, in PrintRunInfo() 305 PetscCall(PetscPrintf(comm, in PrintRunInfo() 332 PetscCall(PetscPrintf(comm, in PrintRunInfo() [all …]
|
| H A D | cloptions.c | 39 PetscErrorCode ISO8601TimeDurationToSeconds(MPI_Comm comm, const char *string, time_t *duration) { in ISO8601TimeDurationToSeconds() argument 49 …PetscCheck(isdigit(string[i]) || string[i] == ':', comm, PETSC_ERR_SUP, "Time duration may only in… in ISO8601TimeDurationToSeconds() 64 …SETERRQ(comm, PETSC_ERR_SUP, "Recieved %d ':' delimited entries, expect either 1, 2, or 3", num_it… in ISO8601TimeDurationToSeconds() 72 MPI_Comm comm = honee->comm; in ProcessCommandLineOptions() local 86 PetscCall(PetscPrintf(comm, "\n########################################\n\n")); in ProcessCommandLineOptions() 87 … PetscCall(PetscPrintf(comm, "HONEE documentation may be found at https://honee.phypid.org\n\n")); in ProcessCommandLineOptions() 88 PetscCall(PetscPrintf(comm, "########################################\n\n")); in ProcessCommandLineOptions() 91 PetscCheck(num_options > 0, comm, PETSC_ERR_USER_INPUT, in ProcessCommandLineOptions() 100 PetscOptionsBegin(comm, NULL, "HONEE - High-Order Navier-stokes Equation Evaluator", NULL); in ProcessCommandLineOptions() 138 MPI_Comm_rank(comm, &rank); in ProcessCommandLineOptions() [all …]
|
| H A D | honee.c | 25 PetscErrorCode HoneeInit(MPI_Comm comm, Honee *honee) { in HoneeInit() argument 30 …PetscCall(PetscHeaderCreate(honee_, HONEE_CLASSID, "Honee", "HONEE", "Honee", comm, HoneeDestroy, … in HoneeInit() 33 honee_->comm = PETSC_COMM_WORLD; in HoneeInit() 54 MPI_Comm comm = honee_->comm; in HoneeDestroy() local 81 …PetscCheck(CeedDestroy(&ceed) == CEED_ERROR_SUCCESS, comm, PETSC_ERR_LIB, "Destroying Ceed object … in HoneeDestroy()
|
| H A D | setupdm.c | 16 MPI_Comm comm = honee->comm; in CreateDM() local 22 PetscCall(HoneeCheckFilenameExtension(comm, filename, ".cgns", &isCGNS)); in CreateDM() 26 PetscCall(DMPlexCreateFromFile(comm, filename, "HONEE", PETSC_TRUE, dm)); in CreateDM() 29 PetscCall(DMCreate(comm, dm)); in CreateDM() 46 PetscCallMPI(MPI_Comm_size(comm, &num_ranks)); in CreateDM()
|
| H A D | setupts.c | 76 MPI_Comm comm = PetscObjectComm((PetscObject)dm); in Surface_Forces_NS() local 110 …CallMPI(MPIU_Allreduce(MPI_IN_PLACE, reaction_force, dim * num_walls, MPIU_SCALAR, MPI_SUM, comm)); in Surface_Forces_NS() 245 PetscCall(PetscViewerBinaryOpen(honee->comm, file_path, FILE_MODE_WRITE, &viewer)); in WriteOutput() 332 MPI_Comm comm; in TSPostStep_MaxWallTime() local 339 PetscCall(PetscObjectGetComm((PetscObject)ts, &comm)); in TSPostStep_MaxWallTime() 340 PetscCallMPI(MPI_Comm_rank(comm, &rank)); in TSPostStep_MaxWallTime() 343 PetscCallMPI(MPI_Bcast(&is_wall_time_exceeded, 1, MPI_C_BOOL, 0, comm)); in TSPostStep_MaxWallTime() 378 MPI_Comm comm = PetscObjectComm((PetscObject)ts); in HoneeTSEvaluationSolutions() local 395 PetscCall(PetscOptionsCreateViewer(comm, NULL, NULL, option, &viewer, &format, &is_viewer_set)); in HoneeTSEvaluationSolutions() 397 …PetscCall(PetscPrintf(comm, "\n\nWARNING: Viewer not set for TSEvaluationSolutions. Set %s to save… in HoneeTSEvaluationSolutions() [all …]
|
| H A D | dm-utils.c | 384 PetscErrorCode PetscFECreateLagrangeFromOptions(MPI_Comm comm, PetscInt dim, PetscInt num_comp, Pet… in PetscFECreateLagrangeFromOptions() argument 394 PetscCall(PetscSpaceCreate(comm, &fe_space)); in PetscFECreateLagrangeFromOptions() 406 PetscCall(PetscSpaceCreate(comm, &fe_space_end)); in PetscFECreateLagrangeFromOptions() 412 PetscCall(PetscSpaceCreate(comm, &fe_space_side)); in PetscFECreateLagrangeFromOptions() 428 PetscCall(PetscSpaceCreate(comm, &fe_space)); in PetscFECreateLagrangeFromOptions() 447 PetscCall(PetscDualSpaceCreate(comm, &fe_dual_space)); in PetscFECreateLagrangeFromOptions() 467 PetscOptionsBegin(comm, NULL, "Uniform quadrature check", NULL); in PetscFECreateLagrangeFromOptions() 470 …PetscCheck(!use_uniform || (is_tensor || dim == 1), comm, PETSC_ERR_SUP, "Can only use uniform qua… in PetscFECreateLagrangeFromOptions() 688 static inline PetscErrorCode Create1DTabulation_Tensor(MPI_Comm comm, PetscFE fe, PetscTabulation *… in Create1DTabulation_Tensor() argument 712 PetscCall(PetscFECreateLagrangeFromOptions(comm, 1, num_comp, !is_tensor, order, in Create1DTabulation_Tensor() [all …]
|
| H A D | monitor_totalkineticenergy.c | 135 MPI_Comm comm; in TSMonitor_TotalKineticEnergy() local 144 PetscCall(PetscObjectGetComm((PetscObject)ts, &comm)); in TSMonitor_TotalKineticEnergy() 145 PetscCallMPI(MPI_Comm_rank(comm, &rank)); in TSMonitor_TotalKineticEnergy() 156 …uce(MPI_IN_PLACE, monitor_ctx->sum_values, monitor_ctx->num_comps, MPIU_SCALAR, MPI_SUM, 0, comm)); in TSMonitor_TotalKineticEnergy() 157 …_ctx->sum_values, monitor_ctx->sum_values, monitor_ctx->num_comps, MPIU_SCALAR, MPI_SUM, 0, comm)); in TSMonitor_TotalKineticEnergy() 179 …PetscCheck(strftime(buf, sizeof(buf), "%FT%T%z", localtime(&t)), comm, PETSC_ERR_SYS, "strftime() … in TSMonitor_TotalKineticEnergy() 181 …PetscCheck(strftime(buf, sizeof(buf), "%Z", localtime(&t)), comm, PETSC_ERR_SYS, "strftime() call … in TSMonitor_TotalKineticEnergy()
|
| H A D | diff_flux_projection.c | 89 …SETERRQ(honee->comm, PETSC_ERR_ARG_WRONG, "Should not reach here with div_diff_flux_projection_met… in DivDiffFluxProjectionCreate() 127 MPI_Comm comm = PetscObjectComm((PetscObject)projection->dm); in DivDiffFluxProjectionSetup_Direct() local 133 PetscCheck(diff_flux_proj->CreateRHSOperator_Direct, comm, PETSC_ERR_ARG_WRONGSTATE, in DivDiffFluxProjectionSetup_Direct() 165 PetscCall(KSPCreate(comm, &projection->ksp)); in DivDiffFluxProjectionSetup_Direct() 201 MPI_Comm comm = PetscObjectComm((PetscObject)projection->dm); in DivDiffFluxProjectionSetup_Indirect() local 215 PetscCheck(diff_flux_proj->CreateRHSOperator_Indirect, comm, PETSC_ERR_ARG_WRONGSTATE, in DivDiffFluxProjectionSetup_Indirect() 237 PetscCall(KSPCreate(comm, &projection->ksp)); in DivDiffFluxProjectionSetup_Indirect() 284 PetscCheck(qf_calc_divergence, comm, PETSC_ERR_SUP, in DivDiffFluxProjectionSetup_Indirect()
|
| H A D | monitor_cfl.c | 141 MPI_Comm comm; in TSMonitor_Cfl() local 149 PetscCall(PetscObjectGetComm((PetscObject)ts, &comm)); in TSMonitor_Cfl() 159 PetscCall(PetscGlobalMinMaxReal(comm, part_minmax, global_minmax)); in TSMonitor_Cfl() 176 …PetscCheck(strftime(buf, sizeof(buf), "%FT%T%z", localtime(&t)), comm, PETSC_ERR_SYS, "strftime() … in TSMonitor_Cfl() 178 …PetscCheck(strftime(buf, sizeof(buf), "%Z", localtime(&t)), comm, PETSC_ERR_SYS, "strftime() call … in TSMonitor_Cfl()
|
| /honee/problems/ |
| H A D | stg_shur14.c | 31 PetscErrorCode CalcCholeskyDecomp(MPI_Comm comm, PetscInt nprofs, const CeedScalar Rij[6][nprofs], … in CalcCholeskyDecomp() argument 41 PetscCheck(!isnan(Cij[0][i]) && !isnan(Cij[1][i]) && !isnan(Cij[2][i]), comm, PETSC_ERR_FP, in CalcCholeskyDecomp() 59 static PetscErrorCode ReadStgInflow(const MPI_Comm comm, const char path[PETSC_MAX_PATH_LEN], StgSh… in ReadStgInflow() argument 68 PetscCall(PhastaDatFileOpen(comm, path, char_array_len, dims, &fp)); in ReadStgInflow() 77 PetscCall(PetscSynchronizedFGets(comm, fp, char_array_len, line)); in ReadStgInflow() 79 PetscCheck(ndims == dims[1], comm, PETSC_ERR_FILE_UNEXPECTED, in ReadStgInflow() 95 …PetscCheck(wall_dist[i] >= 0, comm, PETSC_ERR_FILE_UNEXPECTED, "Distance to wall in %s cannot be n… in ReadStgInflow() 96 …PetscCheck(lt[i] >= 0, comm, PETSC_ERR_FILE_UNEXPECTED, "Turbulent length scale in %s cannot be ne… in ReadStgInflow() 97 …PetscCheck(eps[i] >= 0, comm, PETSC_ERR_FILE_UNEXPECTED, "Turbulent dissipation in %s cannot be ne… in ReadStgInflow() 101 PetscCall(CalcCholeskyDecomp(comm, stg_ctx->nprofs, rij, cij)); in ReadStgInflow() [all …]
|
| H A D | shocktube.c | 15 MPI_Comm comm = honee->comm; in PRINT_SHOCKTUBE() local 18 PetscCall(PetscPrintf(comm, in PRINT_SHOCKTUBE() 30 MPI_Comm comm = honee->comm; in NS_SHOCKTUBE() local 64 PetscOptionsBegin(comm, NULL, "Options for SHOCKTUBE problem", NULL); in NS_SHOCKTUBE() 74 … PetscCall(PetscPrintf(comm, "Warning! -stab supg not implemented for the shocktube problem. \n")); in NS_SHOCKTUBE() 77 PetscCall(PetscPrintf(comm, "Warning! -yzb only implemented for explicit timestepping. \n")); in NS_SHOCKTUBE()
|
| H A D | eulervortex.c | 15 MPI_Comm comm = honee->comm; in PRINT_EULER_VORTEX() local 21 PetscCall(PetscPrintf(comm, in PRINT_EULER_VORTEX() 58 MPI_Comm comm = honee->comm; in NS_EULER_VORTEX() local 99 PetscOptionsBegin(comm, NULL, "Options for EULER_VORTEX problem", NULL); in NS_EULER_VORTEX() 116 PetscCall(PetscPrintf(comm, "Warning! Use -stab supg only with -implicit\n")); in NS_EULER_VORTEX() 119 …PetscCall(PetscPrintf(comm, "Warning! Background velocity vector for -euler_test t1 and -euler_tes… in NS_EULER_VORTEX()
|
| H A D | gaussianwave.c | 17 MPI_Comm comm = honee->comm; in NS_GAUSSIAN_WAVE() local 31 PetscOptionsBegin(comm, NULL, "Options for GAUSSIAN_WAVE problem", NULL); in NS_GAUSSIAN_WAVE() 34 …PetscCheck(narray == 3, comm, PETSC_ERR_ARG_SIZ, "-epicenter should recieve array of size 3, inste… in NS_GAUSSIAN_WAVE()
|
| H A D | taylorgreen.c | 14 MPI_Comm comm = honee->comm; in NS_TAYLOR_GREEN() local 33 PetscOptionsBegin(comm, NULL, "Options for TAYLOR-GREEN problem", NULL); in NS_TAYLOR_GREEN()
|
| H A D | advection.c | 22 MPI_Comm comm = honee->comm; in PRINT_ADVECTION() local 32 PetscCall(PetscPrintf(comm, in PRINT_ADVECTION() 45 …PetscCall(PetscPrintf(comm, " Background Wind : %f,%f\n", wind[0], wind[1])); in PRINT_ADVECTION() 48 …PetscCall(PetscPrintf(comm, " Background Wind : %f,%f,%f\n", wind[0], wind[1… in PRINT_ADVECTION() 53 …PetscCall(PetscPrintf(comm, " Initial Condition Type : %s\n", AdvDifICTypes[setup_c… in PRINT_ADVECTION() 61 …PetscCall(PetscPrintf(comm, " Bubble Continuity : %s\n", AdvDifBubbleContinuit… in PRINT_ADVECTION() 64 …PetscCall(PetscPrintf(comm, " Wave Type : %s\n", AdvDifWaveTypes[setup… in PRINT_ADVECTION() 188 …PetscCheck(qf_rhs_volume, honee->comm, PETSC_ERR_SUP, "%s not valid for DM of dimension %" PetscIn… in DivDiffFluxProjectionCreateRHS_Direct_AdvDif() 336 …PetscCheck(qf_rhs, honee->comm, PETSC_ERR_SUP, "%s not valid for DM of dimension %" PetscInt_FMT, … in DivDiffFluxProjectionCreateRHS_Indirect_AdvDif() 389 MPI_Comm comm = honee->comm; in NS_ADVECTION() local [all …]
|
| H A D | channel.c | 40 MPI_Comm comm = honee->comm; in NS_CHANNEL() local 59 PetscOptionsBegin(comm, NULL, "Options for CHANNEL problem", NULL); in NS_CHANNEL() 120 …PetscCall(PetscPrintf(comm, "WARNING! Channel flow with Inflow and Outflow is currently broken.\n"… in NS_CHANNEL() 132 …PetscCall(PetscPrintf(comm, "WARNING! Channel flow with Inflow and Outflow is currently broken.\n"… in NS_CHANNEL()
|
| H A D | blasius.c | 134 MPI_Comm comm = honee->comm; in NS_BLASIUS() local 158 PetscOptionsBegin(comm, NULL, "Options for BLASIUS problem", NULL); in NS_BLASIUS() 170 …PetscCheck(3 <= N && N <= BLASIUS_MAX_N_CHEBYSHEV, comm, PETSC_ERR_ARG_OUTOFRANGE, "-n_chebyshev %… in NS_BLASIUS() 221 PetscCall(SetupStg(comm, dm, problem, honee, weakT, S_infty.Y.temperature, S_infty.Y.pressure)); in NS_BLASIUS() 225 …_var == STATEVAR_CONSERVATIVE) || (honee->app_ctx->test_type == TESTTYPE_DIFF_FILTER), honee->comm, in NS_BLASIUS()
|
| H A D | densitycurrent.c | 16 MPI_Comm comm = honee->comm; in NS_DENSITY_CURRENT() local 44 PetscOptionsBegin(comm, NULL, "Options for DENSITY_CURRENT problem", NULL); in NS_DENSITY_CURRENT()
|
| /honee/include/ |
| H A D | honee-file.h | 13 PetscErrorCode HoneeCheckFilenameExtension(MPI_Comm comm, const char filename[], const char extensi… 19 PetscErrorCode PhastaDatFileOpen(const MPI_Comm comm, const char path[], const PetscInt char_array_… 20 PetscErrorCode PhastaDatFileGetNRows(const MPI_Comm comm, const char path[], PetscInt *nrows); 21 PetscErrorCode PhastaDatFileReadToArrayReal(const MPI_Comm comm, const char path[], PetscReal array…
|
| /honee/examples/ |
| H A D | navierstokes.c | 22 MPI_Comm comm; in main() local 25 comm = PETSC_COMM_WORLD; in main() 27 PetscCall(HoneeInit(comm, &honee)); in main() 45 …PetscCheck(CeedInit(app_ctx->ceed_resource, &ceed) == CEED_ERROR_SUCCESS, comm, PETSC_ERR_LIB, "Ce… in main() 47 …PetscCheck(CeedSetErrorHandler(ceed, CeedErrorStore) == CEED_ERROR_SUCCESS, comm, PETSC_ERR_LIB, "… in main()
|
| /honee/src/smartsim/ |
| H A D | smartsim.c | 31 PetscCallMPI(MPI_Comm_rank(honee->comm, &rank)); in SmartSimTrainingSetup() 52 PetscOptionsBegin(honee->comm, NULL, "Options for SmartSim integration", NULL); in SmartSimSetup() 58 PetscCallMPI(MPI_Comm_rank(honee->comm, &rank)); in SmartSimSetup()
|
| H A D | solution.c | 16 MPI_Comm comm = PetscObjectComm((PetscObject)ts); in TSMonitor_SmartSimSolutionSetup() local 23 PetscOptionsBegin(comm, NULL, "SmartSim Solution Writing", NULL); in TSMonitor_SmartSimSolutionSetup() 57 PetscCallMPI(MPI_Comm_rank(honee->comm, &rank)); in TSMonitor_SmartSimSolution()
|
| H A D | sgs_dd_training.c | 189 PetscOptionsBegin(honee->comm, NULL, "SGS Data-Driven Training Options", NULL); in SGS_DD_TrainingSetup() 217 PetscCallMPI(MPI_Comm_rank(honee->comm, &rank)); in SGS_DD_TrainingSetup() 218 PetscCallMPI(MPI_Comm_size(honee->comm, &num_ranks)); in SGS_DD_TrainingSetup() 228 PetscCall(PetscGlobalMinMaxInt(honee->comm, local_min_max, global_min_max)); in SGS_DD_TrainingSetup() 302 PetscCallMPI(MPI_Comm_rank(honee->comm, &rank)); in TSMonitor_SGS_DD_Training() 408 …PetscCall(PetscPrintf(honee->comm, "-- Simulation stopped by 'check-run' tensor in Redis database\… in TSPostStep_SGS_DD_Training()
|