/*
 *   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 "NewClipAlgebraicOcta.h"

#include "gui_config.h"


void NewClipAlgebraicOcta::init( void )
{
    // ----------------------------
    //  four polynomials:
    //  p[0] = (r+(x-cx)+(y-cy)+(z-cz))*(r-(x-cx)-(y-cy)-(z-cz))
    //  p[1] = (r-(x-cx)+(y-cy)+(z-cz))*(r+(x-cx)-(y-cy)-(z-cz))
    //  p[2] = (r+(x-cx)-(y-cy)+(z-cz))*(r-(x-cx)+(y-cy)-(z-cz))
    //  p[3] = (r+(x-cx)+(y-cy)-(z-cz))*(r-(x-cx)-(y-cy)+(z-cz))
    //
    //  r  = clip_numeric.radius
    //  cx = clip_numeric.center_x
    //  cy = clip_numeric.center_y
    //  cz = clip_numeric.center_z
    // ----------------------------

    int i;
    polyxyz h[8];

    for( i=0; i<8; i++ )
    {
        new_coeff_polyxyz( &h[i],4 );
        h[i].m[1].kx = h[i].m[2].ky = h[i].m[3].kz = 1;
    }

    h[0].m[0].a = clip_numeric.radius - clip_numeric.center_x - clip_numeric.center_y - clip_numeric.center_z;
    h[1].m[0].a = clip_numeric.radius + clip_numeric.center_x - clip_numeric.center_y - clip_numeric.center_z;
    h[2].m[0].a = clip_numeric.radius - clip_numeric.center_x + clip_numeric.center_y - clip_numeric.center_z;
    h[3].m[0].a = clip_numeric.radius - clip_numeric.center_x - clip_numeric.center_y + clip_numeric.center_z;

    h[0].m[1].a = 1.0; h[0].m[2].a = 1.0; h[0].m[3].a = 1.0;
    h[1].m[1].a =-1.0; h[1].m[2].a = 1.0; h[1].m[3].a = 1.0;
    h[2].m[1].a = 1.0; h[2].m[2].a =-1.0; h[2].m[3].a = 1.0;
    h[3].m[1].a = 1.0; h[3].m[2].a = 1.0; h[3].m[3].a =-1.0;

    h[4].m[0].a = clip_numeric.radius + clip_numeric.center_x + clip_numeric.center_y + clip_numeric.center_z;
    h[5].m[0].a = clip_numeric.radius - clip_numeric.center_x + clip_numeric.center_y + clip_numeric.center_z;
    h[6].m[0].a = clip_numeric.radius + clip_numeric.center_x - clip_numeric.center_y + clip_numeric.center_z;
    h[7].m[0].a = clip_numeric.radius + clip_numeric.center_x + clip_numeric.center_y - clip_numeric.center_z;

    h[4].m[1].a =-1.0; h[4].m[2].a =-1.0; h[4].m[3].a =-1.0;
    h[5].m[1].a = 1.0; h[5].m[2].a =-1.0; h[5].m[3].a =-1.0;
    h[6].m[1].a =-1.0; h[6].m[2].a = 1.0; h[6].m[3].a =-1.0;
    h[7].m[1].a =-1.0; h[7].m[2].a =-1.0; h[7].m[3].a = 1.0;

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

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

    for( i=0; i<8; i++ )
    {
        delete_coeff_polyxyz( &h[i] );
    }
}