//--------------------------------------------------------------------------
// Copyright (C) 2001  Geoffrey Furnish
// Copyright (C) 2002-2014 Alan W. Irwin
// 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 18 in Java.
//--------------------------------------------------------------------------

package plplot.examples;

import plplot.core.*;
import static plplot.core.plplotjavacConstants.*;

import java.lang.Math;
import java.text.*;
class x18 {
    static final int NPTS  = 1000;
    static int       opt[] = { 1, 0, 1, 0 };
    static double    alt[] = { 20.0, 35.0, 50.0, 65.0 };
    static double    az[]  = { 30.0, 40.0, 50.0, 60.0 };
    PLStream         pls   = new PLStream();

    double THETA( int a )
    {
        return 2. * Math.PI * a / 20.;
    }
    double PHI( int a )
    {
        return Math.PI * a / 20.1;
    }
    void test_poly( int k )
    {
        int i, j;
        boolean [][] draw = {
            { true,  true,  true,  true  },
            { true,  false, true,  false },
            { false, true,  false, true  },
            { true,  true,  false, false }
        };

        double [] x = new double [5];
        double [] y = new double [5];
        double [] z = new double [5];

        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( 1 );
        pls.w3d( 1.0, 1.0, 1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, 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 );

        // x = r sin(phi) cos(theta)
        // y = r sin(phi) sin(theta)
        // z = r cos(phi)
        // r = 1 :=)

        for ( i = 0; i < 20; i++ )
        {
            for ( j = 0; j < 20; j++ )
            {
                x[0] = Math.sin( PHI( j ) ) * Math.cos( THETA( i ) );
                y[0] = Math.sin( PHI( j ) ) * Math.sin( THETA( i ) );
                z[0] = Math.cos( PHI( j ) );

                x[1] = Math.sin( PHI( j + 1 ) ) * Math.cos( THETA( i ) );
                y[1] = Math.sin( PHI( j + 1 ) ) * Math.sin( THETA( i ) );
                z[1] = Math.cos( PHI( j + 1 ) );

                x[2] = Math.sin( PHI( j + 1 ) ) * Math.cos( THETA( i + 1 ) );
                y[2] = Math.sin( PHI( j + 1 ) ) * Math.sin( THETA( i + 1 ) );
                z[2] = Math.cos( PHI( j + 1 ) );

                x[3] = Math.sin( PHI( j ) ) * Math.cos( THETA( i + 1 ) );
                y[3] = Math.sin( PHI( j ) ) * Math.sin( THETA( i + 1 ) );
                z[3] = Math.cos( PHI( j ) );

                x[4] = Math.sin( PHI( j ) ) * Math.cos( THETA( i ) );
                y[4] = Math.sin( PHI( j ) ) * Math.sin( THETA( i ) );
                z[4] = Math.cos( PHI( j ) );

                pls.poly3( x, y, z, draw[k], true );
            }
        }

        pls.col0( 3 );
        pls.mtex( "t", 1.0, 0.5, 0.5, "unit radius sphere" );
    }
    // 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 x18( args );
    }

    public x18( String[] args )
    {
        NumberFormat nf = NumberFormat.getNumberInstance();

        int          i, k;
        double       r;

        // Parse and process command line arguments.
        pls.parseopts( args, PL_PARSE_FULL | PL_PARSE_NOPROGRAM );

        // Initialize plplot.
        pls.init();

        for ( k = 0; k < 4; k++ )
            test_poly( k );

        double[] x = new double[NPTS];
        double[] y = new double[NPTS];
        double[] z = new double[NPTS];

        // From the mind of a sick and twisted physicist...

        for ( i = 0; i < NPTS; i++ )
        {
            z[i] = -1. + 2. * i / NPTS;

            // Pick one ...

            // r    = 1. - ( (double) i / (double) NPTS );
            r = z[i];

            x[i] = r * Math.cos( 2. * Math.PI * 6. * i / NPTS );
            y[i] = r * Math.sin( 2. * Math.PI * 6. * i / NPTS );
        }

        for ( k = 0; k < 4; k++ )
        {
            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( 1 );
            pls.w3d( 1.0, 1.0, 1.0, -1.0, 1.0, -1.0, 1.0, -1.0, 1.0, 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 ( opt[k] > 0 )
                pls.line3( x, y, z );
            else
                // U+22C5 DOT OPERATOR.
                pls.string3( x, y, z, "⋅" );

            pls.col0( 3 );
            String title = "#frPLplot Example 18 - Alt=" +
                           nf.format( (int) alt[k] ) +
                           ", Az=" +
                           nf.format( (int) az[k] );
            pls.mtex( "t", 1.0, 0.5, 0.5, title );
        }

        pls.end();
    }
}

//--------------------------------------------------------------------------
//                              End of x18.java
//--------------------------------------------------------------------------