Commits

commits@lists.einsteintoolkit.org

2 participants
6997 discussions

[svn:einsteintoolkit] IDAxiOddBrillBH/trunk/src/ (Rev. 70)
by schnetter＠cct.lsu.edu 30 Oct '13

30 Oct '13

User: eschnett Date: 2013/10/30 10:19 AM Modified: /trunk/src/ IDAxiOddBrillBH.F, ckerr.x, gij.x, kerr.x, ksq_bs.x, qfunc_even.x Log: Change real*8 to CCTK_REAL File Changes: Directory: /trunk/src/ ====================== File [modified]: IDAxiOddBrillBH.F Delta lines: +11 -11 =================================================================== --- trunk/src/IDAxiOddBrillBH.F 2013-05-27 22:32:45 UTC (rev 69) +++ trunk/src/IDAxiOddBrillBH.F 2013-10-30 15:19:25 UTC (rev 70) @@ -34,26 +34,26 @@ DECLARE_CCTK_PARAMETERS DECLARE_CCTK_FUNCTIONS - real*8 deta,dq - real*8, allocatable :: ac(:,:),ae(:,:),aw(:,:),an(:,:),as(:,:), + CCTK_REAL deta,dq + CCTK_REAL, allocatable :: ac(:,:),ae(:,:),aw(:,:),an(:,:),as(:,:), $ rhs(:,:),qfetaeta(:,:),qfqq(:,:),ksq(:,:), $ psi2dv(:,:),dpsi2dv(:,:),ddpsi2dv(:,:), $ detapsisph(:,:),dqpsisph(:,:),detaetapsisph(:,:), $ detaqpsisph(:,:),dqqpsisph(:,:), $ etagrd(:),qgrd(:) - real*8 o1,o2,o3,o4,o5,o6,o7,o8,o9,o10, + CCTK_REAL o1,o2,o3,o4,o5,o6,o7,o8,o9,o10, $ o11,o12,o13,o14,o15,o16,o17,o18,o19, $ o20,o21,o22,o23,o24,o25,o26,o27,o28,o29, $ o30,o31,o32,o33,o34,o35,o36,o37,o38,o39, $ o40,o41,o42,o43,o44,o45,o46,o47,o48,o49, $ o50,o51 - real*8 t1,t2,t3,t4 - real*8 gtil,dngtil,dnngtil,dnnngtil,dnnnngtil,dnnnnngtil, + CCTK_REAL t1,t2,t3,t4 + CCTK_REAL gtil,dngtil,dnngtil,dnnngtil,dnnnngtil,dnnnnngtil, $ dnnnnnngtil,dnnnnnnngtil - real*8 rhsmax,rmax,odd_get2d,adm,Jmom,mass,a,r_iso,rBL,rbh_eps - real*8,parameter :: rbh_tol = 1.0d-7 + CCTK_REAL rhsmax,rmax,odd_get2d,adm,Jmom,mass,a,r_iso,rBL,rbh_eps + CCTK_REAL,parameter :: rbh_tol = 1.0d-7 integer,parameter :: itmax = 100 - real*8 pi,zero,one,two + CCTK_REAL pi,zero,one,two CCTK_REAL :: xmin,xmax,ymin,ymax,zmin,zmax integer :: ne, nq integer :: nx,ny,nz @@ -419,10 +419,10 @@ c c Compute eta,q,phi at the each points of cartesian grid - eta = 0.5d0 * dlog(x**2 + y**2 + z**2) + eta = 0.5 * log(x**2 + y**2 + z**2) abseta = abs (eta) - q = datan2 (sqrt (x**2 + y**2),z) - phi = datan2 (y, x) + q = atan2 (sqrt (x**2 + y**2),z) + phi = atan2 (y, x) do k = 1,nz do j = 1,ny File [modified]: ckerr.x Delta lines: +3 -3 =================================================================== --- trunk/src/ckerr.x 2013-05-27 22:32:45 UTC (rev 69) +++ trunk/src/ckerr.x 2013-10-30 15:19:25 UTC (rev 70) @@ -1,7 +1,7 @@ - mass = sqrt(2.00000000000000d0 + sqrt(4.0000000000000d0 + byJ**2 + mass = sqrt(2.0 + sqrt(4.0 + byJ**2 & )) a = byJ/mass - rBL = 5.0000000000000d-1*exp(eta(i,j,k))*(1.00000000000000d0 + (1.d0*(a + rBL = 5.0e-1*exp(eta(i,j,k))*(1.0 + (1.0*(a & + mass)*exp(-eta(i,j,k)))/sqrt(-a**2 + mass**2))*sqrt(-a**2 + mass**2) - & *(1.00000000000000d0 + (1.d0*(a + mass)*exp(-eta(i,j,k)))/sqrt(a**2 + m + & *(1.0 + (1.*(a + mass)*exp(-eta(i,j,k)))/sqrt(a**2 + m & ass**2)) File [modified]: gij.x Delta lines: +6 -6 =================================================================== --- trunk/src/gij.x 2013-05-27 22:32:45 UTC (rev 69) +++ trunk/src/gij.x 2013-10-30 15:19:25 UTC (rev 70) @@ -1,9 +1,9 @@ - o1 = 2.d0*phi(i,j,k) + o1 = 2.0*phi(i,j,k) o2 = cos(o1) o3 = -xi o4 = cos(q(i,j,k)) o5 = o4*xi - o6 = 1.d0 + o3 + o5 + o6 = 1.0 + o3 + o5 o7 = -eta0_me o8 = eta(i,j,k) + o7 o9 = o8**2 @@ -18,12 +18,12 @@ o18 = o13 + o17 o19 = sin(q(i,j,k)) o20 = o19**n_me - o21 = 2.d0*amp_me*o18*o20*o6 + o21 = 2.0*amp_me*o18*o20*o6 o22 = exp(o21) - o23 = -1.d0 + o22 - gxx(i,j,k) = 5.d-1*(1.d0 + o22 + o2*o23) + o23 = -1.0 + o22 + gxx(i,j,k) = 5.d-1*(1.0 + o22 + o2*o23) gxy(i,j,k) = o23*cos(phi(i,j,k))*sin(phi(i,j,k)) gxz(i,j,k) = 0 - gyy(i,j,k) = 5.d-1*(1.d0 + o22 - o2*o23) + gyy(i,j,k) = 5.d-1*(1.0 + o22 - o2*o23) gyz(i,j,k) = 0 gzz(i,j,k) = o22 File [modified]: kerr.x Delta lines: +3 -3 =================================================================== --- trunk/src/kerr.x 2013-05-27 22:32:45 UTC (rev 69) +++ trunk/src/kerr.x 2013-10-30 15:19:25 UTC (rev 70) @@ -1,7 +1,7 @@ - mass = sqrt(2.00000000000000d0 + sqrt(4.0000000000000d0 + byJ**2 + mass = sqrt(2.0 + sqrt(4.0 + byJ**2 & )) a = byJ/mass - rBL = 5.0000000000000d-1*exp(etagrd(i))*(1.00000000000000d0 + (1.d0*(a + rBL = 5.0e-1*exp(etagrd(i))*(1.0 + (1.0*(a & + mass)*exp(-etagrd(i)))/sqrt(-a**2 + mass**2))*sqrt(-a**2 + mass**2) - & *(1.00000000000000d0 + (1.d0*(a + mass)*exp(-etagrd(i)))/sqrt(a**2 + m + & *(1.0 + (1.0*(a + mass)*exp(-etagrd(i)))/sqrt(a**2 + m & ass**2)) File [modified]: ksq_bs.x Delta lines: +5 -5 =================================================================== --- trunk/src/ksq_bs.x 2013-05-27 22:32:45 UTC (rev 69) +++ trunk/src/ksq_bs.x 2013-10-30 15:19:25 UTC (rev 70) @@ -1,12 +1,12 @@ o1 = cos(qgrd(j)) - o2 = 1.0d0 + o2 = 1.0 o3 = o2**2 o4 = 1/o3 o5 = sin(qgrd(j)) - o6 = 2.0d0 + n + o6 = 2.0 + n o7 = o5**2 o8 = -(o6*o7) - o9 = 1.0d0 + n + o8 + o9 = 1.0 + n + o8 o10 = o3**2 o11 = o10**2 o12 = dngtil**2 @@ -15,9 +15,9 @@ o15 = o9**2 c exc33(i,j) = 0 c exc32(i,j) = -(dngtil*o1*o4*o5**n) -c exc31(i,j) = gtil*o4*o5**(-1.00000000000000d0 + n)*o9 +c exc31(i,j) = gtil*o4*o5**(-1.0 + n)*o9 c exc22(i,j) = 0 c exc21(i,j) = 0 c exc11(i,j) = 0 - ksq(i,j) = (2.0d0*o5**(2*( + ksq(i,j) = (2.0*o5**(2*( $ -2 + n))*(o14*o15 + o12*o13*o7))/(o10*o11) File [modified]: qfunc_even.x Delta lines: +5 -5 =================================================================== --- trunk/src/qfunc_even.x 2013-05-27 22:32:45 UTC (rev 69) +++ trunk/src/qfunc_even.x 2013-10-30 15:19:25 UTC (rev 70) @@ -1,7 +1,7 @@ o1 = -xi o2 = cos(qgrd(j)) o3 = xi*o2 - o4 = 1.d0 + o1 + o3 + o4 = 1.0 + o1 + o3 o5 = -eta0_me o6 = etagrd(i) + o5 o7 = o6**2 @@ -19,8 +19,8 @@ o19 = o18**n_me o20 = o2**2 o21 = o13 + o17 - qfetaeta(i,j) = amp_me*o19*o4*(-2.d0*o11*o13 - 2.d0*o11*o17 + 4.d0*o10*o15*o - & 17 + 4.d0*o10*o13*o7) - qfqq(i,j) = -(amp_me*xi*n_me*o19*o2*o21) - amp_me*xi*(1.d0 + n_me)*o19*o2*o21 - - & amp_me*n_me*o19*o21*o4 + amp_me*(-1.d0 + n_me)*n_me*o18**(-2 + n)*o20*o21*o + qfetaeta(i,j) = amp_me*o19*o4*(-2.0*o11*o13 - 2.0*o11*o17 + 4.0*o10*o15*o + & 17 + 4.0*o10*o13*o7) + qfqq(i,j) = -(amp_me*xi*n_me*o19*o2*o21) - amp_me*xi*(1.0 + n_me)*o19*o2*o21 - + & amp_me*n_me*o19*o21*o4 + amp_me*(-1.0 + n_me)*n_me*o18**(-2 + n)*o20*o21*o & 4

1 0

[svn:einsteintoolkit] EOSG_IdealFluid/trunk/src/ (Rev. 12)
by schnetter＠cct.lsu.edu 30 Oct '13

30 Oct '13

User: eschnett Date: 2013/10/30 10:18 AM Modified: /trunk/src/ EOS_GIF_Setup.F90 Log: Reduce "large" number File Changes: Directory: /trunk/src/ ====================== File [modified]: EOS_GIF_Setup.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/EOS_GIF_Setup.F90 2011-10-25 18:45:09 UTC (rev 11) +++ trunk/src/EOS_GIF_Setup.F90 2013-10-30 15:18:36 UTC (rev 12) @@ -21,6 +21,6 @@ eos_if_gamma_local = eos_ideal_fluid_gamma mean_molecular_weight_local = mean_molecular_weight - k_boltzmann = 1.d80 + k_boltzmann = 1.d30 end subroutine EOS_GIF_Setup

1 0

[svn:einsteintoolkit] Extract/trunk/src/ (Rev. 127)
by schnetter＠cct.lsu.edu 30 Oct '13

30 Oct '13

User: eschnett Date: 2013/10/30 10:17 AM Modified: /trunk/src/ spin_integrand3D.F Log: Abort on error File Changes: Directory: /trunk/src/ ====================== File [modified]: spin_integrand3D.F Delta lines: +1 -1 =================================================================== --- trunk/src/spin_integrand3D.F 2013-05-27 22:31:52 UTC (rev 126) +++ trunk/src/spin_integrand3D.F 2013-10-30 15:17:55 UTC (rev 127) @@ -60,7 +60,7 @@ ip = 1 CASE DEFAULT - WRITE(*,*) "This value of Psi_power is not supported" + CALL CCTK_ERROR("This value of Psi_power is not supported") END SELECT

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 586)
by schnetter＠cct.lsu.edu 29 Oct '13

29 Oct '13

User: eschnett Date: 2013/10/29 07:53 PM Modified: /trunk/src/ GRHydro_BvecfromAvec.F90, GRHydro_ENOReconstruct_drv.F90, GRHydro_MP5Reconstruct_drv.F90, GRHydro_PPMMReconstruct_drv.F90, GRHydro_Reconstruct.F90, GRHydro_ReconstructPoly.F90, GRHydro_TVDReconstruct_drv.F90, GRHydro_WENOReconstruct_drv.F90, Utils.F90 Log: Declare loop variables OpenMP private File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_BvecfromAvec.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_BvecfromAvec.F90 2013-10-29 22:27:49 UTC (rev 585) +++ trunk/src/GRHydro_BvecfromAvec.F90 2013-10-30 00:53:23 UTC (rev 586) @@ -102,7 +102,7 @@ call CCTK_WARN(1,"Bvec from Avec start."); - !$OMP PARALLEL DO PRIVATE( sdet, isdet, local_spatial_order, & + !$OMP PARALLEL DO PRIVATE( i, j, k, sdet, isdet, local_spatial_order, & !$OMP Az_y, Ay_z, Ax_z, Az_x, Ay_x, Ax_y) do k = GRHydro_stencil, cctk_lsh(3)-GRHydro_stencil+1 do j = GRHydro_stencil, cctk_lsh(2)-GRHydro_stencil+1 File [modified]: GRHydro_ENOReconstruct_drv.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_ENOReconstruct_drv.F90 2013-10-29 22:27:49 UTC (rev 585) +++ trunk/src/GRHydro_ENOReconstruct_drv.F90 2013-10-30 00:53:23 UTC (rev 586) @@ -114,7 +114,7 @@ nz = cctk_lsh(3) !initialize trivial_rp to false - !$OMP PARALLEL DO + !$OMP PARALLEL DO PRIVATE(i,j,k) do k=1,nz do j=1,ny do i=1,nx File [modified]: GRHydro_MP5Reconstruct_drv.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_MP5Reconstruct_drv.F90 2013-10-29 22:27:49 UTC (rev 585) +++ trunk/src/GRHydro_MP5Reconstruct_drv.F90 2013-10-30 00:53:23 UTC (rev 586) @@ -129,7 +129,7 @@ nz = cctk_lsh(3) !initialize trivial_rp to false - !$OMP PARALLEL DO + !$OMP PARALLEL DO PRIVATE(i,j,k) do k=1,nz do j=1,ny do i=1,nx File [modified]: GRHydro_PPMMReconstruct_drv.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_PPMMReconstruct_drv.F90 2013-10-29 22:27:49 UTC (rev 585) +++ trunk/src/GRHydro_PPMMReconstruct_drv.F90 2013-10-30 00:53:23 UTC (rev 586) @@ -130,7 +130,7 @@ nz = cctk_lsh(3) !initialize trivial_rp to false - !$OMP PARALLEL DO + !$OMP PARALLEL DO PRIVATE(i,j,k) do k=1,nz do j=1,ny do i=1,nx File [modified]: GRHydro_Reconstruct.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_Reconstruct.F90 2013-10-29 22:27:49 UTC (rev 585) +++ trunk/src/GRHydro_Reconstruct.F90 2013-10-30 00:53:23 UTC (rev 586) @@ -54,7 +54,7 @@ end if ! Initialize plus and minus states - !$OMP PARALLEL DO + !$OMP PARALLEL DO PRIVATE(i,j,k) do k=1,cctk_lsh(3) do j=1,cctk_lsh(2) do i=1,cctk_lsh(1) File [modified]: GRHydro_ReconstructPoly.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_ReconstructPoly.F90 2013-10-29 22:27:49 UTC (rev 585) +++ trunk/src/GRHydro_ReconstructPoly.F90 2013-10-30 00:53:23 UTC (rev 586) @@ -127,7 +127,7 @@ !!$ Should change this. ! Initialize plus and minus states - !$OMP PARALLEL DO + !$OMP PARALLEL DO PRIVATE(i,j,k) do k=1,cctk_lsh(3) do j=1,cctk_lsh(2) do i=1,cctk_lsh(1) File [modified]: GRHydro_TVDReconstruct_drv.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_TVDReconstruct_drv.F90 2013-10-29 22:27:49 UTC (rev 585) +++ trunk/src/GRHydro_TVDReconstruct_drv.F90 2013-10-30 00:53:23 UTC (rev 586) @@ -143,7 +143,7 @@ nz = cctk_lsh(3) !initialize trivial_rp to false - !$OMP PARALLEL DO + !$OMP PARALLEL DO PRIVATE(i,j,k) do k=1,nz do j=1,ny do i=1,nx File [modified]: GRHydro_WENOReconstruct_drv.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/GRHydro_WENOReconstruct_drv.F90 2013-10-29 22:27:49 UTC (rev 585) +++ trunk/src/GRHydro_WENOReconstruct_drv.F90 2013-10-30 00:53:23 UTC (rev 586) @@ -188,7 +188,7 @@ nz = cctk_lsh(3) !initialize trivial_rp to false - !$OMP PARALLEL DO + !$OMP PARALLEL DO PRIVATE(i,j,k) do k=1,nz do j=1,ny do i=1,nx File [modified]: Utils.F90 Delta lines: +1 -1 =================================================================== --- trunk/src/Utils.F90 2013-10-29 22:27:49 UTC (rev 585) +++ trunk/src/Utils.F90 2013-10-30 00:53:23 UTC (rev 586) @@ -107,7 +107,7 @@ ny = cctk_lsh(2) nz = cctk_lsh(3) - !$OMP PARALLEL DO + !$OMP PARALLEL DO PRIVATE(i,j,k) do k=1,nz do j=1,ny do i=1,nx

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 585)
by rhaas＠tapir.caltech.edu 29 Oct '13

29 Oct '13

User: rhaas Date: 2013/10/29 05:27 PM Modified: /trunk/src/ GRHydro_Con2Prim.F90, GRHydro_Con2PrimAM.F90, GRHydro_Con2PrimHot.F90, GRHydro_Con2PrimM.F90, GRHydro_ENOReconstruct_drv.F90, GRHydro_EoSChangeGamma.F90, GRHydro_HLLE.F90, GRHydro_HLLEM.F90, GRHydro_HLLE_AM.F90, GRHydro_MP5Reconstruct_drv.F90, GRHydro_PPMMReconstruct_drv.F90, GRHydro_PPMReconstruct_drv.F90, GRHydro_PPMReconstruct_drv_opt.F90, GRHydro_Prim2Con.F90, GRHydro_Prim2ConM.F90, GRHydro_Reconstruct.F90, GRHydro_Source.F90, GRHydro_SourceAM.F90, GRHydro_SourceM.F90, GRHydro_TVDReconstruct_drv.F90, GRHydro_UpdateMask.F90, GRHydro_UpdateMaskM.F90, GRHydro_WENOReconstruct_drv.F90, Utils.F90 Log: GRHydro: use ashX instead of lsh(X) when declaring array sizes 1.) one should use ash when declaring array sizes 2.) intel11 chokes on lsh(X), this was already in the first Cray pointer patch File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Con2Prim.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Con2PrimAM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Con2PrimHot.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Con2PrimM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_ENOReconstruct_drv.F90 Delta lines: <Unavailable> File [modified]: GRHydro_EoSChangeGamma.F90 Delta lines: <Unavailable> File [modified]: GRHydro_HLLE.F90 Delta lines: <Unavailable> File [modified]: GRHydro_HLLEM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_HLLE_AM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_MP5Reconstruct_drv.F90 Delta lines: <Unavailable> File [modified]: GRHydro_PPMMReconstruct_drv.F90 Delta lines: <Unavailable> File [modified]: GRHydro_PPMReconstruct_drv.F90 Delta lines: <Unavailable> File [modified]: GRHydro_PPMReconstruct_drv_opt.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Prim2Con.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Prim2ConM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Reconstruct.F90 Delta lines: <Unavailable> File [modified]: GRHydro_Source.F90 Delta lines: <Unavailable> File [modified]: GRHydro_SourceAM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_SourceM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_TVDReconstruct_drv.F90 Delta lines: <Unavailable> File [modified]: GRHydro_UpdateMask.F90 Delta lines: <Unavailable> File [modified]: GRHydro_UpdateMaskM.F90 Delta lines: <Unavailable> File [modified]: GRHydro_WENOReconstruct_drv.F90 Delta lines: <Unavailable> File [modified]: Utils.F90 Delta lines: <Unavailable>

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 584)
by schnetter＠cct.lsu.edu 29 Oct '13

29 Oct '13

User: eschnett Date: 2013/10/29 03:39 PM Modified: /trunk/src/ GRHydro_Source.F90 Log: Add parentheses around macro definitions File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Source.F90 Delta lines: +9 -12 =================================================================== --- trunk/src/GRHydro_Source.F90 2013-10-29 15:19:01 UTC (rev 583) +++ trunk/src/GRHydro_Source.F90 2013-10-29 20:39:34 UTC (rev 584) @@ -9,23 +9,20 @@ ! Second order f.d. -#define DIFF_X_2(q) 0.5d0 * (q(i+1,j,k) - q(i-1,j,k)) * ida -#define DIFF_Y_2(q) 0.5d0 * (q(i,j+1,k) - q(i,j-1,k)) * idb -#define DIFF_Z_2(q) 0.5d0 * (q(i,j,k+1) - q(i,j,k-1)) * idc +#define DIFF_X_2(q) (0.5d0 * (q(i+1,j,k) - q(i-1,j,k)) * ida) +#define DIFF_Y_2(q) (0.5d0 * (q(i,j+1,k) - q(i,j-1,k)) * idb) +#define DIFF_Z_2(q) (0.5d0 * (q(i,j,k+1) - q(i,j,k-1)) * idc) ! Fourth order f.d. #if(1) -#define DIFF_X_4(q) ((q(i-2,j,k)-q(i+2,j,k)) + 8.d0 * (q(i+1,j,k) - q(i-1,j,k))) / 12.d0 * ida -#define DIFF_Y_4(q) ((q(i,j-2,k)-q(i,j+2,k)) + 8.d0 * (q(i,j+1,k) - q(i,j-1,k))) / 12.d0 * idb -#define DIFF_Z_4(q) ((q(i,j,k-2)-q(i,j,k+2)) + 8.d0 * (q(i,j,k+1) - q(i,j,k-1))) / 12.d0 * idc +#define DIFF_X_4(q) (((q(i-2,j,k)-q(i+2,j,k)) + 8.d0 * (q(i+1,j,k) - q(i-1,j,k))) / 12.d0 * ida) +#define DIFF_Y_4(q) (((q(i,j-2,k)-q(i,j+2,k)) + 8.d0 * (q(i,j+1,k) - q(i,j-1,k))) / 12.d0 * idb) +#define DIFF_Z_4(q) (((q(i,j,k-2)-q(i,j,k+2)) + 8.d0 * (q(i,j,k+1) - q(i,j,k-1))) / 12.d0 * idc) #else -#define DIFF_X_4(q) (-q(i+2,j,k) + 8.d0 * q(i+1,j,k) - 8.d0 * q(i-1,j,k) + \ - q(i-2,j,k)) / 12.d0 * ida -#define DIFF_Y_4(q) (-q(i,j+2,k) + 8.d0 * q(i,j+1,k) - 8.d0 * q(i,j-1,k) + \ - q(i,j-2,k)) / 12.d0 * idb -#define DIFF_Z_4(q) (-q(i,j,k+2) + 8.d0 * q(i,j,k+1) - 8.d0 * q(i,j,k-1) + \ - q(i,j,k-2)) / 12.d0 * idc +#define DIFF_X_4(q) ((-q(i+2,j,k) + 8.d0 * q(i+1,j,k) - 8.d0 * q(i-1,j,k) + q(i-2,j,k)) / 12.d0 * ida) +#define DIFF_Y_4(q) ((-q(i,j+2,k) + 8.d0 * q(i,j+1,k) - 8.d0 * q(i,j-1,k) + q(i,j-2,k)) / 12.d0 * idb) +#define DIFF_Z_4(q) ((-q(i,j,k+2) + 8.d0 * q(i,j,k+1) - 8.d0 * q(i,j,k-1) + q(i,j,k-2)) / 12.d0 * idc) #endif #include "cctk.h"

1 0

[svn:einsteintoolkit] GRHydro/trunk/src/ (Rev. 583)
by rhaas＠tapir.caltech.edu 29 Oct '13

29 Oct '13

User: rhaas Date: 2013/10/29 10:19 AM Modified: /trunk/src/ GRHydro_Wrappers.F90 Log: GRHydro: correct usage of sqrt(detg) in prim2con wrappers found by Ian Hawke File Changes: Directory: /trunk/src/ ====================== File [modified]: GRHydro_Wrappers.F90 Delta lines: +3 -3 =================================================================== --- trunk/src/GRHydro_Wrappers.F90 2013-10-28 16:36:15 UTC (rev 582) +++ trunk/src/GRHydro_Wrappers.F90 2013-10-29 15:19:01 UTC (rev 583) @@ -193,7 +193,7 @@ sdetg = sqrt(det) - call prim2con(handle, gxx, gxy, gxz, gyy, gyz, gzz, det, dens,& + call prim2con(handle, gxx, gxy, gxz, gyy, gyz, gzz, sdetg, dens,& sx, sy, sz, tau, rho, velx, vely, velz, epsilon,& press, w_lorentz) end subroutine @@ -219,7 +219,7 @@ sdetg = sqrt(det) - call prim2conpolytype(handle, gxx, gxy, gxz, gyy, gyz, gzz, det, dens,& + call prim2conpolytype(handle, gxx, gxy, gxz, gyy, gyz, gzz, sdetg, dens,& sx, sy, sz, tau, rho, velx, vely, velz, epsilon,& press, w_lorentz) end subroutine @@ -318,7 +318,7 @@ sdetg = sqrt(det) - call prim2conpolytypeM(handle, gxx, gxy, gxz, gyy, gyz, gzz, det,& + call prim2conpolytypeM(handle, gxx, gxy, gxz, gyy, gyz, gzz, sdetg,& dens, sx, sy, sz, tau, Bconsx, Bconsy, Bconsz,& rho, velx, vely, velz, epsilon, press, Bvecx,& Bvecy, Bvecz, w_lorentz)

1 0

[McLachlan]McLachlan branch, master, updated. ET_2013_05_v0-14-g18fe149
by Carpet git repositories 29 Oct '13

29 Oct '13

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "McLachlan". The branch, master has been updated via 18fe14998dac0391740f14a2bc5a8a2f014eb36b (commit) from 924d268009089a6a8b178997f56df248268dcfc6 (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 18fe14998dac0391740f14a2bc5a8a2f014eb36b Author: Frank Löffler <knarf(a)cct.lsu.edu> Date: Tue Oct 29 03:52:57 2013 -0500 par/CCZ4: remove parfiles because they don't work and never worked. See trac #1387. ----------------------------------------------------------------------- Summary of changes: par/CCZ4/bbh_ccz4.par | 493 ------------------------------------------------ par/CCZ4/gw1d_ccz4.par | 283 --------------------------- par/CCZ4/rbst_ccz4.par | 230 ---------------------- 3 files changed, 1006 deletions(-) delete mode 100644 par/CCZ4/bbh_ccz4.par delete mode 100644 par/CCZ4/gw1d_ccz4.par delete mode 100644 par/CCZ4/rbst_ccz4.par hooks/post-receive -- McLachlan

1 0

[svn:einsteintoolkit] NullConstr/trunk/src/ (Rev. 11)
by rhaas＠tapir.caltech.edu 28 Oct '13

28 Oct '13

User: rhaas Date: 2013/10/28 03:12 PM Modified: /trunk/src/ NullConstr_R00.F90, NullConstr_R01.F90, NullConstr_R0A.F90 Log: replace array operations by explicit loops File Changes: Directory: /trunk/src/ ====================== File [modified]: NullConstr_R00.F90 Delta lines: <Unavailable> File [modified]: NullConstr_R01.F90 Delta lines: <Unavailable> File [modified]: NullConstr_R0A.F90 Delta lines: <Unavailable>

1 0

[McLachlan]McLachlan branch, master, updated. ET_2013_05_v0-13-g924d268
by Carpet git repositories 28 Oct '13

28 Oct '13

This is an automated email from the git hooks/post-receive script. It was generated because a ref change was pushed to the repository containing the project "McLachlan". The branch, master has been updated via 924d268009089a6a8b178997f56df248268dcfc6 (commit) from caadfb1214b142be4512d86994fdfdac9a45206b (commit) Those revisions listed above that are new to this repository have not appeared on any other notification email; so we list those revisions in full, below. - Log ----------------------------------------------------------------- commit 924d268009089a6a8b178997f56df248268dcfc6 Author: Ian Hinder <ian.hinder(a)aei.mpg.de> Date: Mon Oct 28 17:51:56 2013 +0100 Regenerate code (Kranc commit 20e7a8feb560e494efe6b00db07af8b350fc9783) ----------------------------------------------------------------------- Summary of changes: ML_ADM/interface.ccl | 6 +++--- ML_ADMConstraints/interface.ccl | 4 ++-- ML_ADMQuantities/interface.ccl | 4 ++-- ML_BSSN/interface.ccl | 6 +++--- ML_CCZ4/interface.ccl | 6 +++--- ML_WaveToy/interface.ccl | 6 +++--- ML_WaveToy_CL/interface.ccl | 6 +++--- 7 files changed, 19 insertions(+), 19 deletions(-) hooks/post-receive -- McLachlan

1 0

Jump to page:

2026

2025

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

Commits