/*
 *   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 "NewClipSphereParallel.h"
#include "def.h"

int NewClipSphereParallel::clip_user_y (double uy)
{
	// return  fabs( uy - cy ) <= r;
	if (fabs( uy - cy ) > r)
		return FALSE;

	A0 = uy - cy;
	A0 *= A0;
	return TRUE;
}

int NewClipSphereParallel::clip_user_xy (double ux, double uy, double &zmin, double &zmax)
{
	double  g = ux - cx;
	double  h = g*g + A0;

	if( h <= rsq ) {
		h = sqrt( rsq - h );
		zmin = cz - h;
		zmax = cz + h;
		// fixme
                // mach ich doch ...
		zmin = max( zmin,cb /*-10.0*/ /*cb*/);
		zmax = min( zmax,cf /*10.0*/ /*cf*/ );
		
		return  TRUE;
	} else {
		return  FALSE;
	}
}

int NewClipSphereParallel::clip_user_xyz (double uz)
{
    return  TRUE;
}