[cig-commits] commit: Fix a bug where we were using the wrong dp_jump in an intermediate calculation.

[Mercurial]

changeset:   128:22520c38035a
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Mar 28 05:50:40 2012 -0700
files:       compute_jumps.cxx
description:
Fix a bug where we were using the wrong dp_jump in an intermediate calculation.


diff -r cd3b7d3189e4 -r 22520c38035a compute_jumps.cxx
--- a/compute_jumps.cxx	Wed Mar 28 05:15:13 2012 -0700
+++ b/compute_jumps.cxx	Wed Mar 28 05:50:40 2012 -0700
@@ -42,7 +42,7 @@ void compute_jumps(const double &eta_jum
 
   FTensor::Tensor3_christof<double,2,2> ddz_jump_nt;
   ddz_jump_nt(a,t,t)=eta_jump*ddv(a);
-  ddz_jump_nt(a,n,n)=-ddz_jump_nt(a,t,t) + dp_jump(a);
+  ddz_jump_nt(a,n,n)=-ddz_jump_nt(a,t,t) + dp_jump_nt(a);
   ddz_jump_nt(a,n,t)=-eta_jump*ddv(b)*ddel(a,b);
 
   /* It would be nice if there were a way to do this directly instead