aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
...
* Added a small ifdef to make it work with PUGH (as it was the derivs were ↵tiglio2006-05-18
| | | | | | | | | | 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
* Fixed stupid mistake in having the assignment of stored coefficients anddiener2006-04-21
| | | | | | | 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
* Fixed bug when running on one processor with ghost size less than 3.diener2006-04-21
| | | | | | | | 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
* Another file I forgot previously.diener2006-04-13
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@65 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Modify the interface to allow for spatial dependence on the upwindingdiener2006-04-13
| | | | | | | | 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
* First attempt at upwinding SBP operators. Don't use yet. They will changediener2006-04-10
| | | | | | | soon to be more flexible. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@62 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* 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
* Aliased routine for returning the diagonal norm coefficients.diener2006-02-24
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@60 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Add new parameters derivatives_zero_y and derivatives_zero_z in order todiener2006-02-14
| | | | | | | | | | | | | | 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
* 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
* Cosmetic change.diener2006-02-03
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@56 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 a new boolean parameter "scale_with_h" that allows the user todiener2005-12-06
| | | | | | | | 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
* Remove reference to unused file.diener2005-12-06
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@51 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Remove unused file.diener2005-12-06
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@50 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Call the same 2-1 and 4-2 operators regardless of whether the user requestsdiener2005-12-02
| | | | | | | "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
* Remove declaration of unused variables nghostzones, bbox, lsh and ierr.diener2005-11-29
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@48 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Remove declaration of unused variables nwt and xmax.diener2005-11-29
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@47 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Remove declaration of unused variables i,j,k.diener2005-11-29
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@46 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Remove commented debug statements.diener2005-11-29
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@44 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Remove commented lines.diener2005-11-29
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@43 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Realized that in the CVS version the optimized restricted full operatorsdiener2005-11-29
| | | | | | | | | | 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
* Remove some commented debugging print statements.diener2005-11-29
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@41 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Remove an unused subroutine.diener2005-11-29
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@40 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Clean up for publication.diener2005-11-29
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@39 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
* Fixed a typo/bug that for some reason didn't cause any problems yet (evendiener2005-05-26
| | | | | | | 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
* 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
* New 4-3 and 6-5 operators that are optimised for accuracy. Added a newdiener2005-05-12
| | | | | | | 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
* Analyse the options table more gracefully, and abort with a nice errorschnetter2005-02-27
| | | | | | | message if something is wrong. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@29 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Analyse the options table more gracefully, and abort with a nice errorschnetter2005-02-27
| | | | | | | message if something is wrong. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@28 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Added a table handle to the derivative interface.diener2005-02-25
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@27 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Forgot to commit this.diener2005-02-24
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@26 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Routines for getting the coefficients for the various derivatives. By ↵diener2005-02-23
| | | | | | request from Bela. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@23 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
* Remove unused routine SetNMask.schnetter2005-01-31
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@20 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Use correct routine name "SBP_SetNormMask" instead ofschnetter2005-01-31
| | | | | | | "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
* I forgot the actual finite difference routine the first time.diener2005-01-07
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@18 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* 6th order interior, 5th order at the boundary finite differences (restricteddiener2005-01-07
| | | | | | | | 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
* Routine to check if grid sizes and ghost_zones are enough for the order ofdiener2004-12-01
| | | | | | | 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
* Make declaration of Diff_gv and Diff_gf const-correct.schnetter2004-11-10
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@15 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Move scalar product weight mask and its calculation fromschnetter2004-10-25
| | | | | | | MultiPatchWaveToy to SummationByParts. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@12 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Pass Fortran arguments by referenceschnetter2004-10-22
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@11 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Added 4-3 Full restricted norm operators. Added parameter to choose the typediener2004-10-11
| | | | | | | | | 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