| Commit message (Collapse) | Author | Age |
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@130 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
| |
where centered finite differences are not performed. This is necessary
in order to mix SBP and hard coded stencils in CTGamma.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@120 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
|
|
|
|
|
|
|
|
| |
MultiPatch_GetBoundarySpecification instead of reading shiftout
parameters directly from the CoordBase thorn.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@113 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
|
| |
(+) Added get_grid_offsets to thorn interface as GetFDGridOffsets
(+) Added another thorn interface function GetLSHIndexRanges
to inquire actual local index ranges (in Fortran notation),
which account for both the ghostzones and shiftout values
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@109 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@100 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
| |
FourierDifferencing, providing a routine to return mixed derivatives obtained
by applying two succesive 1st derivatives in different directions. Not
tested yet, but compiles.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@87 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
|
|
| |
them in one place. Added routines to return coefficients in all the
missing cases. Also added routines to return 2 derivatives (only 2-1 and
4-2 so far). The routines calculating derivatives for whole grid functions
have not been tested yet, so if you are using them, you will probably want
to hold off on updating until I give the clear signal.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@85 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
| |
Initialise normmask in basegrid instead of initial. Mark it so that
it is not checkpointed.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@68 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
| |
direction for the upwinding derivatives. The directions are specified in a
real 3d array with the sign choosing the upwinding direction.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@64 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@63 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
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@60 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@27 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
| |
request from Bela.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@23 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@14 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
MultiPatchWaveToy to SummationByParts.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@12 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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
(after Erik explained to me that CCTK_ARGUMENTS can be accessed even though
they are not passed into the routine) that is called with the variables
themselves and not their names. These are provided by the aliased function
Diff_gf, so now the x-derivatives of rho can be found with the call:
Diff_gv ( cctkGH, 0, rho, drhox );
Also made the determination of the gridspacing use CCTK_DELTA_SPACE (now that I know how to do it) so the routines should now also work with mesh-refinement.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@4 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@2 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|