//-------------------------------------------------------------------------- // $Id: x08.java 12240 2012-10-04 22:10:12Z andrewross $ //-------------------------------------------------------------------------- //-------------------------------------------------------------------------- // Copyright (C) 2001 Geoffrey Furnish // Copyright (C) 2001, 2002, 2003, 2004 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 { static final int XPTS = 35; static final int YPTS = 46; PLStream pls = new PLStream(); static double alt[] = { 60.0, 20.0 }; static double az[] = { 30.0, 60.0 }; static String[] title = { "#frPLplot Example 8 - Alt=60, Az=30", "#frPLplot Example 8 - Alt=20, Az=60", }; // 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; // Parse and process command line arguments. pls.parseopts( args, PL_PARSE_FULL | PL_PARSE_NOPROGRAM ); boolean rosen = true; // Initialize plplot. pls.init(); for ( i = 0; i < XPTS; i++ ) { x[i] = (double) ( i - ( XPTS / 2 ) ) / (double) ( XPTS / 2 ); if ( rosen ) x[i] *= 1.5; } for ( j = 0; j < YPTS; j++ ) { y[j] = (double) ( j - ( YPTS / 2 ) ) / (double) ( YPTS / 2 ); 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 may 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]; } } 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 < 4; 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 ); if ( ifshade == 0 ) // diffuse light surface plot { cmap1_init( 1 ); // with new interface haven't yet made null work so have // to put in specific zero-length array. pls.surf3d( x, y, z, 0, clev_null ); } else if ( ifshade == 1 ) // magnitude colored plot { cmap1_init( 0 ); pls.surf3d( x, y, z, MAG_COLOR, clev_null ); } else if ( ifshade == 2 ) // magnitude colored plot with faceted squares { cmap1_init( 0 ); pls.surf3d( x, y, z, MAG_COLOR | FACETED, clev_null ); } else // magnitude colored plot with contours { cmap1_init( 0 ); pls.surf3d( x, y, z, MAG_COLOR | SURF_CONT | BASE_CONT, clevel ); } } } pls.end(); } } //-------------------------------------------------------------------------- // End of x08.java //--------------------------------------------------------------------------