[Plplot-cvs] SF.net SVN: plplot:[12766] trunk/examples/java/x22.java

[<and...@us...>]

Revision: 12766
          http://sourceforge.net/p/plplot/code/12766
Author:   andrewross
Date:     2013-11-27 12:16:50 +0000 (Wed, 27 Nov 2013)
Log Message:
-----------
Update java version of example 22.

Modified Paths:
--------------
    trunk/examples/java/x22.java

Modified: trunk/examples/java/x22.java
===================================================================
--- trunk/examples/java/x22.java	2013-11-27 12:12:14 UTC (rev 12765)
+++ trunk/examples/java/x22.java	2013-11-27 12:16:50 UTC (rev 12766)
@@ -33,6 +33,18 @@
 
 import java.lang.Math;
 
+class Transform implements PLCallbackCT {
+    public void coordTransform( double x, double y, double[] xt, double[] yt, Object data )
+    {
+        double xmax;
+
+        xmax = 10.0;
+
+        xt[0] = x;
+        yt[0] = y / 4.0 * ( 3 - Math.cos( Math.PI * x / xmax ) );
+    }
+}
+
 class x22 {
     double MIN( double x, double y ) { return ( x < y ? x : y ); }
     double MAX( double x, double y ) { return ( x > y ? x : y ); }
@@ -43,7 +55,7 @@
     double[][] v  = null;
     double[][] xg = null;
     double[][] yg = null;
-    int nx, ny;
+    int nx, ny, nc, nseg;
 
     // Vector plot of the circulation about the origin
     void circulation()
@@ -84,7 +96,7 @@
     }
 
     // Vector plot of flow through a constricted pipe
-    void constriction()
+    void constriction( int astyle )
     {
         int    i, j;
         double dx, dy, x, y;
@@ -112,7 +124,7 @@
                 if ( Math.abs( y ) < b )
                 {
                     dbdx = ymax / 4.0 * Math.sin( Math.PI * x / xmax ) *
-                           y / b;
+                           Math.PI / xmax * y / b;
                     u[i][j] = Q * ymax / b;
                     v[i][j] = dbdx * u[i][j];
                 }
@@ -125,12 +137,70 @@
         }
 
         pls.env( xmin, xmax, ymin, ymax, 0, 0 );
-        pls.lab( "(x)", "(y)", "#frPLplot Example 22 - constriction" );
+        pls.lab( "(x)", "(y)", "#frPLplot Example 22 - constriction (arrow style " + astyle + ")" );
         pls.col0( 2 );
-        pls.vect( u, v, -0.5, xg, yg );
+        pls.vect( u, v, -1.0, xg, yg );
         pls.col0( 1 );
     }
 
+//
+// Vector plot of flow through a constricted pipe
+// with a coordinate transform
+//
+    void constriction2( )
+    {
+        int   i, j;
+        double dx, dy, x, y;
+        double xmin, xmax, ymin, ymax;
+        double Q, b, dbdx;
+        double[] clev = new double[nc];
+        Transform transform = new Transform();
+
+        dx = 1.0;
+        dy = 1.0;
+
+        xmin = -nx / 2 * dx;
+        xmax = nx / 2 * dx;
+        ymin = -ny / 2 * dy;
+        ymax = ny / 2 * dy;
+
+        pls.stransform( transform, xmax );
+
+        Q = 2.0;
+        for ( i = 0; i < nx; i++ )
+        {
+            x = ( i - nx / 2 + 0.5 ) * dx;
+            for ( j = 0; j < ny; j++ )
+            {
+                y = ( j - ny / 2 + 0.5 ) * dy;
+                xg[i][j] = x;
+                yg[i][j] = y;
+                b       = ymax / 4.0 * ( 3 - Math.cos( Math.PI * x / xmax ) );
+                u[i][j] = Q * ymax / b;
+                v[i][j] = 0.0;
+            }
+        }
+
+        for ( i = 0; i < nc; i++ )
+        {
+            clev[i] = Q + i * Q / ( nc - 1 );
+        }
+
+        pls.env( xmin, xmax, ymin, ymax, 0, 0 );
+        pls.lab( "(x)", "(y)", "#frPLplot Example 22 - constriction with plstransform" );
+        pls.col0( 2 );
+        pls.shades( u, xmin + dx / 2, xmax - dx / 2, 
+            ymin + dy / 2, ymax - dy / 2,
+            clev, 0, 1, 1.0, 0, xg, yg );
+        pls.vect( u, v, -1.0, xg, yg );
+        // Plot edges using plpath (which accounts for coordinate transformation) rather than plline
+        pls.path( nseg, xmin, ymax, xmax, ymax );
+        pls.path( nseg, xmin, ymin, xmax, ymin );
+        pls.col0( 1 );
+
+        pls.stransform( null, null );
+}
+
     // Vector plot of the gradient of a shielded potential (see example 9)
     void potential()
     {
@@ -270,6 +340,8 @@
 
         nx = 20;
         ny = 20;
+        nc = 11;
+        nseg = 20;
 
         // Allocate arrays
         u  = new double[nx][ny];
@@ -284,14 +356,18 @@
         // Set arrow style using arrow_x and arrow_y then
         // plot uMath.sing these arrows.
         pls.svect( arrow_x, arrow_y, fill );
-        constriction();
+        constriction( 1 );
 
         // Set arrow style using arrow2_x and arrow2_y then
         // plot using these filled arrows.
         fill = true;
         pls.svect( arrow2_x, arrow2_y, fill );
-        constriction();
+        constriction( 2 );
 
+        constriction2();
+
+        pls.svect(null,null,false);
+
         potential();
 
         pls.end();

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