aboutsummaryrefslogtreecommitdiff
Commit message (Collapse)AuthorAge
* Move check for domain size to the right place.diener2007-10-30
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@91 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* This operator is now implemented correctly.diener2007-10-30
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@90 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Added 6-3 minimal bandwidth 2nd derivative operators. Added 4-2 minimaldiener2007-10-19
| | | | | | | | bandwidth 2nd derivative operators. Moved the 4-2 optimised operators. Fixed missing break statements get_up_coeffs.c. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@89 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Fix error message.diener2007-10-18
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@88 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Make 2nd derivative interface consistent with FiniteDifferencing anddiener2007-10-18
| | | | | | | | | 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
* The second derivative in the 4-2 case is not unique, so I added the minimaldiener2007-10-18
| | | | | | | bandwidth operator in addition to the ABTE operator already there. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@86 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Moved definitions of coefficients to a single module in order to only havediener2007-10-16
| | | | | | | | | | | 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
* Add a deallocate statement for the local allocatable arrays in case wediener2007-05-14
| | | | | | | encounter a F90 only compiler. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@84 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Require 2, not 1 ghost zones.schnetter2007-02-12
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@83 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Do not abort during paramcheck if there are not enough ghost zones;schnetter2007-02-12
| | | | | | | print a warning only. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@82 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Check that the boundary regions (where the stencils need to beschnetter2006-10-05
| | | | | | | | | one-sided) do not overlap. Check the number of ghost zones. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@81 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Correct one-off error in the check that the boundary regions do notschnetter2006-10-05
| | | | | | | | | overlap. Check the number of ghost zones as well. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@80 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Ensure that the boundary regions where the stencil is off-centered doschnetter2006-10-05
| | | | | | | not overlap. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@79 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* 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
generated by cgit v1.2.3 (git 2.39.1) at 2024-06-30 04:25:16 +0000