// Mark routines and markers written by Philippe Ivaldi.
// http://www.piprime.fr/

marker operator * (transform T, marker m)
{
  marker M=new marker;
  M.f=T*m.f;
  M.above=m.above;
  M.markroutine=m.markroutine;
  return M;
}

// Add n frames f midway (in arclength) between n+1 uniformly spaced marks.
markroutine markinterval(int n=1, frame f, bool rotated=false)
{
  return new void(picture pic=currentpicture, frame mark, path g) {
    markuniform(n+1,rotated)(pic,mark,g);
    markuniform(centered=true,n,rotated)(pic,f,g);
  };
}

// Return a frame containing n copies of the path g shifted by space
// drawn with pen p.
frame duplicate(path g, int n=1, pair space=0, pen p=currentpen)
{
  if(space == 0) space=dotsize(p);
  frame f;
  int pos=0;
  int sign=1;
  int m=(n+1) % 2;
  for(int i=1; i <= n; ++i) {
    draw(f,shift(space*(pos-0.5*m))*g,p);
    pos += i*sign;
    sign *= -1;
  }
  return f;
}

real tildemarksizefactor=5;
real tildemarksize(pen p=currentpen)
{
  static real golden=(1+sqrt(5))/2;
  return (1mm+tildemarksizefactor*sqrt(linewidth(p)))/golden;
}
frame tildeframe(int n=1, real size=0, pair space=0,
                 real angle=0, pair offset=0, pen p=currentpen)
{
  size=(size == 0) ? tildemarksize(p) : size;
  space=(space == 0) ? 1.5*size : space;
  path g=yscale(1.25)*((-1.5,-0.5)..(-0.75,0.5)..(0,0)..(0.75,-0.5)..(1.5,0.5));
  return duplicate(shift(offset)*rotate(angle)*scale(size)*g,n,space,p);
}

frame tildeframe=tildeframe();

real stickmarkspacefactor=4;
real stickmarksizefactor=10;
real stickmarksize(pen p=currentpen)
{
  return 1mm+stickmarksizefactor*sqrt(linewidth(p));
}
real stickmarkspace(pen p=currentpen)
{
  return stickmarkspacefactor*sqrt(linewidth(p));
}
frame stickframe(int n=1, real size=0, pair space=0, real angle=0,
                 pair offset=0, pen p=currentpen)
{
  if(size == 0) size=stickmarksize(p);
  if(space == 0) space=stickmarkspace(p);
  return duplicate(shift(offset)*rotate(angle)*scale(0.5*size)*(N--S),n,
                   space,p);
}

frame stickframe=stickframe();

real circlemarkradiusfactor=stickmarksizefactor/2;
real circlemarkradius(pen p=currentpen)
{
  static real golden=(1+sqrt(5))/2;
  return (1mm+circlemarkradiusfactor*sqrt(linewidth(p)))/golden;
}
real barmarksizefactor=stickmarksizefactor;
real barmarksize(pen p=currentpen)
{
  return 1mm+barmarksizefactor*sqrt(linewidth(p));
}
frame circlebarframe(int n=1, real barsize=0,
                     real radius=0,real angle=0,
                     pair offset=0, pen p=currentpen,
                     filltype filltype=NoFill, bool above=false)
{
  if(barsize == 0) barsize=barmarksize(p);
  if(radius == 0) radius=circlemarkradius(p);
  frame opic;
  path g=circle(offset,radius);
  frame f=stickframe(n,barsize,space=2*radius/(n+1),angle,offset,p);
  if(above) {
    add(opic,f);
    filltype.fill(opic,g,p);
  } else {
    filltype.fill(opic,g,p);
    add(opic,f);
  }
  return opic;
}

real crossmarksizefactor=5;
real crossmarksize(pen p=currentpen)
{
  return 1mm+crossmarksizefactor*sqrt(linewidth(p));
}
frame crossframe(int n=3, real size=0, pair space=0,
                 real angle=0, pair offset=0, pen p=currentpen)
{
  if(size == 0) size=crossmarksize(p);
  frame opic;
  draw(opic,shift(offset)*rotate(angle)*scale(size)*cross(n),p);
  return opic;
}

real markanglespacefactor=4;
real markangleradiusfactor=8;
real markangleradius(pen p=currentpen)
{
  return 8mm+markangleradiusfactor*sqrt(linewidth(p));
}
real markangleradius=markangleradius();
real markanglespace(pen p=currentpen)
{
  return markanglespacefactor*sqrt(linewidth(p));
}
real markanglespace=markanglespace();
// Mark the oriented angle AOB counterclockwise with optional Label, arrows, and markers.
// With radius < 0, AOB-2pi is marked clockwise.
void markangle(picture pic=currentpicture, Label L="",
               int n=1, real radius=0, real space=0,
               pair A, pair O, pair B, arrowbar arrow=None,
               pen p=currentpen, filltype filltype=NoFill,
               margin margin=NoMargin, marker marker=nomarker)
{
  if(space == 0) space=markanglespace(p);
  if(radius == 0) radius=markangleradius(p);
  picture lpic,phantom;
  frame ff;
  path lpth;
  p=squarecap+p;
  pair OB=unit(B-O), OA=unit(A-O);
  real xoa=degrees(OA,false);
  real gle=degrees(acos(dot(OA,OB)));
  if((conj(OA)*OB).y < 0) gle *= -1;
  bool ccw=radius > 0;
  if(!ccw) radius=-radius;
  bool drawarrow = !arrow(phantom,arc((0,0),radius,xoa,xoa+gle,ccw),p,margin);
  if(drawarrow && margin == NoMargin) margin=TrueMargin(0,0.5linewidth(p));
  if(filltype != NoFill) {
    lpth=margin(arc((0,0),radius+(n-1)*space,xoa,xoa+gle,ccw),p).g;
    pair p0=relpoint(lpth,0), p1=relpoint(lpth,1);
    pair ac=p0-p0-A+O, bd=p1-p1-B+O, det=(conj(ac)*bd).y;
    pair op=(det == 0) ? O : p0+(conj(p1-p0)*bd).y*ac/det;
    filltype.fill(ff,op--lpth--relpoint(lpth,1)--cycle,p);
    add(lpic,ff);
  }
  for(int i=0; i < n; ++i) {
    lpth=margin(arc((0,0),radius+i*space,xoa,xoa+gle,ccw),p).g;
    draw(lpic,lpth,p=p,arrow=arrow,margin=NoMargin,marker=marker);
  }
  Label lL=L.copy();
  real position=lL.position.position.x;
  if(lL.defaultposition) {lL.position.relative=true; position=0.5;}
  if(lL.position.relative) position=reltime(lpth,position);
  if(lL.align.default) {
    lL.align.relative=true;
    lL.align.dir=unit(point(lpth,position));
  }
  label(lpic,lL,point(lpth,position),align=NoAlign, p=p);
  add(pic,lpic,O);
}

marker StickIntervalMarker(int i=2, int n=1, real size=0, real space=0,
                           real angle=0, pair offset=0, bool rotated=true,
                           pen p=currentpen, frame uniform=newframe,
                           bool above=true)
{
  return marker(uniform,markinterval(i,stickframe(n,size,space,angle,offset,p),
                                     rotated),above);
}


marker CrossIntervalMarker(int i=2, int n=3, real size=0, real space=0,
                           real angle=0, pair offset=0, bool rotated=true,
                           pen p=currentpen, frame uniform=newframe,
                           bool above=true)
{
  return marker(uniform,markinterval(i,crossframe(n,size,space,angle,offset,p),
                                     rotated=rotated),above);
}

marker CircleBarIntervalMarker(int i=2, int n=1, real barsize=0, real radius=0,
                               real angle=0, pair offset=0, bool rotated=true,
                               pen p=currentpen, filltype filltype=NoFill,
                               bool circleabove=false, frame uniform=newframe,
                               bool above=true)
{
  return marker(uniform,markinterval(i,circlebarframe(n,barsize,radius,angle,
                                                      offset,p,filltype,
                                                      circleabove),
                                     rotated),above);
}

marker TildeIntervalMarker(int i=2, int n=1, real size=0, real space=0,
                           real angle=0, pair offset=0, bool rotated=true,
                           pen p=currentpen, frame uniform=newframe,
                           bool above=true)
{
  return marker(uniform,markinterval(i,tildeframe(n,size,space,angle,offset,p),
                                     rotated),above);
}