[Teem-commit] SF.net SVN: teem:[4582] teem/trunk/src/ten/tendEllipse.c

[<kin...@us...>]

Revision: 4582
          http://teem.svn.sourceforge.net/teem/?rev=4582&view=rev
Author:   kindlmann
Date:     2010-02-18 20:04:58 +0000 (Thu, 18 Feb 2010)

Log Message:
-----------
hacking tend ellipse

Modified Paths:
--------------
    teem/trunk/src/ten/tendEllipse.c

Modified: teem/trunk/src/ten/tendEllipse.c
===================================================================
--- teem/trunk/src/ten/tendEllipse.c	2010-02-18 20:04:24 UTC (rev 4581)
+++ teem/trunk/src/ten/tendEllipse.c	2010-02-18 20:04:58 UTC (rev 4582)
@@ -117,16 +117,27 @@
     }
     conf = tdata[0];
     if (conf > cthresh) {
+      double eval0, eval1, dd;
       Dxx = tdata[1];
       Dxy = tdata[2];
       Dyy = tdata[3];
+      dd = Dxx - Dyy;
+      eval0 = 0.5*(-Dxx + sqrt(4*Dxy*Dxy + dd*dd) - Dyy);
+      eval1 = 0.5*(-Dxx - sqrt(4*Dxy*Dxy + dd*dd) - Dyy);
       fprintf(file, "gsave\n");
       fprintf(file, "matrix currentmatrix\n");
       fprintf(file, "[%g %g %g %g %g %g] concat\n",
               Dxx, -Dxy, -Dxy, Dyy, px, py);
       fprintf(file, "0 0 %g 0 360 arc closepath\n", gscale);
       fprintf(file, "setmatrix\n");
-      fprintf(file, "fill\n");
+      if (eval0 * eval1 < 0) {
+        fprintf(file, "gsave\n");
+        fprintf(file, "0.15 setgray\n");
+        fprintf(file, "fill\n");
+        fprintf(file, "grestore\n");
+      } else {
+        fprintf(file, "fill\n");
+      }
       fprintf(file, "grestore\n");
     }
     tdata += 4;


This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.