fluids: petsc-ops accept CEED_VECTOR_NONE, NULL CeedVector

fluids: Add CeedOperatorCreateLocalVecs

fluids: Change if...SETERRQ to PetscCheck

fluids: Enable TSSetErrorIfStepFails when testing

- Misc error fixes too
- And style change to match index sizes for similar loops
- Other similar nested for loops have `j<5` and `k<3`, so cha

fluids: Enable TSSetErrorIfStepFails when testing

- Misc error fixes too
- And style change to match index sizes for similar loops
- Other similar nested for loops have `j<5` and `k<3`, so changed
this one to match

show more ...

fluids: Add flag for turning on SGS model

fluids: Apply nodal SGS tensor to IFunction

fluids: Project grid anisotropy tensor

fluids: Perform velocity gradient projection

fluids: Setup velocity gradient projection

fluids: Create DM for velocity gradient projection

fluids: Read in files for SGS Data-Driven Model

fluids: Move PHASTA file functions to misc.c

fluids: Refactor matops -> petsc_ops

fluids: Add VecCopyP2C function

IWYU fixes (#1182)

* iwyu - include fixes

* iwyu - silence some iwyu output

* minor - clearer macro names

* iwyu - fix suggestion of "ceed/ceed.h" externally

* iwyu - lighter petsc heade

IWYU fixes (#1182)

* iwyu - include fixes

* iwyu - silence some iwyu output

* minor - clearer macro names

* iwyu - fix suggestion of "ceed/ceed.h" externally

* iwyu - lighter petsc headers

* iwyu - ceed/ceed.h -> ceed.h

* iwyu - cuda/hip include fixes

show more ...

Ken/idl (#1138)

* inflow sponge pushing p towards P0 compiles and runs but unknown yet whether it is effective or not

* looks like my first attempt at this is missing something as there is neglig

Ken/idl (#1138)

* inflow sponge pushing p towards P0 compiles and runs but unknown yet whether it is effective or not

* looks like my first attempt at this is missing something as there is negligible effect from Jacobian addition

* examples/fluids: IDL fixes and demo

* Fix IDL sign error and rename variables/clang-format
* Let newtonianwave.yaml use -reference_pressure and temperature
* Propagate P0 = reference.pressure

* Small tweak to add IDL to the 3 ramp parameters.

* fluids: Extract fluids binary vector load function

- Needs to be used in multiple places (for ICs and for testing data)

* fluids: Correct typo in previous commits

* fluids: Add IDL and newtonian wave test

* fluids: Rename Newtonian Wave --> Gaussian Wave

* fluids: Move RampCoefficient into utils

- Also change name to LinearRampCoefficient, as we may have different
kinds of ramps in the future

* fluids: Change IDL commandline options

* docs(fluids): Document internal damping layer

Co-authored-by: Jed Brown <jed@jedbrown.org>

* fluids: Add ScaleN QF utility

* fluids: Extract InternalDampingLayer QFunction helper

---------

Co-authored-by: Jed Brown <jed@jedbrown.org>
Co-authored-by: KennethEJansen <jansenke@colorado.edu>
Co-authored-by: James Wright <james@jameswright.xyz>

show more ...

fluids: Use ContextGetDouble for time and dt label updates

Merge pull request #1174 from CEED/jeremy/operator-context

CeedOperatorContext* -> CeedOperator*Context*

style - rename CeedOperatorContext functions for consistency

fluids - fix QFunctionContext access errors

Fix for cubic elements in surface forces for lift/drag coefficients (#1158)

* update for surface forces to return dof

* update surface forces

* update on surface forces

* PetscCall section

Fix for cubic elements in surface forces for lift/drag coefficients (#1158)

* update for surface forces to return dof

* update surface forces

* update on surface forces

* PetscCall section get field components

* Update examples/fluids/src/setupts.c

Co-authored-by: Jed Brown <jed@jedbrown.org>

* Update examples/fluids/src/setupts.c

Co-authored-by: Jed Brown <jed@jedbrown.org>

* examples/fluids: lift query from inner loop

* examples/fluids: enable coarser meshes

* update yaml file

* fix path to creating figure

---------

Co-authored-by: Jed Brown <jed@jedbrown.org>

show more ...

fluids: Local vector option for matops

- Passing NULL to the global vector parameter means that the
GlobalToLocal/LocalToGlobal operation is skipped and the context's
local PETSc Vec is used dir

fluids: Local vector option for matops

- Passing NULL to the global vector parameter means that the
GlobalToLocal/LocalToGlobal operation is skipped and the context's
local PETSc Vec is used directly

show more ...

fluids: Remove unused spanstats lumped mass matrix

fluids: Add ComputeL2Projection function

fluids: Add VecP2C, VecC2P, and friends