[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