| Commit message (Collapse) | Author | Age |
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@134 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@119 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
| |
switch to progressively lower order centered stencil as the boundery is
approached, similar to what is done with the second order and upwind stencils.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@118 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@117 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
|
|
| |
containing the delta's used in the Kreiss-Oliger type dissipation. This is
activated with a parameter (use_variable_deltas) and a routine
(SBP_DeltaInitial) is scheduled at BASEGRID to set the default values
initially. The grid functions can be modified later by a scheduled user
function if so desired.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@116 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
operators but does not satisfy any SBP property.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@114 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
|
|
|
| |
- poison_derivatives
- poison_dissipation
- poison_value
The first two set values of derivative GF and RHS respectively
to poison_value in boundary_shiftout zone, which is specified by
CoordBase shiftout parameters.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@111 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@104 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
| |
onesided_interpatch_boundaries. These decide whether the user wants
one-sided differencing stencil near outer and near multi-patch
boundaries. Both default to yes for backward compatibility.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@96 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
part is an optimisation. Previously the coefficients were constructed every
time the dissipation was applied. This turned out to have a significant
overhead, especially for large domains. Now the coefficients are constructed
only once for each patch on each processor and stored. The thorn now uses
the aliased functions MultiPatch_GetMap and MultiPatch_GetMaps (provided by
thorn MultiPatch) to figure out how many patches are present and which patch
is currently being worked on. These dissipation operators will not work unless
the MultiPatch thorn is used. Right now this should not be a problem (but can
be changed if it is), unless somebody uses dissipation operators outside of
the multipatch infrastructure.
The second part is an update aimed at making the restricted full norm
dissipation operators better behaved on grids with very different number of
grid points in different directions. For that reason the parameter:
diss_fraction
has been made vector valued (of length 3), so that the width of the transition
region from the boundary to the interior operators can be controlled
separately in the 3 dimensions. Unfortunately this means that any parameter
files, that explicitly sets this parameter, will have to be modified. I don't
think too many people are using the restricted full dissipation operators yet,
so it shouldn't cause too many problems. In additiona new parameter:
h_scaling
has been introduced. This parameter is also vector valued of length 3 and
allows the local grid spacing to be rescaled to more closely resemble the
physical grid spacing in any dimension.
Note, these changes only affect the non-diagonal dissipation operators.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@61 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
ecplictly setting any derivatives in the y- and/or z-direction to zero
when simulating 1- or 2-D problems with a 3D code. This avoids issues with
errors at roundoff level in the y- and z-directions. Should be used with care
and together with the parameter check_grid_sizes to avoid complaints about
too small grids in the y- and z-directions.
Has only been tested in a few cases. If problems are caused by this commit,
please report them immediately.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@59 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
|
|
| |
operator suitably modified near the boundary. Since the dissipation operator
in the interior uses a stencil that is one point wider than the derivative
operator, an additional check on the number of ghostzones in multi processor
runs have been added. Only implemented for the 8-4 and optimized 6-5
operators.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@53 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
| |
choose to scale the dissipation with the grid spacing h. This results
formally in the loss of one order of convergence. The default value is "no".
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@52 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
8-4 is minimal spectral radius instead of minimal bandwidth.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@45 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
|
| |
and ghost zones. This is useful for Burkhard, when he wants to use the
operators in lower dimensional cases. It could be done in a more selective
way, so the checks were done where it makes sense, but this will have to do
for now.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@37 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
still in progress, but it looks like the convergence order is right.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@33 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
parameter to choose between minimim bandwidth and optimized operators.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@31 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
finite difference operators. Not fully tested yet. The MoL dissipation
interface was shamelessly stolen from Eriks Dissipation thorn (in AEIThorns).
Things to do:
Add a test for the presence of MoLQueryEvolvedRHS (since the use of
dissipation is optional, I didn't want to require the presence of
this routine).
Add compatible dissipation operators to the other finite difference operators.
More testing.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@21 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@13 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
| |
of differences (diagonal, full restricted). Added routine to return the
coefficient of the norm at the boundary. Added routine to return the
mask for the norm calculation.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@10 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
has to be modified in other ways as well.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@9 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@6 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@2 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|