| Commit message (Collapse) | Author | Age |
... | |
|
|
|
|
|
|
|
|
|
| |
wrong near real boundaries if
using PUGH). One can do something better but I did not want to mess up the MP case.
It will not work with reflection symmetries and PUGH as it is, though, so this is a temporary solution only.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@69 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
deallocation of the temporary array within the wrong if statement.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@67 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
| |
Also fix stupid mistake in having the assignment of stored coefficients and
deallocation of the temporary array within the wrong if statement.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@66 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@65 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
|
|
|
|
|
|
|
| |
soon to be more flexible.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@62 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
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
| |
the problems we have with periodic boundary conditions. Also added the missing
declarations of fortran routines with Kreiss-Oliger type dissipation for
certain cases.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@58 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@57 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@56 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
6-3 derivative operators.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@55 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
the diagonal 2-1 derivative operators.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@54 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
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@51 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@50 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
"Minimal Bandwidth" or "Optimized" operators, since they are unique.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@49 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@48 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@47 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@46 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@44 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@43 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
|
|
| |
were using the non-optimized values for the norm at the boundary (the
element of the norm used in the penalty term). The results in the SBP paper
were all computed with versions of the code using the right values. It
shouldn't make much difference anyway.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@42 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@41 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@40 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@39 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
diagonal norm operators.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@38 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@36 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
Don't use yet. Still being tested.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@35 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
though I think it should have).
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@34 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
|
|
|
|
|
|
|
| |
Tests of these are on going.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@32 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
|
|
|
|
|
|
|
| |
message if something is wrong.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@29 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
message if something is wrong.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@28 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@27 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@26 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@20 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
"sbp_set_norm_mask" (which has too many underscores).
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@19 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@18 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
|
| |
full norm). Thanks to José María Martín García for the mathematica notebook
that provided the coefficients. Not yot tested extensively, but seems to work.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@17 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
|
|
|
| |
finite differencing requires. Has not been tested thoroughly.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@16 f69c4107-0314-4c4f-9ad4-17e986b73f4a
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@15 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
|
|
|
|
| |
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@11 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
|