/*
 *   surf - visualizing algebraic curves and algebraic surfaces
 *   Copyright (C) 1996-1997 Friedrich-Alexander-Universitaet
 *                           Erlangen-Nuernberg
 *                 1997-2000 Johannes Gutenberg-Universitaet Mainz
 *   Authors: Stephan Endrass, Hans Huelf, Ruediger Oertel,
 *            Kai Schneider, Ralf Schmitt, Johannes Beigel
 *
 *   This program is free software; you can redistribute it and/or modify
 *   it under the terms of the GNU General Public License as published by
 *   the Free Software Foundation; either version 2 of the License, or
 *   (at your option) any later version.
 *
 *   This program 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 General Public License for more details.
 *
 *   You should have received a copy of the GNU General Public License
 *   along with this program; if not, write to the Free Software
 *   Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
 *
 */


#include "stdio.h"
#include "NewClipAlgebraic.h"
#include "NewClipAlgebraicCube.h"

#include "gui_config.h"


void NewClipAlgebraicCube::init( void )
{
    // ----------------------------
    //  three polynomials:
    //  p[0] = r^2-(x-cx)^2
    //  p[1] = r^2-(y-cy)^2
    //  p[2] = r^2-(z-cz)^2
    //
    //  r  = clip_numeric.radius
    //  cx = clip_numeric.center_x
    //  cy = clip_numeric.center_y
    //  cz = clip_numeric.center_z
    // ----------------------------

    int i,j;

    for( i=0; i<3; i++ )
    {
        new_coeff_polyxyz( &p[i],4 );

        p[i].m[0].a = clip_numeric.radius*clip_numeric.radius;
        p[i].m[1].a = -1.0;
    }

    p[0].m[1].kx = p[1].m[1].ky = p[2].m[1].kz = 2;
    p[0].m[2].kx = p[1].m[2].ky = p[2].m[2].kz = 1;

    p[0].m[2].a = 2.0*clip_numeric.center_x;
    p[1].m[2].a = 2.0*clip_numeric.center_y;
    p[2].m[2].a = 2.0*clip_numeric.center_z;

    p[0].m[3].a = -clip_numeric.center_x*clip_numeric.center_x;
    p[1].m[3].a = -clip_numeric.center_y*clip_numeric.center_y;
    p[2].m[3].a = -clip_numeric.center_z*clip_numeric.center_z;

    for( i=0; i<3; i++ )
    {
        polyxyz_collect( &p[i] );
        polyxyz_set_degree( &p[i] );
        polyxyz_adjust( &p[i] );
	hp[i]  = new hornerpolyxyz ( p[i] );
    }

    for( i=3; i<MAIN_CLIP_AMOUNT_NUM; i++ )
    {
        p[i].n = 0;
    }
}