From 58180b9b54b8705bd089a49def078d6cd6b41269 Mon Sep 17 00:00:00 2001 From: reisswig Date: Fri, 2 Nov 2007 17:57:30 +0000 Subject: removed...now in setup.cc git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/SphericalSurface/trunk@30 40f6ab95-0e4f-0410-8daa-ee8d7420be1d --- src/setup.c | 122 ------------------------------------------------------------ 1 file changed, 122 deletions(-) delete mode 100644 src/setup.c diff --git a/src/setup.c b/src/setup.c deleted file mode 100644 index d4cd599..0000000 --- a/src/setup.c +++ /dev/null @@ -1,122 +0,0 @@ -/* $Header$ */ - -#include -#include - -#include "cctk.h" -#include "cctk_Arguments.h" -#include "cctk_Parameters.h" - - - -void SphericalSurface_Setup (CCTK_ARGUMENTS) -{ - DECLARE_CCTK_ARGUMENTS; - DECLARE_CCTK_PARAMETERS; - - CCTK_REAL const pi = 3.1415926535897932384626433832795028841971693993751; - - int group; - cGroup groupinfo; - cGroupDynamicData groupdata; - - int n; - int ierr; - - - - if (nsurfaces == 0) return; - - - - group = CCTK_GroupIndex ("SphericalSurface::sf_radius"); - assert (group>=0); - ierr = CCTK_GroupData (group, &groupinfo); - assert (!ierr); - ierr = CCTK_GroupDynamicData (cctkGH, group, &groupdata); - assert (!ierr); - - - - for (n=0; n= ntheta[n]); - assert (groupdata.gsh[1] >= nphi[n]); - - assert (ntheta[n] >= 3*nghoststheta[n] && ntheta[n] <= maxntheta); - assert (nphi[n] >= 3*nghostsphi[n] && nphi[n] <= maxnphi); - - - - /* copy parameters into grid functions */ - sf_ntheta[n] = ntheta[n]; - sf_nphi[n] = nphi[n]; - sf_nghoststheta[n] = nghoststheta[n]; - sf_nghostsphi[n] = nghostsphi[n]; - - - - /* coordinates in the theta direction */ - /* avoid_sf_origin_theta = 1 */ - if (symmetric_z[n]) { - - /* upper hemisphere: z>=0, theta in (0, pi/2) */ - sf_delta_theta[n] = pi/2 / (ntheta[n] - 2*nghoststheta[n] - 0.5); - sf_origin_theta[n] = - (nghoststheta[n] - 0.5) * sf_delta_theta[n]; - - } else { - - /* both hemispheres: theta in (0, pi) */ - sf_delta_theta[n] = pi / (ntheta[n] - 2*nghoststheta[n]); - sf_origin_theta[n] = - (nghoststheta[n] - 0.5) * sf_delta_theta[n]; - - } - - - - /* coordinates in the phi direction */ - /* avoid_sf_origin_phi = 0 */ - if (symmetric_x[n]) { - if (symmetric_y[n]) { - - /* one quadrant: x>=0, y>=0, phi in [0, pi/2] */ - assert (nphi[n] - 2*nghostsphi[n] >= 1); - sf_delta_phi[n] = pi/2 / (nphi[n] - 2*nghostsphi[n] - 1); - sf_origin_phi[n] = - nghostsphi[n] * sf_delta_phi[n]; - - } else { - - /* two quadrants: x>=0, phi in [-pi/2, pi/2] */ - assert (nphi[n] - 2*nghostsphi[n] >= 2); - sf_delta_phi[n] = pi / (nphi[n] - 2*nghostsphi[n] - 1); - sf_origin_phi[n] = - pi/2 - nghostsphi[n] * sf_delta_phi[n]; - - } - } else { - if (symmetric_y[n]) { - - /* two quadrants: y>=0, phi in [0, pi] */ - assert (nphi[n] - 2*nghostsphi[n] >= 2); - sf_delta_phi[n] = pi / (nphi[n] - 2*nghostsphi[n] - 1); - sf_origin_phi[n] = - nghostsphi[n] * sf_delta_phi[n]; - - } else { - - /* all quadrants: phi in [0, 2pi) */ - assert (nphi[n] - 2*nghostsphi[n] >= 4); - sf_delta_phi[n] = 2*pi / (nphi[n] - 2*nghostsphi[n]); - sf_origin_phi[n] = - nghostsphi[n] * sf_delta_phi[n]; - - } - } - - - - /* mark surface as uninitialised */ - sf_active[n] = 0; - sf_valid[n] = 0; - - } /* for n */ -} -- cgit v1.2.3