diff options
| author | lanfer <lanfer@10716dce-81a3-4424-a2c8-48026a0d3035> | 2000-05-12 14:02:41 +0000 |
|---|---|---|
| committer | lanfer <lanfer@10716dce-81a3-4424-a2c8-48026a0d3035> | 2000-05-12 14:02:41 +0000 |
| commit | 27dec4c0b1c40a77c639071bcedfde881ea4e14a (patch) | |
| tree | 7880b2367a7a37698dc90f4e287f7e196544fe71 /src | |
| parent | 3f803d9c288cafde42592e6166ee289c98c7daff (diff) | |
adding downsampling
git-svn-id: http://svn.cactuscode.org/arrangements/CactusPUGH/PUGHSlab/trunk@15 10716dce-81a3-4424-a2c8-48026a0d3035
Diffstat (limited to 'src')
| -rw-r--r-- | src/GetLocalLine.c | 12 |
1 files changed, 6 insertions, 6 deletions
diff --git a/src/GetLocalLine.c b/src/GetLocalLine.c index 9130fb3..f533a45 100644 --- a/src/GetLocalLine.c +++ b/src/GetLocalLine.c @@ -76,8 +76,8 @@ int GetLocalLine(cGH *GH, printf("\n+++++++++\n"); printf("S_l [XYZ] : %d %d %d \n", S_l[0], S_l[1], S_l[2]); - printf("u_l [XYZ] : %d %d %d \n", - u_l[0], u_l[1], u_l[2]); + printf("u_l [XYZ] : %d %d %d ds %d \n", + u_l[0], u_l[1], u_l[2], ds); for (idim=0;idim<vardim;idim++) { @@ -140,7 +140,7 @@ int GetLocalLine(cGH *GH, for (idim=0;idim<vardim;idim++) { tmp1+= (S_l[idim]-P_s[idim][iface])*n_s[idim][iface]; - tmp2+= (u_l[idim]*n_s[idim][iface]); + tmp2+= (ds*u_l[idim]*n_s[idim][iface]); } if (tmp2==0) { @@ -189,7 +189,7 @@ int GetLocalLine(cGH *GH, if (notcoll[iface]==1 && lmd_s[iface]!=BAD) { for (idim=0;idim<vardim;idim++) - t_pnt1[idim]=S_l[idim]+lmd_s[iface]*u_l[idim]; + t_pnt1[idim]=S_l[idim]+lmd_s[iface]*ds*u_l[idim]; printf("testing for local intersec. (face %d lamda %+4.2f): %4.2f %4.2f %4.2f\n", iface, lmd_s[iface], t_pnt1[0],t_pnt1[1],t_pnt1[2]); @@ -289,8 +289,8 @@ int GetLocalLine(cGH *GH, nlocpoints[0]= ABS(lmd[0]-lmd[1])+1; for (idim=0;idim<vardim;idim++) { - locstart[idim] = (S_l[idim]+lmd[0]*u_l[idim]); - locend[idim] = (S_l[idim]+lmd[1]*u_l[idim]); + locstart[idim] = (S_l[idim]+lmd[0]*ds*u_l[idim]); + locend[idim] = (S_l[idim]+lmd[1]*ds*u_l[idim]); } } else |