[cig-commits] commit: Remove dd_norm etc. because we should just compute it for each
Mercurial
hg at geodynamics.org
Fri Apr 13 16:41:19 PDT 2012
changeset: 153:dd661cbf8cc4
user: Walter Landry <wlandry at caltech.edu>
date: Fri Apr 13 16:34:40 2012 -0700
files: compute_coefficients/Interface.hxx
description:
Remove dd_norm etc. because we should just compute it for each
interface in compute_v_interface anyway.
diff -r 8945333f6754 -r dd661cbf8cc4 compute_coefficients/Interface.hxx
--- a/compute_coefficients/Interface.hxx Fri Apr 13 15:46:54 2012 -0700
+++ b/compute_coefficients/Interface.hxx Fri Apr 13 16:34:40 2012 -0700
@@ -10,12 +10,8 @@ class Interface
class Interface
{
public:
- FTensor::Tensor1<double,2> grid_pos,
- dd_pos[2][2], dd_norm[2][2], dd_tangent[2][2],
- sides_pos[2][2], sides_norm[2][2], sides_tangent[2][2],
- corner_pos[2][2], corner_norm[2][2], corner_tangent[2][2],
- anticorner_pos[2], anticorner_norm[2], anticorner_tangent[2],
- dp_pos[2], dp_norm[2], dp_tangent[2];
+ FTensor::Tensor1<double,2> grid_pos, dd_pos[2][2], sides_pos[2][2],
+ corner_pos[2][2], anticorner_pos[2], dp_pos[2];
bool intersect_dd[2][2], intersect_sides[2][2], intersect_corner[2][2],
intersect_anticorner[2][2], intersect_dp[2];
@@ -45,7 +41,7 @@ public:
static double intersection(const double &u_m, const double &u_0,
const double &u_p, const double &H)
{
- /* This assumes that there is only one intersection between u_p and u_m */
+ /* This returns the lowest intersection between u_m and u_p */
double du=(u_p-u_m)/2;
double ddu=u_p-2*u_0+u_m;
More information about the CIG-COMMITS
mailing list