/*****************************************************************************
 * Copyright 2005 Kevin Hobbs <hobbsk@ohiou.edu>                             *
 *                                                                           *
 * This file is part of Balloon.                                             *
 *                                                                           *
 * Balloon 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.                                       *
 *                                                                           *
 * Balloon 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 Balloon; if not, write to the Free Software                    *
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA *
 *****************************************************************************/


#include <stddef.h>
#include <math.h>
#include "balloon.h"

int in_out( const double * x, double * y, const size_t ndims)
{
	double sum = 0;
	double r = 6;
	double r2 = 2.0;
	double at = 3.0;
	double sum2 = 0;

	size_t i;

	for (i = 0; i < ndims; i++)
		sum += x[i] * x[i];
	sum = sqrt(sum);

	for (i = 0; i < ndims; i++)
		sum2 += (x[i] - at) * (x[i] - at);
	sum2 = sqrt(sum2);
	
	*y = sum;

	if ((sum < r)  && (sum2 > r2))
		return FUNCTION_INSIDE;

	return FUNCTION_OUTSIDE;
}