aboutsummaryrefslogtreecommitdiff
Commit message (Collapse)AuthorAge
...
* Call SBP_determine_onesided_stencil instead of determining theschnetter2006-10-05
| | | | | | | | boundary type independently. This may correct an error on multiple processors. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@78 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Simplify some expressions calculating dissipation coefficients. Thisschnetter2006-08-26
| | | | | | | | | | makes it compile with certain versions of gfortran. Change type of argument "part" to set_dmatrix_half from CCTK_INT to integer since it is called with integer (not CCTK_INT) values. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@77 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Correct Fortran wrapperschnetter2006-08-03
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@76 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Introduce a function SBP_determine_onesided_stencil which determinesschnetter2006-07-05
| | | | | | | | | | | | which faces should use one-sided stencils, depending on which boundaries are inter-processor boundaries, symmetry boundaries, and multi-patch boundaries. Use this function everywhere. Remove the previous mechinisms; some were not in all cases correct. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@75 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Replace numbers of the form "0_wp" with "0.0_wp". The former definesschnetter2006-06-11
| | | | | | | integer constants, whereas real constants are intended. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@74 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Replace the temporary solution, which checks whether PUGH is active,schnetter2006-06-11
| | | | | | | | with a better way to determine which boundaries are outer boundaries, symmetry boundaries, interprocessor boundaries, or penalty boundaries. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@73 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* *** empty log message ***tiglio2006-05-18
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@72 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* *** empty log message ***tiglio2006-05-18
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@71 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* A C version (as opposed to C++) of the hack.tiglio2006-05-18
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@70 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* 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
* Do not inherit from grid. This is not necessary.schnetter2006-05-13
| | | | | | | | 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
* 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
* Declare the aliased function to call the upwinding operators.diener2006-04-10
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@63 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
* Changed the default operator type to "Optimized". Added in the comment thatdiener2005-11-29
| | | | | | | 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
* 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
* Added a boolean parameter "check_grid_sizes" to turn off checking of grid sizesdiener2005-06-10
| | | | | | | | | | 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
* 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
* Additional documentation. Update with the table_handle interface.diener2005-03-01
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@30 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
generated by cgit v1.2.3 (git 2.39.1) at 2024-05-20 00:51:25 +0000