/*
* 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 <math.h>
#include <stdio.h>
#include "NewClipAlgebraic.h"
#include "NewClipAlgebraicIco.h"
#include "gui_config.h"
void NewClipAlgebraicIco::init( void )
{
// ------------------------------------------------------------------------
// ten polynomials:
//
// p[0] = (2*r+0*(x-cx)+a*(y-cy)+b*(z-cz))*(r+0*(x-cx)-a*(y-cy)-b*(z-cz))
// p[1] = (2*r+b*(x-cx)+0*(y-cy)+a*(z-cz))*(r-b*(x-cx)+0*(y-cy)-a*(z-cz))
// p[2] = (2*r+a*(x-cx)+b*(y-cy)+0*(z-cz))*(r-a*(x-cx)-b*(y-cy)+0*(z-cz))
//
// p[3] = (2*r+0*(x-cx)+a*(y-cy)-b*(z-cz))*(r+0*(x-cx)-a*(y-cy)+b*(z-cz))
// p[4] = (2*r-b*(x-cx)+0*(y-cy)+a*(z-cz))*(r+b*(x-cx)+0*(y-cy)-a*(z-cz))
// p[5] = (2*r+a*(x-cx)-b*(y-cy)+0*(z-cz))*(r-a*(x-cx)+b*(y-cy)+0*(z-cz))
//
// p[6] = (2*r+c*(x-cx)+c*(y-cy)+c*(z-cz))*(r-c*(x-cx)-c*(y-cy)-c*(z-cz))
// p[7] = (2*r+c*(x-cx)-c*(y-cy)-c*(z-cz))*(r-c*(x-cx)+c*(y-cy)+c*(z-cz))
// p[8] = (2*r-c*(x-cx)+c*(y-cy)-c*(z-cz))*(r+c*(x-cx)-c*(y-cy)+c*(z-cz))
// p[9] = (2*r-c*(x-cx)-c*(y-cy)+c*(z-cz))*(r+c*(x-cx)+c*(y-cy)-c*(z-cz))
//
// r = clip_numeric.radius
// cx = clip_numeric.center_x
// cy = clip_numeric.center_y
// cz = clip_numeric.center_z
//
// a = sqrt(10-2*sqrt(5))
// b = sqrt(50-22*sqrt(5))
// c = 2*sqrt(5-2*sqrt(5))
// ------------------------------------------------------------------------
double a = sqrt(10.0-2.0*sqrt(5.0));
double b = sqrt(50.0-22.0*sqrt(5.0));
double c = 2.0*sqrt(5.0-2.0*sqrt(5.0));
polyxyz h[20];
int i;
for( i=0; i<20; 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 = 2*clip_numeric.radius - c*clip_numeric.center_x - c*clip_numeric.center_y - c*clip_numeric.center_z;
h[1].m[0].a = 2*clip_numeric.radius - c*clip_numeric.center_x + c*clip_numeric.center_y + c*clip_numeric.center_z;
h[2].m[0].a = 2*clip_numeric.radius + c*clip_numeric.center_x - c*clip_numeric.center_y + c*clip_numeric.center_z;
h[3].m[0].a = 2*clip_numeric.radius + c*clip_numeric.center_x + c*clip_numeric.center_y - c*clip_numeric.center_z;
h[0].m[1].a = c; h[0].m[2].a = c; h[0].m[3].a = c;
h[1].m[1].a = c; h[1].m[2].a = -c; h[1].m[3].a = -c;
h[2].m[1].a = -c; h[2].m[2].a = c; h[2].m[3].a = -c;
h[3].m[1].a = -c; h[3].m[2].a = -c; h[3].m[3].a = c;
h[10].m[0].a = 2*clip_numeric.radius + c*clip_numeric.center_x + c*clip_numeric.center_y + c*clip_numeric.center_z;
h[11].m[0].a = 2*clip_numeric.radius + c*clip_numeric.center_x - c*clip_numeric.center_y - c*clip_numeric.center_z;
h[12].m[0].a = 2*clip_numeric.radius - c*clip_numeric.center_x + c*clip_numeric.center_y - c*clip_numeric.center_z;
h[13].m[0].a = 2*clip_numeric.radius - c*clip_numeric.center_x - c*clip_numeric.center_y + c*clip_numeric.center_z;
h[10].m[1].a = -c; h[10].m[2].a = -c; h[10].m[3].a = -c;
h[11].m[1].a = -c; h[11].m[2].a = c; h[11].m[3].a = c;
h[12].m[1].a = c; h[12].m[2].a = -c; h[12].m[3].a = c;
h[13].m[1].a = c; h[13].m[2].a = c; h[13].m[3].a = -c;
h[4].m[0].a = 2*clip_numeric.radius - 0*clip_numeric.center_x - a*clip_numeric.center_y - b*clip_numeric.center_z;
h[5].m[0].a = 2*clip_numeric.radius - b*clip_numeric.center_x - 0*clip_numeric.center_y - a*clip_numeric.center_z;
h[6].m[0].a = 2*clip_numeric.radius - a*clip_numeric.center_x - b*clip_numeric.center_y - 0*clip_numeric.center_z;
h[7].m[0].a = 2*clip_numeric.radius - 0*clip_numeric.center_x - a*clip_numeric.center_y + b*clip_numeric.center_z;
h[8].m[0].a = 2*clip_numeric.radius + b*clip_numeric.center_x - 0*clip_numeric.center_y - a*clip_numeric.center_z;
h[9].m[0].a = 2*clip_numeric.radius - a*clip_numeric.center_x + b*clip_numeric.center_y - 0*clip_numeric.center_z;
h[4].m[1].a = 0; h[4].m[2].a = a; h[4].m[3].a = b;
h[5].m[1].a = b; h[5].m[2].a = 0; h[5].m[3].a = a;
h[6].m[1].a = a; h[6].m[2].a = b; h[6].m[3].a = 0;
h[7].m[1].a = 0; h[7].m[2].a = a; h[7].m[3].a = -b;
h[8].m[1].a = -b; h[8].m[2].a = 0; h[8].m[3].a = a;
h[9].m[1].a = a; h[9].m[2].a = -b; h[9].m[3].a = 0;
h[14].m[0].a = 2*clip_numeric.radius + 0*clip_numeric.center_x + a*clip_numeric.center_y + b*clip_numeric.center_z;
h[15].m[0].a = 2*clip_numeric.radius + b*clip_numeric.center_x + 0*clip_numeric.center_y + a*clip_numeric.center_z;
h[16].m[0].a = 2*clip_numeric.radius + a*clip_numeric.center_x + b*clip_numeric.center_y + 0*clip_numeric.center_z;
h[17].m[0].a = 2*clip_numeric.radius + 0*clip_numeric.center_x + a*clip_numeric.center_y - b*clip_numeric.center_z;
h[18].m[0].a = 2*clip_numeric.radius - b*clip_numeric.center_x + 0*clip_numeric.center_y + a*clip_numeric.center_z;
h[19].m[0].a = 2*clip_numeric.radius + a*clip_numeric.center_x - b*clip_numeric.center_y + 0*clip_numeric.center_z;
h[14].m[1].a = 0; h[14].m[2].a = -a; h[14].m[3].a = -b;
h[15].m[1].a = -b; h[15].m[2].a = 0; h[15].m[3].a = -a;
h[16].m[1].a = -a; h[16].m[2].a = -b; h[16].m[3].a = 0;
h[17].m[1].a = 0; h[17].m[2].a = -a; h[17].m[3].a = b;
h[18].m[1].a = b; h[18].m[2].a = 0; h[18].m[3].a = -a;
h[19].m[1].a = -a; h[19].m[2].a = b; h[19].m[3].a = 0;
for( i=0; i<10; i++ )
{
p[i] = polyxyz_mult( &h[i],&h[i+10] );
polyxyz_collect( &p[i] );
polyxyz_set_degree( &p[i] );
polyxyz_adjust( &p[i] );
hp[i] = new hornerpolyxyz ( p[i] );
}
for( i=10; i<MAIN_CLIP_AMOUNT_NUM; i++ )
{
p[i].n = 0;
}
for( i=0; i<20; i++ )
{
delete_coeff_polyxyz( &h[i] );
}
}
syntax highlighted by Code2HTML, v. 0.9.1