diff options
| author | diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c> | 2002-05-23 15:27:10 +0000 |
|---|---|---|
| committer | diener <diener@2a26948c-0e4f-0410-aee8-f1d3e353619c> | 2002-05-23 15:27:10 +0000 |
| commit | cc6fff6b1a0158cc10ddc5e66e0003921b5fddee (patch) | |
| tree | 922aeb8ca09a985959b588bf5858287a2a710177 /src/include | |
| parent | 11b26b0e33a56e1dbc45757226176f1b1472321f (diff) | |
Added updates of logical varibles storing symmetries. Might not be
necessary anyway. The jury is still out.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/EHFinder/trunk@23 2a26948c-0e4f-0410-aee8-f1d3e353619c
Diffstat (limited to 'src/include')
| -rw-r--r-- | src/include/physical_part.h | 62 |
1 files changed, 50 insertions, 12 deletions
diff --git a/src/include/physical_part.h b/src/include/physical_part.h index 44166ea..c14439a 100644 --- a/src/include/physical_part.h +++ b/src/include/physical_part.h @@ -15,6 +15,8 @@ ngz = cctk_nghostzones(3) ! ixl and ixr is set to the minimum and maximum local cell index in the ! x-direction. The corresponding is done for the y- and z-directions. +symx = .false.; symy = .false.; symz = .false.; + ixl = 1 ixr = nx jyl = 1 @@ -40,9 +42,18 @@ if ( cctk_bbox(6) .eq. 0 ) kzr = kzr - ngz ! minimum and maximum cell index appropriately with the size of the ghostzone. if ( CCTK_Equals(domain,'octant') ) then - if (cctk_lbnd(1) .eq. 0 ) ixl = ixl + ngx - if (cctk_lbnd(2) .eq. 0 ) jyl = jyl + ngy - if (cctk_lbnd(3) .eq. 0 ) kzl = kzl + ngz + if (cctk_lbnd(1) .eq. 0 ) then + ixl = ixl + ngx + symx = .true. + end if + if (cctk_lbnd(2) .eq. 0 ) then + jyl = jyl + ngy + symy = .true. + end if + if (cctk_lbnd(3) .eq. 0 ) then + kzl = kzl + ngz + symz = .true. + end if end if ! In quadrant I do the same, except here I have to check in which direction @@ -50,16 +61,34 @@ end if if ( CCTK_Equals(domain,'quadrant') ) then if ( CCTK_Equals(quadrant_direction,'x') ) then - if (cctk_lbnd(2) .eq. 0 ) jyl = jyl + ngy - if (cctk_lbnd(3) .eq. 0 ) kzl = kzl + ngz + if (cctk_lbnd(2) .eq. 0 ) then + jyl = jyl + ngy + symy = .true. + end if + if (cctk_lbnd(3) .eq. 0 ) then + kzl = kzl + ngz + symz = .true. + end if end if if ( CCTK_Equals(quadrant_direction,'y') ) then - if (cctk_lbnd(1) .eq. 0 ) ixl = ixl + ngx - if (cctk_lbnd(3) .eq. 0 ) kzl = kzl + ngz + if (cctk_lbnd(1) .eq. 0 ) then + ixl = ixl + ngx + symx = .true. + end if + if (cctk_lbnd(3) .eq. 0 ) then + kzl = kzl + ngz + symz = .true. + end if end if if ( CCTK_Equals(quadrant_direction,'z') ) then - if (cctk_lbnd(1) .eq. 0 ) ixl = ixl + ngx - if (cctk_lbnd(2) .eq. 0 ) jyl = jyl + ngy + if (cctk_lbnd(1) .eq. 0 ) then + ixl = ixl + ngx + symx = .true. + end if + if (cctk_lbnd(2) .eq. 0 ) then + jyl = jyl + ngy + symy = .true. + end if end if end if @@ -67,12 +96,21 @@ end if if ( CCTK_Equals(domain,'bitant') ) then if ( CCTK_Equals(bitant_plane,'xy') ) then - if (cctk_lbnd(3) .eq. 0 ) kzl = kzl + ngz + if (cctk_lbnd(3) .eq. 0 ) then + kzl = kzl + ngz + symz = .true. + end if endif if ( CCTK_Equals(bitant_plane,'xz') ) then - if (cctk_lbnd(2) .eq. 0 ) jyl = jyl + ngy + if (cctk_lbnd(2) .eq. 0 ) then + jyl = jyl + ngy + symy = .true. + end if endif if ( CCTK_Equals(bitant_plane,'yz') ) then - if (cctk_lbnd(1) .eq. 0 ) ixl = ixl + ngx + if (cctk_lbnd(1) .eq. 0 ) then + ixl = ixl + ngx + symx = .true. + end if endif end if |