aboutsummaryrefslogtreecommitdiff
path: root/src/dissipation.c
Commit message (Collapse)AuthorAge
* Two part update of the restricted full norm dissipation operators. The firstdiener2006-03-10
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | 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
* Only use one sided finite differences near multipatch boundaries. Should fixdiener2006-02-13
| | | | | | | | | 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
* Kreiss-Oliger type dissipation for the restricted full norm 4-3 operator.diener2006-02-03
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@57 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Kreiss-Oliger type dissipation operators compatible with the diagonal 4-2 anddiener2006-02-02
| | | | | | | 6-3 derivative operators. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@55 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Implementation of a Kreiss-Oliger type dissipation operator compatible withdiener2006-02-01
| | | | | | | 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
* Introduce a new kind of dissipation operators. Kreiss-Oliger type dissipationdiener2005-12-15
| | | | | | | | | | | 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
* Added an optimized 6-3 operator and dissipation operators for the remainingdiener2005-06-15
| | | | | | | diagonal norm operators. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@38 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Correct prototype for Fortran routine.schnetter2005-06-04
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@36 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Dissipation operator for the optimized 6-5 restricted full norm operator.diener2005-05-26
| | | | | | | 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
* Dissipation operator for the optimized 4-3 restricted full norm case. Testingdiener2005-05-25
| | | | | | | 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
* Optimised 8-4 derivative operators. Experimental 8-4 dissipation operators.diener2005-05-19
| | | | | | | 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
* Implementation of 6th order dissipation operator compatible with the 6-5diener2005-02-08
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