Hole2 update with invert

/*
  JWildfire - an image and animation processor written in Java
  Copyright (C) 1995-2011 Andreas Maschke

  This is free software; you can redistribute it and/or modify it under the terms of the GNU Lesser
  General Public License as published by the Free Software Foundation; either version 2.1 of the
  License, or (at your option) any later version.

  This software 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
  Lesser General Public License for more details.

  You should have received a copy of the GNU Lesser General Public License along with this software;
  if not, write to the Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA
  02110-1301 USA, or see the FSF site: http://www.fsf.org.
 */
package org.jwildfire.create.tina.variation;

import static org.jwildfire.base.mathlib.MathLib.*;

import org.jwildfire.base.Tools;
import org.jwildfire.create.tina.base.XForm;
import org.jwildfire.create.tina.base.XYZPoint;

public class Hole2Func extends VariationFunc {
   private static final long serialVersionUID = 1L;
   private static final String PARAM_A = "a";
   private static final String PARAM_B = "b";
   private static final String PARAM_C = "c";
   private static final String PARAM_D = "d";
   public static final String PARAM_INSIDE = "inside";
   public static final String PARAM_SHAPE = "shape";
   public static final String PARAM_INVERT = "invert";
   private static final String[] paramNames = { PARAM_A, PARAM_B, PARAM_C, PARAM_D, PARAM_INSIDE, PARAM_SHAPE, PARAM_INVERT };

   private double a = 1.0;
   private double b = 2.0;
   private double c = 1.0;
   private double d = 1.0;

   private int inside = 0;
   private int shape = 0;
   private int invert = 0;

   @Override
   public void transform(FlameTransformationContext pContext, XForm pXForm,
         XYZPoint pAffineTP, XYZPoint pVarTP, double pAmount) {
      // Hole3 by Michael Faber, Brad Stefanov, and Rick Sidwell
	    if (invert == 0) {

	  double rhosq = pAffineTP.getPrecalcSumsq();
      double theta = pAffineTP.getPrecalcAtanYX() * d;
      double delta = pow(theta / M_PI + 1.0, a) * c;
      double r1 = 1;



      switch (shape) {
      case 0: // hole
          r1 = sqrt(rhosq) + delta;
          break;
      case 1:// hole1
          r1 = sqrt(rhosq + delta);
          break;
      case 2:// double hole
          r1 = sqrt(rhosq + sin(b * theta) + delta);
          break;
      case 3:// heart
          r1 = sqrt(rhosq + sin(theta) + delta);
          break;
      case 4:// heart2
          r1 = sqrt(rhosq + sqr(theta) - delta + 1);
          break;
      case 5:
          r1 = sqrt(rhosq + fabs(tan(theta)) + delta);
          break;
      case 6:
    	  r1 = sqrt(rhosq * (1 + sin(b * theta)) + delta);
    	  break;
      case 7:
    	  r1 = sqrt(rhosq + fabs(sin(0.5 * b * theta)) + delta);
    	  break;
      case 8:
    	  r1 = sqrt(rhosq + sin(M_PI * sin(b * theta)) + delta);
    	  break;
      case 9:
    	  r1 = sqrt(rhosq + (sin(b * theta) + sin(2 * b * theta + M_PI_2)) / 2 + delta);
    	  break;
      }

      if (inside != 0)
    	  r1 = pAmount / r1;
      else
    	  r1 = pAmount * r1;

      pVarTP.x += r1 * cos(theta);
      pVarTP.y += r1 * sin(theta);
      if (pContext.isPreserveZCoordinate()) {
          pVarTP.z += pAmount * pAffineTP.z;
      }
      }
	    else /*if (invert == 1)*/{
	  	  double rhosq = -pAffineTP.getPrecalcSumsq();
	      double theta = -pAffineTP.getPrecalcAtanYX() * d;
	      double delta = pow(theta / M_PI - 1.0, a) * c;
	      double r1 = 1;


	      switch (shape) {
	      case 0: // hole
	          r1 = sqrt(rhosq) + delta;
	          break;
	      case 1:// hole1
	          r1 = sqrt(rhosq + delta);
	          break;
	      case 2:// double hole
	          r1 = sqrt(rhosq + sin(b * theta) + delta);
	          break;
	      case 3:// heart
	          r1 = sqrt(rhosq + sin(theta) + delta);
	          break;
	      case 4:// heart2
	          r1 = sqrt(rhosq + sqr(theta) - delta + 1);
	          break;
	      case 5:
	          r1 = sqrt(rhosq + fabs(tan(theta)) + delta);
	          break;
	      case 6:
	    	  r1 = sqrt(rhosq * (1 + sin(b * theta)) + delta);
	    	  break;
	      case 7:
	    	  r1 = sqrt(rhosq + fabs(sin(0.5 * b * theta)) + delta);
	    	  break;
	      case 8:
	    	  r1 = sqrt(rhosq + sin(M_PI * sin(b * theta)) + delta);
	    	  break;
	      case 9:
	    	  r1 = sqrt(rhosq + (sin(b * theta) + sin(2 * b * theta + M_PI_2)) / 2 + delta);
	    	  break;
	      }

	      if (inside != 0)
	    	  r1 = pAmount / r1;
	      else
	    	  r1 = pAmount * r1;

	      pVarTP.x += r1 * cos(theta);
	      pVarTP.y += r1 * sin(theta);
	      if (pContext.isPreserveZCoordinate()) {
	          pVarTP.z += pAmount * pAffineTP.z;
	      }
	    }

   }

   @Override
   public String[] getParameterNames() {
      return paramNames;
   }

   @Override
   public Object[] getParameterValues() {
      return new Object[] { a, b, c, d, inside, shape, invert };
   }

   @Override
   public void setParameter(String pName, double pValue) {
      if (PARAM_A.equalsIgnoreCase(pName))
         a = pValue;
      else if (PARAM_B.equalsIgnoreCase(pName))
         b = pValue;
      else if (PARAM_C.equalsIgnoreCase(pName))
         c = pValue;
      else if (PARAM_D.equalsIgnoreCase(pName))
         d = pValue;
      else if (PARAM_INSIDE.equalsIgnoreCase(pName))
         inside = limitIntVal(Tools.FTOI(pValue), 0, 1);

      else if (PARAM_SHAPE.equalsIgnoreCase(pName))
         shape = limitIntVal(Tools.FTOI(pValue), 0, 9);
      else if (PARAM_INVERT.equalsIgnoreCase(pName))
      invert = Tools.FTOI(pValue);

    else
         throw new IllegalArgumentException(pName);
   }

   @Override
   public String getName() {
      return "hole2";
   }

}