//--------------------------------------------------------------------------
// Copyright (C) 2001 Geoffrey Furnish
// Copyright (C) 2001-2014 Alan W. Irwin
// Copyright (C) 2002 Maurice LeBrun
// Copyright (C) 2002 Joao Cardoso
// Copyright (C) 2004 Andrew Ross
//
// This file is part of PLplot.
//
// PLplot is free software; you can redistribute it and/or modify
// it under the terms of the GNU Library General Public License as published by
// the Free Software Foundation; version 2 of the License.
//
// PLplot is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU Library General Public License for more details.
//
// You should have received a copy of the GNU Library General Public License
// along with PLplot; if not, write to the Free Software
// Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
//--------------------------------------------------------------------------
//--------------------------------------------------------------------------
// Implementation of PLplot example 8 in Java.
//--------------------------------------------------------------------------
package plplot.examples;
import plplot.core.*;
import static plplot.core.plplotjavacConstants.*;
import java.lang.Math;
class x08 {
// These values must be odd, for the middle
// of the index range to be an integer, and thus
// to correspond to the exact floating point centre
// of the sombrero.
static final int XPTS = 35;
static final int YPTS = 45;
PLStream pls = new PLStream();
static double alt[] = { 60.0, 40.0 };
static double az[] = { 30.0, -30.0 };
static String[] title =
{
"#frPLplot Example 8 - Alt=60, Az=30",
"#frPLplot Example 8 - Alt=40, Az=-30",
};
// cmap1_init1
// Initializes color map 1 in HLS space.
// Basic grayscale variation from half-dark (which makes more interesting
// looking plot compared to dark) to light.
// An interesting variation on this:
// s[1] = 1.0
void cmap1_init( int gray )
{
double [] i = new double[2];
double [] h = new double[2];
double [] l = new double[2];
double [] s = new double[2];
boolean [] rev = new boolean[2];
i[0] = 0.0; // left boundary
i[1] = 1.0; // right boundary
if ( gray == 1 )
{
h[0] = 0.0; // hue -- low: red (arbitrary if s=0)
h[1] = 0.0; // hue -- high: red (arbitrary if s=0)
l[0] = 0.5; // lightness -- low: half-dark
l[1] = 1.0; // lightness -- high: light
s[0] = 0.0; // minimum saturation
s[1] = 0.0; // minimum saturation
}
else
{
h[0] = 240; // blue -> green -> yellow ->
h[1] = 0; // -> red
l[0] = 0.6;
l[1] = 0.6;
s[0] = 0.8;
s[1] = 0.8;
}
rev[0] = false; // interpolate on front side of colour wheel.
rev[1] = false; // interpolate on front side of colour wheel.
pls.scmap1n( 256 );
pls.scmap1l( false, i, h, l, s, rev );
}
// Does a series of 3-d plots for a given data set, with different viewing
// options in each plot.
public static void main( String[] args )
{
new x08( args );
}
public x08( String[] args )
{
int i, j, k;
final int LEVELS = 10;
double[] x = new double[ XPTS ];
double[] y = new double[ YPTS ];
double[][] z = new double[XPTS][YPTS];
double clevel[] = new double[LEVELS];
double clev_null [] = new double[0];
double xx, yy, r;
double zmin = Double.MAX_VALUE, zmax = Double.MIN_VALUE;
int ifshade;
int indexxmin = 0;
int indexxmax = XPTS;
int[] indexymin = new int[XPTS];
int[] indexymax = new int[XPTS];
double[][] zlimited = new double[XPTS][YPTS];
double dx = 2. / ( XPTS - 1 );
double dy = 2. / ( YPTS - 1 );
// parameters of ellipse (in x, y index coordinates) that limits the data.
// x0, y0 correspond to the exact floating point centre of the index
// range.
double x0 = 0.5 * ( XPTS - 1 );
double a = 0.9 * x0;
double y0 = 0.5 * ( YPTS - 1 );
double b = 0.7 * y0;
double square_root;
// Parse and process command line arguments.
pls.parseopts( args, PL_PARSE_FULL | PL_PARSE_NOPROGRAM );
boolean rosen = false;
// Initialize plplot.
pls.init();
for ( i = 0; i < XPTS; i++ )
{
x[i] = -1. + dx * i;
if ( rosen )
x[i] *= 1.5;
}
for ( j = 0; j < YPTS; j++ )
{
y[j] = -1. + dy * j;
if ( rosen )
y[j] += 0.5;
}
for ( i = 0; i < XPTS; i++ )
{
xx = x[i];
for ( j = 0; j < YPTS; j++ )
{
yy = y[j];
if ( rosen )
{
z[i][j] = Math.pow( 1. - xx, 2. ) + 100 *
Math.pow( yy - Math.pow( xx, 2. ), 2. );
// The log argument might be zero for just the right grid.
if ( z[i][j] > 0. )
z[i][j] = Math.log( z[i][j] );
else
z[i][j] = -5.; // -MAXFLOAT would mess-up up the scale
}
else
{
r = Math.sqrt( xx * xx + yy * yy );
z[i][j] = Math.exp( -r * r ) * Math.cos( 2.0 * Math.PI * r );
}
if ( zmin > z[i][j] )
zmin = z[i][j];
if ( zmax < z[i][j] )
zmax = z[i][j];
}
}
for ( i = indexxmin; i < indexxmax; i++ )
{
square_root = Math.sqrt( 1. - Math.min( 1., Math.pow( ( i - x0 ) / a, 2. ) ) );
// Add 0.5 to find nearest integer and therefore preserve symmetry
// with regard to lower and upper bound of y range.
indexymin[i] = Math.max( 0, (int) ( 0.5 + y0 - b * square_root ) );
// indexymax calculated with the convention that it is 1
// greater than highest valid index.
indexymax[i] = Math.min( YPTS, 1 + (int) ( 0.5 + y0 + b * square_root ) );
for ( j = indexymin[i]; j < indexymax[i]; j++ )
zlimited[i][j] = z[i][j];
}
double step = ( zmax - zmin ) / ( LEVELS + 1 );
for ( i = 0; i < LEVELS; i++ )
clevel[i] = zmin + step * ( i + 1 );
pls.lightsource( 1., 1., 1. );
for ( k = 0; k < 2; k++ )
{
for ( ifshade = 0; ifshade < 5; ifshade++ )
{
pls.adv( 0 );
pls.vpor( 0.0, 1.0, 0.0, 0.9 );
pls.wind( -1.0, 1.0, -0.9, 1.1 );
pls.col0( 3 );
pls.mtex( "t", 1.0, 0.5, 0.5, title[k] );
pls.col0( 1 );
if ( rosen )
pls.w3d( 1.0, 1.0, 1.0, -1.5, 1.5, -0.5, 1.5, zmin, zmax,
alt[k], az[k] );
else
pls.w3d( 1.0, 1.0, 1.0, -1.0, 1.0, -1.0, 1.0, zmin, zmax,
alt[k], az[k] );
pls.box3( "bnstu", "x axis", 0.0, 0,
"bnstu", "y axis", 0.0, 0,
"bcdmnstuv", "z axis", 0.0, 0 );
pls.col0( 2 );
switch ( ifshade )
{
case 0: // diffuse light surface plot
cmap1_init( 1 );
pls.surf3d( x, y, z, 0, clev_null );
break;
case 1: // magnitude colored plot
cmap1_init( 0 );
pls.surf3d( x, y, z, MAG_COLOR, clev_null );
break;
case 2: // magnitude colored plot with faceted squares
cmap1_init( 0 );
pls.surf3d( x, y, z, MAG_COLOR | FACETED, clev_null );
break;
case 3: // magnitude colored plot with contours
cmap1_init( 0 );
pls.surf3d( x, y, z, MAG_COLOR | SURF_CONT | BASE_CONT, clevel );
break;
case 4: // magnitude colored plot with contours and index limits.
cmap1_init( 0 );
pls.surf3dl( x, y, zlimited, MAG_COLOR | SURF_CONT | BASE_CONT, clevel, indexxmin, indexymin, indexymax );
}
}
}
pls.end();
}
}
//--------------------------------------------------------------------------
// End of x08.java
//--------------------------------------------------------------------------