/*
    Apophysis Plugin

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

// Must define this structure before we include apoplugin.h
typedef struct
{
		double hole_a;
		int hole_inside;
} Variables;

#define _USE_MATH_DEFINES
#include "apoplugin.h"

// Set the name of this plugin
APO_PLUGIN("hole");

// Define the Variables
APO_VARIABLES(
	VAR_REAL(hole_a, 1.0),
	VAR_INTEGER_RANGE( hole_inside, 0, 1, 0)
);

// You must call the argument "vp".
int PluginVarPrepare(Variation* vp)
{
    return TRUE; // Always return TRUE.
}

// You must call the argument "vp".
int PluginVarCalc(Variation* vp)
{

/*	double a = atan2(FTy, FTx);
	double r = VVAR * sqrt(sqr(FTx) + sqr(FTy) + pow( a/M_PI + 1.0, VAR(hole_a)));
	double s, c;
	fsincos( a, &s, &c);
	
	FPx += r * c;
	FPy += r * s;
	*/
	
	double a = atan2(FTy, FTx);
	double delta = pow( a/M_PI + 1.0, VAR(hole_a));
	double r;
	
	if(VAR(hole_inside))
	{
       r = VVAR * delta / ( sqrt(sqr(FTx) + sqr(FTy)) + delta );
    }
    else
    {
	   r = VVAR * sqrt(sqr(FTx) + sqr(FTy)) + delta ;
    }
	double s, c;
	fsincos( a, &s, &c);

	FPx += r * c;
	FPy += r * s;
	

/* heart 2
	double a = atan2(FTy, FTx);
	double r = VVAR * sqrt(sqr(FTx) + sqr(FTy) + sqr(a) + 0.0 );
	double s, c;
	fsincos( a, &s, &c);
	
	FPx += r * c;
	FPy += r * s;


	double a = atan2(FTy, FTx);
	double r = VVAR * sqrt(sqr(FTx) + sqr(FTy) + fabs(tan(a)) );
	double s, c;
	fsincos( a, &s, &c);
	
	FPx += r * c;
	FPy += r * s;


/* heart
	double a = atan2(FTy, FTx);
	double r = VVAR * sqrt(sqr(FTx) + sqr(FTy) + sin(a) + 1.0 );
	double s, c;
	fsincos( a, &s, &c);
	
	FPx += r * c;
	FPy += r * s;
*/
/* two hole
	double a = atan2(FTy, FTx);
	double r = VVAR * sqrt(sqr(FTx) + sqr(FTy) + sin(2.0 * a) + 1.0 );
	double s, c;
	fsincos( a, &s, &c);
	
	FPx += r * c;
	FPy += r * s;
*/
    return TRUE;
}