Fix an algebra error in the E_I term.

diff --git a/equations.frag b/equations.frag
index b8dbdde4216bb58e5381c7f0ed988b6af530e95b..190be60adc3a2dea2652371d55b72075a87af665 100644 (file)
--- a/equations.frag
+++ b/equations.frag
@@ -36,8 +36,8 @@ void main()
        float A11 = k1 * I_x * I_x;
        float A12 = k1 * I_x * I_y;
        float A22 = k1 * I_y * I_y;
-       float b1 = -k1 * I_t;
-       float b2 = -k1 * I_t;
+       float b1 = -k1 * I_t * I_x;
+       float b2 = -k1 * I_t * I_y;
 
        // Compute the second derivatives. First I_xx and I_xy.
        vec2 I_x_y_m2 = textureOffset(I_x_y_tex, tc, ivec2(-2,  0)).xy;

diff --git a/variational_refinement.txt b/variational_refinement.txt
index b94f7981636e59c3600a7b643d80cad3bff03dd3..29f99e6b480e6924a8cc17967cd9322005c6985c 100644 (file)
--- a/variational_refinement.txt
+++ b/variational_refinement.txt
@@ -153,8 +153,8 @@ iteration)
 
 we have the much more manageable
 
-  k1 I_x²    u + k1 I_x I_y v = - k1 I_t
-  k1 I_x I_y u + k1 I_y²    v = - k1 I_t
+  k1 I_x²    u + k1 I_x I_y v = - k1 I_t I_x
+  k1 I_x I_y u + k1 I_y²    v = - k1 I_t I_y
 
 ie., two linear equations in u and v. Now, you will notice two
 immediate problems by this equation set: