// Flower generation - Jeff Clark, Aug 23rd, 2006
// This application is based on the one found online at http://rpi.edu/~mcdonk/code/flower/ .
// The source is hacked together but has some interesting ideas.
// Set bCreateOne = true; in the setup function to get a flower when first run.
// Uncomment the MousePressed() method to make a new one on a click.

int maxle, maxra, lightd, lvar, la, darkd, dvar, da, maxbright, brightvar, sa, basehu, inithuvar, huvar;
float minbend, edgevar, maxedgevar, edgemem, stepinc, wot, hot, steple, edgememp1, petalsperlevel, dHairSparseness;
float edgeWidth;
float fExternalHue;
boolean bCreateOne;
float fHue1, fHue2, fHue3;

float[] toRect(float ra, float th)
{
float[] ou={sin(th)*ra,cos(th)*ra};
  return ou;
}

void setup()
{
  int iMrg = 30; 
  int iHgt = 20; 
  framerate(5); 

  randomSeed(0);
  background(255);
  size(400,400);
  setParms();
  bCreateOne = false;
  fExternalHue = 128;
  fHue1 = 0;
  fHue2 = 160;
  fHue3 = 40;
}

void setParms()
{
  wot=width/2;
  hot=height/2;

  colorMode(HSB);

  maxle=(int)random(3,10); //layers of petals
  maxra=(int)random(120,160); //initial radius

  edgevar=random(50,200); //velocity on random edge variance, low values -> smoother edges
  maxedgevar=random(10,60); //max edge deviation
  edgemem=random(300,600); //memory for the edge randomness, low values -> really jagged but depends on edgevar

  lightd=(int)random(5,15); //light distance
  lvar=(int)random(5,20); //light randomness variance
  la=(int)random(5,20); //light alpha multiplier
  darkd=(int)random(10,20); //(int)random(10,20); //dark distance
  dvar=(int)random(10,20); //dark randomness variance
  da=(int)random(8,18); //dark alpha multiplier
  dHairSparseness = 0;//.75;//random(0,1); // lower values - > less hair

  maxbright=240;//(int)random(200,240); //maximum brightness of petals
  brightvar=0;//120;//150;//(int)random(50,150); //maxbright-brightvar is the brightness at the center
  sa=240;//(int)random(180,255); //saturation
  basehu=128; //color center
  inithuvar=128; //color variance
  huvar=(int)random(5,25); //step variance

  stepinc=.0015;

  steple=maxra/maxle;
  edgememp1=edgemem+1; //for speed
  
  edgeWidth = random(.80,1);
  if(random(0,1) < .8)
    edgeWidth = 1;
  edgeWidth = .95;  
  //minbend=random(.05,.80); //random(.65,.80); //"bend" on petals (1=circles)
  minbend=random(.00,.70); //"bend" on petals (1=circles)
}

void CreateRose(float fh1, float fh2, float fh3, float fSat, float fBrightness, float fHairiness, float fEdgeVariance, float fMinBend, int iLevels, int iRandomseed)
{
  //background(iback);
  randomSeed(iRandomseed);
  setParms();
  fHue1 = fh1;
  fHue2 = fh2;
  fHue3 = fh3;
  sa = (int)fSat;
  brightvar = (int) fBrightness;
  maxle = iLevels;
  dHairSparseness = fHairiness;
  edgevar = 50 + 150 * fEdgeVariance;
  minbend = fMinBend;
  bCreateOne = true;
}

void CreateRose(int iback)
{
  //background(iback);
  setParms();
  fHue1 = iback;
  fHue2 = fHue3 = fHue1;
  bCreateOne = true;
}

void draw()
{
  //if(IsNewText())
  if(bCreateOne)
   {
    background(255);
    rose();
    bCreateOne = false;
  }
}

void mousePressed()
{
  /*
  setParms();
  fHue1 = nonGreenRandomHue();
  fHue2 = nonGreenRandomHue();
  fHue3 = fHue1;
  maxle=(int)random(3,5); //layers of petals
  bCreateOne = true;
  */
}

float nonGreenRandomHue()
{
  float hu = random(0,255);
  // eliminate green
  while(hu > 42 && hu < 110)
  {
    hu = random(0,255);
  }
  return hu;
}

/*
void mouseReleased()
{ 
  noLoop(); 
} 
*/

void rose()
{
  translate(wot,hot);

/*
  float hu = random(basehu-inithuvar,basehu+inithuvar);
  // eliminate green
  while(hu > 42 && hu < 110)
  {
    hu = random(basehu-inithuvar,basehu+inithuvar);;
  }
*/
  float hu = fExternalHue;
  for(float le=0;le<maxle;le++)
  {
    // le == 0 is outermost level
    if(le % 3 != 2 )
    {
      hu = fHue1;
    }
    else
    {
      hu = fHue2;
    }
   float fHairSparseness = dHairSparseness;
   /*
    if(le == 0)
    {
      fHairSparseness = dHairSparseness / 10;
    }
    */
    //normally globals:
    petalsperlevel=(int)random(2,5)+(float).5; //odd numbers /2
    //if(petalsperlevel % 2 == 0) petalsperlevel -=1;
    //if(petalsperlevel < 2) petalsperlevel = 2;
  
    float completed=(le/maxle);
    float rad=maxra-(le*steple);
    float ap=(minbend);//+(le/50));
    float a=ap*rad;
    float b=(1-ap)*rad;
    float an=completed*TWO_PI+random(1);//1);

    float edgepos=0;

    //hu=(hu+random(-huvar/2,huvar))%255; //a;ways move forward
    float br=maxbright-(completed*brightvar);
    float shadea=(1-completed);
    shadea*=shadea;//darker on the outside
    
    float huEdge = (hu+random(-huvar/2,huvar))%255;//(int)random(64,192);
    println("Level:" + le + " Hue:" + hu + " sa=" + sa + " br=" + br);

    float[] tr,trl,trd,tr2,tr3;
    for(float t=0;t<2*TWO_PI;t+=stepinc)
    {

      float r=a+b*cos(petalsperlevel*t+an)+edgepos;
      tr3=toRect(r*1.005,t);
      tr2=toRect(r*edgeWidth,t);
      tr =toRect(r,t);
      if(r>a)
      {
        //println(r);
        stroke(hu,sa,br,100);
        //stroke(0,0,240,100);
        smooth();
        if(true)//le == maxle - 1)
          blendedline(0,0,tr[0],tr[1],hu,fHue2,20,br,(int)r/5); //40);
        else
          line(0,0,tr[0],tr[1]);
        if(edgeWidth != 1 && le % 3 == 0)
        {
          stroke(fHue3,sa,br,190);
          line(tr[0],tr[1],tr2[0],tr2[1]);
        }
        stroke(hu,sa,60,140);
        line(tr[0],tr[1],tr3[0],tr3[1]);
        noSmooth();

        trl=toRect(r-lightd,t);
        stroke(0,0,255,shadea*la);
        line(tr[0],tr[1],trl[0]+random(-lvar,lvar),trl[1]+random(-lvar,lvar));

        if(random(0,1) < fHairSparseness)
        {
          trd=toRect(r+darkd,t);
          //stroke(0,0,0,shadea*da);
          stroke(hu,sa,br/2,shadea*da);
          line(tr[0],tr[1],trd[0]+random(-dvar,dvar),trd[1]+random(-dvar,dvar));
        }
      }
      edgepos=(edgepos*edgemem+random(-edgevar,edgevar))/(edgememp1);
      edgepos=constrain(edgepos,-maxedgevar,maxedgevar);
    }
  }
}

void blendedline(float x1, float y1, float x2, float y2,float fhue1, float fhue2, float fbr1, float fbr2, int iSegments)
{
  if(iSegments == 0) return; // error
  float[] ftr1;
  float[] ftr2;
  ftr1 = toRect(100,.5);
  ftr2 = toRect(100,.5);
  float fSegx = (x2 - x1) / iSegments;
  float fSegy = (y2 - y1) / iSegments;
  float fBrDelta = (fbr2 - fbr1) / iSegments;
  //if(fhue1 > fhue2)
  //{
  //  float ftemp = fhue2;
  //  fhue2 = fhue1;
  //  fhue1 = ftemp;
  //}
  float fHueDelta = (fhue2 - fhue1) / (float)iSegments;
  for(int iSeg = 1; iSeg <= iSegments; iSeg++)
  {
    ftr1[0] = x1 + (iSeg - 1) * fSegx;    
    ftr2[0] = x1 + (iSeg)     * fSegx;    
    ftr1[1] = y1 + (iSeg - 1) * fSegy;    
    ftr2[1] = y1 + (iSeg)     * fSegy; 
    float fbrCur = fbr1 + iSeg * fBrDelta;
    //float fHueCur = fhue1 + iSeg * fHueDelta;
    stroke(fhue1,sa,fbrCur,100);
    line(ftr1[0],ftr1[1],ftr2[0],ftr2[1]);
  }
}