/*****************************************************************************
 * 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 *
 *****************************************************************************/


#define FUNCTION_INSIDE 0
#define FUNCTION_OUTSIDE 1

#define STILL_LONG 1
#include <gsl/gsl_rng.h>

struct point ** 
init_balloon(
		const double * x0, 
		const double * s, 
		const size_t ndims,
		size_t * id,
		gsl_rng * r);

int inflate( 
		struct point * p, 
		double * s, 
		const double * s0,
		const size_t ndims,
		gsl_rng * r);

void direction( const struct point * p,
		double * s, 
		const size_t ndims );

void mid_point( const struct point * p1,
		const struct point * p2,
		double * x,
		const size_t ndims );

int long_segment(
		const struct point ** balloon, 
		const size_t num_points,
		const double * s0,
		const size_t ndims,
		struct point ** p1, 
		struct point ** p2 );
		
		
int in_out( const double * x, double * y, const size_t ndims);