aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
* Improve OpenMP parallelisationHEADmastereschnett2012-09-13
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@136 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Rewrite the interior part of the 8-4 Kreiss-Oliger dissipation operatordiener2012-05-14
| | | | | | | | | | | to use an explicit loop instead of Fortran 90 array statements since the openmp workshare directive is not parallelized by the intel compilers. This is based on an original patch by Christian Reisswig. Do for this is only done for this operator since this is most often used in BBH simulations. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@132 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Declare all loop variables as OpenMP privateeschnett2012-05-05
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@131 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Correct grid extent in call to deriv2_mixedeschnett2011-05-11
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@127 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Add a routine to return the size of the boundary region, i.e. the regiondiener2010-03-19
| | | | | | | | where centered finite differences are not performed. This is necessary in order to mix SBP and hard coded stencils in CTGamma. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@120 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Introduce a boolean parameter to switch of SBP near outer boundary. Insteaddiener2009-03-31
| | | | | | | | switch to progressively lower order centered stencil as the boundery is approached, similar to what is done with the second order and upwind stencils. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@118 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Added functionality to set up public grid functions (sbp_dx, sbp_dy, sbp_dz)diener2008-07-14
| | | | | | | | | | | containing the delta's used in the Kreiss-Oliger type dissipation. This is activated with a parameter (use_variable_deltas) and a routine (SBP_DeltaInitial) is scheduled at BASEGRID to set the default values initially. The grid functions can be modified later by a scheduled user function if so desired. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@116 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Make these routines aware of the offset applied to the underlying patchdiener2008-07-14
| | | | | | | | structure so that integration done with these coefficients will leave out the overlap points. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@115 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Add non-sbp upwinding operators. These have the same order as the centereddiener2008-05-08
| | | | | | | operators but does not satisfy any SBP property. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@114 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* (*) Small modifications to use GetBoundarySpecification orkorobkin2008-05-06
| | | | | | | | MultiPatch_GetBoundarySpecification instead of reading shiftout parameters directly from the CoordBase thorn. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@113 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* (+) Poisoning.F90korobkin2008-05-02
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@112 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* (+) implemented poisoning: added three parameterskorobkin2008-05-02
| | | | | | | | | | | | - poison_derivatives - poison_dissipation - poison_value The first two set values of derivative GF and RHS respectively to poison_value in boundary_shiftout zone, which is specified by CoordBase shiftout parameters. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@111 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* (-) removed debugging outputkorobkin2008-05-01
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@110 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* (*) CCTK_FNAME(get_grid_offsets) -> get_grid_offsetskorobkin2008-05-01
| | | | | | | | | | (+) Added get_grid_offsets to thorn interface as GetFDGridOffsets (+) Added another thorn interface function GetLSHIndexRanges to inquire actual local index ranges (in Fortran notation), which account for both the ghostzones and shiftout values git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@109 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Fixed use of wrong offset in one of the y-directions.diener2008-04-30
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@108 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* (+) get_offset.ckorobkin2008-04-29
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@107 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* (*) The grid offsets of the dissipation and derivative operators are madekorobkin2008-04-29
| | | | | | | | | correlated. If use_shiftout = yes, both offsets are read from the CoordBase::boundary_shiftout_* parameters. Otherwise, both offsets are set to zero. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@106 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* (+) Added the capability to specify grid offsets in derivative and dissipationkorobkin2008-04-28
| | | | | | | | operators, as set by boundary_shiftout* parameters of CoordBase. Only first-order centered derivative operators are affected. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@105 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Add stencils for non-SBP operators.diener2008-04-18
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@104 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Fix wrong initialization of imax.diener2008-04-18
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@103 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Determine correctly which boundaries are outer boundaries whenschnetter2008-04-17
| | | | | | | off-centring the stencil, if multi-patch is used git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@102 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Add OpenMP workshare constructsschnetter2008-04-17
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@101 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Parallelise with OpenMP. Remove function outer_prod_3d, ad write itsschnetter2008-04-08
| | | | | | | | content explicitly, since this is much easier to parallelise (and also shorter). git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@99 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Determine index range correctly: Initialise imax to -1 instead of 0,schnetter2008-04-08
| | | | | | | | | so that empty domains are handled correctly. Begin to add some OpenMP parallelisation statements. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@98 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Remove unused local variablesschnetter2008-04-08
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@97 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Add two new parameters onesided_outer_boundaries andschnetter2008-04-08
| | | | | | | | | onesided_interpatch_boundaries. These decide whether the user wants one-sided differencing stencil near outer and near multi-patch boundaries. Both default to yes for backward compatibility. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@96 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Fix wrong subroutine call.diener2007-10-31
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@94 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Add 8-4 second derivative operators.diener2007-10-31
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@93 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* Fix a couple of typos.diener2007-10-30
| | | | git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@92 f69c4107-0314-4c4f-9ad4-17e986b73f4a
* 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@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