annotate SFDemoUnit.cpp @ 1:f3fd4e19cec0 tip

first binary upload

/*
    Harmonic Visualiser

    An audio file viewer and editor.
    Centre for Digital Music, Queen Mary, University of London.
    This file copyright 2011 Wen Xue.

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

#include <vcl.h>
#pragma hdrstop

#include <Math.hpp>
#include <math.h>
#include "SFDemoUnit.h"
#include "Matrix.h"
#include "Unit1.h"
#include "splines.h"
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TSFDemoForm *SFDemoForm;
//---------------------------------------------------------------------------
const double log10e=Log10(exp(1));
const Bdw=5, Bdwc=2; //border width of charts
const ftrr=5; //filter control point radius


__fastcall TSFDemoForm::TSFDemoForm(TComponent* Owner)
	: TForm(Owner)
{
	ForceUpdate=false;

	double updb=updbEdit->Text.ToDouble(), downdb=downdbEdit->Text.ToDouble();
	dbrange=updb-downdb;

	SUThread=0;
	pThread=0;
	memset(ThreadList, 0, sizeof(TSUThread*)*ThreadCaps);

  WaveAudio1=new TWaveAudio(NULL);
  WaveAudio1->UseMemoryStream=true;
	WaveAudio1->Channels=1;
  WaveView1=new TWaveView(NULL);
  WaveView1->Parent=Panel2;
  WaveView1->Align=alClient;
  WaveView1->WaveAudio=WaveAudio1;
  WaveView1->OnGetOpMode=WaveView1OpMode;
  WaveView1->CreatePanes(1, 1);
  WaveView1->SetContent(0, 0, 0, 1);
  WaveView1->ClickFocus=true;
  WaveView1->DefaultPopupMenu=false;
  WaveView1->CustomInfo=WaveView1CustomInfo;
  WaveView1->InfoLeft=5; WaveView1->InfoTop=5;

  WaveAudio2=new TWaveAudio(NULL);
  WaveAudio2->UseMemoryStream=true;
  WaveAudio2->Channels=1;
  WaveView2=new TWaveView(NULL);
  WaveView2->Parent=Panel3;
  WaveView2->Align=alClient;
  WaveView2->WaveAudio=WaveAudio2;
  WaveView2->CreatePanes(1, 1);
  WaveView2->SetContent(0, 0, 0, 1);
  WaveView2->ClickFocus=true;
  WaveView2->DefaultPopupMenu=false;
  WaveView2->CustomInfo=WaveView1CustomInfo;
  WaveView2->InfoLeft=5; WaveView2->InfoTop=5;
	WaveView2->LoopPlay=CheckBox1->Checked;		


  HS=0;
  cyclefrs=0;
	cyclefs=0;

	MethodListBox->ItemIndex=0;

	analyze2=false;

  AnsiString Path=ExtractFilePath(Application->ExeName)+"*.sf";
  TSearchRec sr;
  if (FindFirst(Path, faAnyFile, sr)==0)
  {
    do
    {
      TmplListBox->Items->Add(sr.Name);
    }
    while (FindNext(sr) == 0);
    FindClose(sr);
  }
  Wave1=0;
  CurrentP=-1;

  memset(&SF, 0, sizeof(TSF));
  memset(&SF0, 0, sizeof(TSF));

  datain=0;
}

__fastcall TSFDemoForm::~TSFDemoForm()
{
	delete WaveAudio1;
	delete WaveAudio2;
  delete WaveView1;
  delete WaveView2;
  delete HS;
  delete[] cyclefrs;
	delete[] cyclefs;
	for (int i=0; i<ThreadCaps; i++) if (ThreadList[i]) delete ThreadList[i];
  delete[] Wave1;
  delete[] datain;
}
//---------------------------------------------------------------------------

void TSFDemoForm::Copydata()
{
  int dst, den;
  double* xrec;
  if (HS)
  {
    xrec=SynthesisHSp(HS, dst, den);
    if (dst<0) dst=0;

    delete[] datain;
    datain=new __int16[den-dst];

    DoubleToInt(datain, sizeof(__int16), xrec, den-dst);

    free8(xrec);
  }
}


void TSFDemoForm::FindNote()
{
}

void TSFDemoForm::Reset()
{
	if (!HS) return;
	if (HS->M<=0 || HS->Fr<=0) return;

	if (MethodListBox->ItemIndex<0) MethodListBox->ItemIndex=0;

	double sps=WaveView1->SamplesPerSec;
	double h=WaveView1->SpecOffst;
	int FSMode=MethodListBox->ItemIndex;
	double FSF=FEdit->Text.ToDouble();
	int FSFScale=FScaleCombo->ItemIndex;
	if (FSFScale==0) FSF/=sps;
	else FSF/=Amel*log(1+sps/700);
	double FStheta=ThetaEdit->Text.ToDouble();

	AnalyzeSF_1(*HS, SF, sps, h);
	AnalyzeSF_2(*HS, SF, cyclefrs, cyclefs, sps, &cyclefrcount, FSMode, FSF, FSFScale, FStheta);

  SaveSF();
  SaveSF0();

  UpdateDisplay();
}

void TSFDemoForm::Synthesize()
{
  int dst, den;
  double* xrec;
  if (HS)
  {
      double t=GetTickCount(); 
      xrec=SynthesisHS(HS, dst, den);
      t=GetTickCount()-t;
      Label13->Caption=t;
	}
  int bips=WaveAudio1->BitsPerSample;
	double max=(1<<(bips-1)); for (int i=dst; i<den; i++) if (xrec[i]>max) max=xrec[i]; else if (xrec[i]<-max) max=-xrec[i];
	if (max>(1<<(bips-1)))
	{
		max=(1<<(bips-1))/max;
		for (int i=dst; i<den; i++) xrec[i]*=max;
	}
  int bps=bips/8;
  WaveAudio2->Clear(NULL);
  WaveAudio2->GetWaveProperties(WaveAudio1);
  if (HS)
  {
    if (dst<0) dst=0;
    void* data=new char[(den-dst)*bps];
    DoubleToInt(data, bps, xrec, den-dst);

    WaveAudio2->Clear(NULL);
    WaveAudio2->WriteSamples(data, den-dst);
    free8(xrec);
    delete[] data;
  }
}

  //for frequency modulator shape
  void DrawsBarMap(Graphics::TBitmap* Bitmap, TRect Rect, int Count, double* data, bool res, int* KX1=0, int* KX2=0)
  {
    TCanvas* Canvas=Bitmap->Canvas;
    Bitmap->Width=Rect.Width();
    Bitmap->Height=Rect.Height();

    TFont* F=Canvas->Font; F->Color=clWhite; F->Height=12; F->Name="Ariel";
    int fh=Canvas->TextHeight("0");
    int DrawCount=res?(Count+1):Count;

    double XX=Rect.right-Bdw*2, YY=Rect.bottom-Bdw-Bdwc-fh;
    Canvas->Brush->Color=clBlack; Canvas->Brush->Style=bsSolid; Canvas->FillRect(Rect);
    Canvas->Brush->Color=clAqua; Canvas->Pen->Color=clAqua; Canvas->Brush->Style=bsSolid;
		int empspace=XX/(DrawCount+1)/4;
    double r=0; for (int i=0; i<Count; i++) r+=data[i]*data[i];

    double sqrtr=sqrt(r);
    double res_amp=(r>1)?0:sqrt(1-r);
    for (int i=0; i<DrawCount; i++)
    {
      double ldata;
      if (res)
      {
        if (i==0) ldata=data[i];
        else if (i<Count) ldata=data[i];
        else ldata=res_amp;
      }
      else ldata=data[i]/sqrtr;
      int X=floor(Bdw+XX*(i+1)/(DrawCount+1)+0.5);
      int Y=floor(Bdw+YY*(1-ldata)+0.5);
      if (empspace>=1)
      {
        Canvas->Brush->Style=bsSolid; Canvas->Brush->Color=TColor(0xFFFF00);
        Canvas->Rectangle(X-empspace, Y, X+empspace, Bdw+YY);
        if (KX1) KX1[i]=X-empspace, KX2[i]=X+empspace-1;
      }
      else
      {
        Canvas->MoveTo(X, Y); Canvas->LineTo(X, Bdw+YY);
        if (KX1) KX1[i]=X-1, KX2[i]=X+1;
      }
      AnsiString S=i; if (i>=Count) S="R"; int fw=Canvas->TextWidth(S);
      Canvas->Brush->Style=bsClear;
      Canvas->TextOut(X-fw/2, Bdw+YY, S);
      S.sprintf("%.3g", 100*ldata*ldata); if (S[1]=='0') S=S.SubString(2, S.Length()-1); fw=Canvas->TextWidth(S);
      if (Y-Bdw>fh) Canvas->TextOut(X-fw/2, Y-fh, S);
      else Canvas->TextOut(X+empspace+1, Y, S);
    }
    Canvas->Brush->Color=clBlack; Canvas->Brush->Style=bsClear; Canvas->Pen->Color=clWhite;
    Canvas->Rectangle(Bdw, Bdw, Bdw+XX, Bdw+YY);
  }

  TColor RotateColors[10]={clWhite, clAqua, clLime, clYellow, clRed, clTeal, clGray, clGreen, clFuchsia, clBlue};
  //draw A-F distribution
	void DrawAF(Graphics::TBitmap* Bitmap, TRect Rect, double xst, double xen, double dbst, double dben, int M, int frst, int fren, atom** Partials, bool MA, bool bars, double sps, int TextHeight=12)
	{
		double log10e=1.0/log(10);
		Bitmap->Width=Rect.Width();
		Bitmap->Height=Rect.Height();
		TCanvas* Canvas=Bitmap->Canvas;
		double idbrange=1.0/(dbst-dben);
		TFont* F=Canvas->Font; F->Color=clWhite; F->Height=TextHeight; F->Name="Ariel";
		int fh=Canvas->TextHeight("0");

		double XX=Rect.right-Bdw*2, YY=Rect.bottom-Bdw-Bdwc-fh;
		Canvas->Brush->Color=clBlack; Canvas->Brush->Style=bsSolid; Canvas->FillRect(Rect);

		Canvas->Pen->Color=clSilver; Canvas->MoveTo(Bdw, Bdw+YY*0.5); Canvas->LineTo(Bdw+XX, Bdw+YY*0.5);

		double *x=new double[(fren-frst)*4], *y=&x[fren-frst], *newa=&x[(fren-frst)*2], *neww=&x[(fren-frst)*3];

		for (int m=0; m<M; m++)
		{
			Canvas->Pen->Color=RotateColors[m%10]; Canvas->Brush->Color=RotateColors[m%10];

			int Sc=0;
			bool visible=false;
			atom* Partialsm=Partials[m];
			for (int fr=frst; fr<fren; fr++)
      {
        if (Partialsm[fr].f>0 && Partialsm[fr].a>0)
        {
          x[Sc]=Partials[m][fr].f;
					y[Sc]=20*(Log10(Partialsm[fr].a));
          if (x[Sc]<xen) visible=true;
          Sc++;
        }
      }
      if (!visible) break;

      if (MA)
      {
        double fw;
        for (int c=0; c<Sc; c++) newa[c]=y[c], neww[c]=1;
        for (int c=0; c<Sc; c++)
        {
          fw=0.002; if (fw>x[c]*0.05) fw=x[c]*0.05;

          for (int c1=0; c1<c; c1++)
          {
            double df=fabs(x[c]-x[c1]);
            if (df<fw)
            {
              double w=0.5+0.5*cos(M_PI*df/fw);
              newa[c]+=y[c1]*w; neww[c]+=w;
              newa[c1]+=y[c]*w; neww[c1]+=w;
            }
          }
        }
        for (int c=0; c<Sc; c++) y[c]=newa[c]/neww[c];
      }

      for (int c=0; c<Sc; c++)
      {
        if (x[c]<xst || x[c]>xen) continue;
				double ldata=(y[c]-dben)*idbrange;
        int X=floor(Bdw+XX*(x[c]-xst)/(xen-xst)+0.5);
        int Y=floor(Bdw+YY*(1-ldata)+0.5);
        if (bars) {Canvas->MoveTo(X, Y); Canvas->LineTo(X, Bdw+YY);}
        else {Canvas->MoveTo(X-1, Y); Canvas->LineTo(X+2, Y); Canvas->MoveTo(X, Y-1); Canvas->LineTo(X, Y+2);}
      }
    }

		AnsiString S=AnsiString().sprintf("%.2gkhz", sps*xst/1000); Canvas->Brush->Style=bsClear; Canvas->TextOut(Bdw, Bdw+YY, S);
		S=AnsiString().sprintf("%.2gkhz", sps*(xst+xen)/2000); int fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+(XX-fw)/2, Bdw+YY, S);
		S=AnsiString().sprintf("%.2gkhz", sps*xen/1000); fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+XX-fw, Bdw+YY, S);
    S.sprintf("%gdB", dbst); fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+XX-fw, Bdw, S);
    S.sprintf("%gdB", (dbst+dben)/2); fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+XX-fw, Bdw+(YY-fh)/2, S);
    S.sprintf("%gdB", dben); fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+XX-fw, Bdw+YY-fh, S);

    Canvas->Brush->Color=clBlack; Canvas->Brush->Style=bsClear; Canvas->Pen->Color=clWhite;
    Canvas->Rectangle(Bdw, Bdw, Bdw+XX, Bdw+YY);

    delete[] x;

  }   //*/
	//for source-filter results
	void DrawSFr(Graphics::TBitmap* Bitmap, TRect Rect, double xst, double xen, double dbst, double dben, int M, int frst, int fren, atom** Partials, double* logA0C, double* b, double F, int K, double* h, int FScaleMode, double Fs, int TextHeight=12)
	{
    Bitmap->Width=Rect.Width();
    Bitmap->Height=Rect.Height();
    TCanvas* Canvas=Bitmap->Canvas;
		TFont* Ft=Canvas->Font; Ft->Color=clWhite; Ft->Height=TextHeight; Ft->Name="Ariel";
    int fh=Canvas->TextHeight("0");

    double XX=Rect.right-Bdw*2, YY=Rect.bottom-Bdw-Bdwc-fh;
    Canvas->Brush->Color=clBlack; Canvas->Brush->Style=bsSolid; Canvas->FillRect(Rect);

    Canvas->Pen->Color=clSilver; Canvas->MoveTo(Bdw, Bdw+YY*0.5); Canvas->LineTo(Bdw+XX, Bdw+YY*0.5);

    double xrange=xen-xst, XXixrange=XX/xrange, xstXXixrange=xst*XXixrange,
					 xrangeiXX=xrange/XX,
           YYidbrange=YY/(dbst-dben);

    int *XRec=new int[XX];
    for (int m=0; m<M; m++)
    {
      bool visible=false;
      atom* Partialsm=Partials[m];
      Canvas->Brush->Color=RotateColors[m%10]; Canvas->Pen->Color=RotateColors[m%10];
      memset(XRec, 0, sizeof(int)*XX);

      for (int fr=frst; fr<fren; fr++)
      {
        try{
        if (Partialsm[fr].f>0 && Partialsm[fr].a>0)
        {
          double xx=Partialsm[fr].f;
          if (xx<xen) visible=true;
					if (xx<xst || xx>xen) continue;
					int X=floor(XXixrange*xx-xstXXixrange+0.5);
					if (X>=0 && X<XX && !XRec[X])
					{
						XRec[X]=true;
						double logaf;
							//calculate LogAF
							if (FScaleMode) xx=log(1+xx*Fs/700)/log(1+Fs/700);
							int k=floor(xx/F); double f_plus=xx/F-k;
							if (k<K+1) logaf=h[k]*(1-f_plus)+h[k+1]*f_plus;
							else logaf=h[K+1];

						double yy=20*log10e*(logaf+b[m]+logA0C[fr]);
            int Y=Bdw+floor((YY-YYidbrange*(yy-dben))+0.5);
            X+=Bdw;
            Canvas->MoveTo(X-1, Y); Canvas->LineTo(X+2, Y);
            Canvas->MoveTo(X, Y-1); Canvas->LineTo(X, Y+2);
          }
        }}catch(...)
        {int k=0;}
      }
      if (!visible) break;
    }
    delete[] XRec;

		AnsiString S=AnsiString().sprintf("%.2gkhz", Fs*xst/1000); Canvas->Brush->Style=bsClear; Canvas->TextOut(Bdw, Bdw+YY, S);
		S=AnsiString().sprintf("%.2gkhz", Fs*(xst+xen)/2000); int fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+(XX-fw)/2, Bdw+YY, S);
		S=AnsiString().sprintf("%.2gkhz", Fs*xen/1000); fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+XX-fw, Bdw+YY, S);
    S.sprintf("%gdB", dbst); fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+XX-fw, Bdw, S);
    S.sprintf("%gdB", (dbst+dben)/2); fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+XX-fw, Bdw+(YY-fh)/2, S);
    S.sprintf("%gdB", dben); fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+XX-fw, Bdw+YY-fh, S);

    Canvas->Brush->Color=clBlack; Canvas->Brush->Style=bsClear; Canvas->Pen->Color=clWhite;
    Canvas->Rectangle(Bdw, Bdw, Bdw+XX, Bdw+YY);
	}

	//*for source-filter model
	int DrawSF(Graphics::TBitmap* Bitmap, TRect Rect, double xst, double xen, double updb, double downdb, double f0, int M, double* b, double F, int K, double* h, int FScaleMode, double Fs=44100, int* SX=0, int* FX=0, int* FY=0, TColor* Colors=0, int TextHeight=12, int FShiftdB=0)
	{
		double FShift=FShiftdB/(20.0*log10e); //allowing the filter to be drawn with a vertical shift, and source an opposite shift
		Bitmap->Width=Rect.Width();
		Bitmap->Height=Rect.Height();
		TCanvas* Canvas=Bitmap->Canvas;
		double *Sx=new double[M];
		for (int i=0; i<M; i++) Sx[i]=20*log10e*(b[i]-FShift);

		int result=0;
		TFont* Ft=Canvas->Font; Ft->Color=clWhite; Ft->Height=TextHeight; Ft->Name="Ariel";
		int fh=Canvas->TextHeight("0");

		double XX=Rect.right-Bdw*2, YY=Rect.bottom-Bdw-Bdwc-fh;
		Canvas->Brush->Color=clBlack; Canvas->Brush->Style=bsSolid; Canvas->FillRect(Rect);
		Canvas->Pen->Color=clSilver; Canvas->MoveTo(Bdw, Bdw+YY*0.5); Canvas->LineTo(Bdw+XX, Bdw+YY*0.5);

		double idbrange=1.0/(updb-downdb), xrange=xen-xst, XXixrange=XX/xrange, xstXXixrange=xst*XXixrange,
			YY20log10eidbrange=YY*20*log10e*idbrange, BdwpYYpYYdbenidbrange=Bdw+YY+YY*downdb*idbrange,
			Bdw_xstXXixrange=Bdw-xstXXixrange, YYidbrange=YY*idbrange;

		double st=xst/F; int ist=ceil(st); if (ist<0) ist=0;
		Canvas->Pen->Color=clWhite;
		int X=Bdw+XX*(F*ist-xst)/xrange+0.5, Y=BdwpYYpYYdbenidbrange-YY20log10eidbrange*(h[ist]+FShift)+0.5;
		if (FX){FX[ist]=X, FY[ist]=Y; Canvas->MoveTo(X-1, Y); Canvas->LineTo(X+2, Y); Canvas->MoveTo(X, Y-1); Canvas->LineTo(X, Y+2);}
		Canvas->MoveTo(X, Y);
		//draw filter contour
		int LineCount=xrange/F;
		if (LineCount<XX)
		{
			for (int i=ist+1; i<K+2; i++)
			{
				double xi=F*i;
				if (FScaleMode)
				{
					double ximel=Amel*log(1+Fs/700)*xi;
					double xihz=700*(exp(ximel/Amel)-1);
					xi=xihz/Fs;
				}
				X=floor(Bdw_xstXXixrange+xi*XXixrange+0.5);
				Y=floor(BdwpYYpYYdbenidbrange-YY20log10eidbrange*(h[i]+FShift)+0.5);
				Canvas->LineTo(X, Y);
				if (FX){FX[i]=X, FY[i]=Y; Canvas->MoveTo(X-1, Y); Canvas->LineTo(X+2, Y); Canvas->MoveTo(X, Y-1); Canvas->LineTo(X, Y+2); Canvas->MoveTo(X, Y);}
			}
		}
		else
		{
			double di=xrange/(XX*F);
			for (int x=0; x<XX; x++)
			{
				int i;
				if (FScaleMode) i=log(1+(xst+xrange*x/XX)*Fs/700)/log(1+Fs/700)/F;
				else i=st+x*di;
				X=Bdw+x;
				Y=floor(BdwpYYpYYdbenidbrange-YY20log10eidbrange*(h[i]+FShift)+0.5);
				Canvas->LineTo(Bdw+x, Y);
				if (FX){FX[i]=X, FY[i]=Y; Canvas->MoveTo(X-1, Y); Canvas->LineTo(X+2, Y); Canvas->MoveTo(X, Y-1); Canvas->LineTo(X, Y+2); Canvas->MoveTo(X, Y);}
			}
		}
    //draw source lines
    for (int i=0; i<M; i++)
    {
			double xi=f0*(i+1);
      if (xi<xst || xi>xen) continue;
			int X=floor(Bdw_xstXXixrange+XXixrange*xi+0.5);
      int Y=floor(BdwpYYpYYdbenidbrange-YYidbrange*Sx[i]+0.5);
      if (Colors) Canvas->Pen->Color=Colors[i];
      else Canvas->Pen->Color=RotateColors[i%10];

      Canvas->MoveTo(X, Y); Canvas->LineTo(X, Bdw+YY);

			if (SX) SX[i]=X, result++;
    }
		AnsiString S=AnsiString().sprintf("%.2gkhz", Fs*xst/1000); Canvas->Brush->Style=bsClear; Canvas->TextOut(Bdw, Bdw+YY, S);
		S=AnsiString().sprintf("%.2gkhz", Fs*(xst+xen)/2000); int fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+(XX-fw)/2, Bdw+YY, S);
    S=AnsiString().sprintf("%.2gkhz", Fs*xen/1000); fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+XX-fw, Bdw+YY, S);
    S.sprintf("%gdB", updb); fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+XX-fw, Bdw, S);
    S.sprintf("%gdB", (updb+downdb)/2); fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+XX-fw, Bdw+(YY-fh)/2, S);
    S.sprintf("%gdB", downdb); fw=Canvas->TextWidth(S); Canvas->TextOut(Bdw+XX-fw, Bdw+YY-fh, S);

    Canvas->Brush->Color=clBlack; Canvas->Brush->Style=bsClear; Canvas->Pen->Color=clWhite;
    Canvas->Rectangle(Bdw, Bdw, Bdw+XX, Bdw+YY);
    delete[] Sx;
		return result;
	}   //*/

	void DrawF0(Graphics::TBitmap* Bitmap, TCanvas* Canvas, TRect Rect, double* lp, int P, double* F0, double f_c, double f_ex, double sps, double F0Overall, int Fr, atom** Partials, double L, bool Captions=true, TColor F0Color=clLime, bool peakmark=false, TColor peakmarkcolor=clYellow, double* frs=0, double* fs=0, int* BX1=0, int* BX2=0, int HLX1=0, int HLX2=0)
	{
    if (f_ex==0) f_ex=0.01;
		Canvas->Brush->Color=cl3DDkShadow; Canvas->FillRect(Rect);
		int Width=Rect.Width(), Height=Rect.Height();
		Bitmap->Width=Width; Bitmap->Height=Height;
		TPen* Pn=Canvas->Pen;
		TFont* F=Canvas->Font;

		if (HLX1<HLX2)
		{
			Canvas->Brush->Color=clDkGray; Canvas->FillRect(TRect(HLX1, 0, HLX2, Height));
			Canvas->Brush->Color=cl3DDkShadow;
		}

		int X, Y, Y1=0, Y2=Height, Y_;
    AnsiString S, as;

		Pn->Color=clBlack; Pn->Style=psDot; F->Color=clBlack; F->Height=12; F->Name="Ariel";
		int lp0=ceil(lp[0]), lpp=ceil(lp[P-1]);
    double f0max=-0.5, f0min=0.5; for (int fr=lp0; fr<lpp; fr++){if (f0max<F0[fr]) f0max=F0[fr]; if (f0min>F0[fr]) f0min=F0[fr];}
    Y1=Height*(0.5-(f0max-f_c)/f_ex); Canvas->MoveTo(0, Y1); Canvas->LineTo(Width, Y1);
		if (Captions){as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(f0max*sps)/C4)); S.sprintf("%.4ghz (%s)", f0max*sps, as.c_str()); Canvas->TextOut(1, Y1-Canvas->TextHeight(S)-1, S);}
    Y2=Height*(0.5-(f0min-f_c)/f_ex); Canvas->MoveTo(0, Y2); Canvas->LineTo(Width, Y2);
    if (Captions){as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(f0min*sps)/C4)); S.sprintf("%.4ghz (%s)", f0min*sps, as.c_str()); Canvas->TextOut(1, Y2+1, S);}

    if (peakmark)
    {
      int Y0=0.5*(Y1+Y2), Y0ex=0.5*(Y1-Y2);
			
      for (int pfr=0; pfr<P; pfr++)
      {
				int fr=floor(lp[pfr]); double frplus=lp[pfr]-fr;
				X=Width*(Partials[1][fr].t*(1-frplus)+Partials[1][fr+1].t*frplus)/L;
				Pn->Color=peakmarkcolor; Pn->Style=psDot; Canvas->MoveTo(X, Y1); Canvas->LineTo(X, Y2);
				if (BX1) BX1[pfr]=X-1;
				if (BX2) BX2[pfr]=X+1;
			}
    }

		if (Captions) {Y_=Height*(0.5-(F0Overall-f_c)/f_ex); Pn->Color=clWhite; Canvas->MoveTo(0, Y_); Canvas->LineTo(Width, Y_);}

    Pn->Style=psSolid; Pn->Color=F0Color;
    for (int fr=0; fr<Fr; fr++)
    {
      if (fr>0) Canvas->MoveTo(X, Y);
      X=Width*Partials[1][fr].t/L, Y=Height*(0.5-(F0[fr]-f_c)/f_ex);
      if (fr>0) Canvas->LineTo(X, Y);
    }

    if (frs)
    {
      Pn->Color=clRed;
      for (int pfr=0; pfr<P-1; pfr++)
      {
        int ifrs=frs[pfr]; double rfrs=frs[pfr]-ifrs;
        X=(Width*Partials[1][ifrs].t*(1-rfrs)+Width*Partials[1][ifrs+1].t*rfrs)/L, Y=Height*(0.5-(fs[pfr]-f_c)/f_ex);
        Canvas->MoveTo(X-4, Y); Canvas->LineTo(X+5, Y);
				Canvas->MoveTo(X, Y-4); Canvas->LineTo(X, Y+5);
      }
		}

    if (Captions)
    {
      as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(F0Overall*sps)/C4)); S.sprintf("%.4ghz (%s)", F0Overall*sps, as.c_str());
      F->Color=clWhite; F->Style=F->Style<<fsBold; Canvas->TextOut(Width-Canvas->TextWidth(S), Y_+1, S); F->Style=F->Style>>fsBold;
      F->Color=clSilver; S="F0"; Canvas->TextOut(1, Height-Canvas->TextHeight(S), S);
      double f=(f_c+f_ex/2)*sps; as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(f)/C4)); S.sprintf("%.4ghz (%s)", f, as.c_str()); Canvas->TextOut(Width-Canvas->TextWidth(S), 0, S);
			f=(f_c-f_ex/2)*sps; as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(f)/C4)); S.sprintf("%.4ghz (%s)", f, as.c_str()); Canvas->TextOut(Width-Canvas->TextWidth(S), Height-Canvas->TextHeight(S), S);
    }
  }

	void DrawF0C(Graphics::TBitmap* Bitmap, TCanvas* Canvas, TRect Rect, double* F0C, double& F0Cmax, double& F0Cmin, double f_c, double f_ex, double sps, double F0Overall, int Fr, atom** Partials, double L, bool clearbackground=true, int frcount=0, double* frs=0, double* fs=0, int* CX1=0, int* CX2=0, int* CY1=0, int* CY2=0)
  {
    if (f_ex==0) f_ex=f_c*0.005;
    //Draw the F0 carrier
    Canvas->Brush->Color=cl3DDkShadow; if (clearbackground) Canvas->FillRect(Rect);
    int Width=Rect.Width(), Height=Rect.Height();
    Bitmap->Width=Width; Bitmap->Height=Height;
    TPen* P=Canvas->Pen;
    TFont* F=Canvas->Font;

    P->Color=clBlack;
		F0Cmax=f_c-f_ex, F0Cmin=f_c+f_ex;
    for (int fr=0; fr<Fr; fr++){if (F0Cmax<F0C[fr]) F0Cmax=F0C[fr]; if (F0Cmin>F0C[fr]) F0Cmin=F0C[fr];}
    P->Style=psDot;
    int X, Y=Height*(0.5-(F0Cmax-f_c)/f_ex);
    Canvas->MoveTo(0, Y); Canvas->LineTo(Width, Y);
    AnsiString as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(F0Cmax*sps)/C4)); AnsiString S; S.sprintf("%.4ghz (%s)", F0Cmax*sps, as.c_str());
    F->Color=clBlack; F->Height=12; F->Name="Ariel"; Canvas->TextOut(1, Y-Canvas->TextHeight(S)-1, S);
    Y=Height*(0.5-(F0Cmin-f_c)/f_ex);
    Canvas->MoveTo(0, Y); Canvas->LineTo(Width, Y);
    as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(F0Cmin*sps)/C4)); S.sprintf("%.4ghz (%s)", F0Cmin*sps, as.c_str());
    Canvas->TextOut(1, Y+1, S);
    P->Color=clWhite;
    as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(F0Overall*sps)/C4)); S.sprintf("%.4ghz (%s)", F0Overall*sps, as.c_str());
    int Y_=Height*(0.5-(F0Overall-f_c)/f_ex);
    Canvas->MoveTo(0, Y_); Canvas->LineTo(Width, Y_);
    F->Color=clWhite; F->Style=F->Style<<fsBold;
    if (Y_+1+Canvas->TextHeight(S)<Y) Canvas->TextOut(Width-Canvas->TextWidth(S), Y_+1, S);
    else Canvas->TextOut(Width-Canvas->TextWidth(S), Y+1, S);
    F->Style=F->Style>>fsBold;

    P->Style=psSolid; P->Color=clLime;
    for (int fr=0; fr<Fr; fr++)
    {
      if (fr>0) Canvas->MoveTo(X, Y);
			X=Width*Partials[1][fr].t/L, Y=Height*(0.5-(F0C[fr]-f_c)/f_ex);
      if (fr>0) Canvas->LineTo(X, Y);
    }

    if (frcount)
    {
			P->Color=clRed;
      for (int p=0; p<frcount; p++)
      {
        int ifrs=frs[p]; double rfrs=frs[p]-ifrs;
        X=(Width*Partials[1][ifrs].t*(1-rfrs)+Width*Partials[1][ifrs+1].t*rfrs)/L, Y=Height*(0.5-(fs[p]-f_c)/f_ex);
        Canvas->MoveTo(X-4, Y); Canvas->LineTo(X+5, Y);
        Canvas->MoveTo(X, Y-4); Canvas->LineTo(X, Y+5);
        if (CX1) CX1[p]=X-4;
        if (CX2) CX2[p]=X+4;
        if (CY1) CY1[p]=Y-4;
        if (CY2) CY2[p]=Y+4;
      }
    }

    double f=(f_c+f_ex/2)*sps;
    F->Color=clSilver;
    as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(f)/C4)); S.sprintf("%.4ghz (%s)", f, as.c_str());
    Canvas->TextOut(Width-Canvas->TextWidth(S), 0, S);
		f=(f_c-f_ex/2)*sps; as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(f)/C4)); S.sprintf("%.4ghz (%s)", f, as.c_str());
    Canvas->TextOut(Width-Canvas->TextWidth(S), Height-Canvas->TextHeight(S), S);
    F->Color=clSilver; S="F0 carrier";
    Canvas->TextOut(1, Height-Canvas->TextHeight(S), S);
  }

  //y=ax+b
	void linefit(int n, double* x, double* y, double* indi, double& a, double& b)
  {
    double sx=0, sy=0, sxy=0, sxx=0, m=0;
    for (int i=0; i<n; i++) if (indi[i]>0) sx+=x[i], sy+=y[i], sxx+=x[i]*x[i], sxy+=x[i]*y[i], m++;
    if (m<=0) {a=-1; b=-1; return;}
    double id=1.0/(m*sxx-sx*sx);
    a=(m*sxy-sx*sy)*id;
    b=(sxx*sy-sx*sxy)*id;
  }



void TSFDemoForm::UpdateDisplay(bool f0, bool f0c, bool f0d, bool sf)
{
	int Fr=HS->Fr;
  atom** Partials=HS->Partials;
	//find the highest point of the 2nd partial
	double f2min=1, f2max=0, fst=0, fen=fmaxEdit->Text.ToDouble();
	for (int fr=0; fr<Fr; fr++)
  {
		if (f2max<Partials[1][fr].f) f2max=Partials[1][fr].f;
    if (f2min>Partials[1][fr].f) f2min=Partials[1][fr].f;
  }
  f_c=(f2min+f2max)/4; f_ex=f2max-f2min;

	if (f0)
	{
		int HLX1=0, HLX2=0;
		if (CurrentP>=0) HLX1=BX1[CurrentP]+1, HLX2=BX1[CurrentP+1]+1;
		DrawF0(Image0->Picture->Bitmap, Image0->Canvas, Image0->ClientRect, SF.lp, SF.P, SF.F0, f_c, f_ex, WaveView1->SamplesPerSec, SF.F0Overall, Fr, Partials, WaveView2->Length, true, clLime, true, clYellow, 0, cyclefs, BX1, BX2, HLX1, HLX2);
	}
	if (f0c)
	{
		memset(CX1, 0xFF, sizeof(int)*128); memset(CX2, 0xFF, sizeof(int)*128); memset(CY1, 0xFF, sizeof(int)*128); memset(CY1, 0xFF, sizeof(int)*128);
		DrawF0C(Image1->Picture->Bitmap, Image1->Canvas, Image1->ClientRect, SF.F0C, SF.F0Cmax, SF.F0Cmin, f_c, f_ex, WaveView1->SamplesPerSec, SF.F0Overall, Fr, Partials, WaveView2->Length, true, cyclefrcount, cyclefrs, cyclefs, CX1, CX2, CY1, CY2);
  }
  if (f0d)
  {
    //Draw the F0 modulator
    DrawModulator();
  }

	if (sf)
	{
		int lp0=ceil(SF.lp[0]), lpp=ceil(SF.lp[SF.P-1]);
		double updb=updbEdit->Text.ToDouble(), downdb=downdbEdit->Text.ToDouble(), dbshift=SFPicShift->Text.ToDouble();
		dbrange=updb-downdb;
		if (SFCheck->Checked)
			DrawSFr(AImage1->Picture->Bitmap, AImage1->ClientRect, fst, fen, updb, downdb, SF.M, lp0, lpp, Partials, useA0?SF.logA0:SF.logA0C, SF.avgb, SF.F, SF.K, SF.avgh, SF.FScaleMode, SF.Fs);
		else
			DrawAF(AImage1->Picture->Bitmap, AImage1->ClientRect, fst, fen, updb, downdb, SF.M, lp0, lpp, Partials, MACheck->Checked, false, WaveView1->SamplesPerSec);
		dbrange=updb-downdb;
		SXc=DrawSF(AImage3->Picture->Bitmap, AImage3->ClientRect, fst, fen, updb+dbshift, downdb+dbshift, SF.F0Overall, SF.M, SF.avgb, SF.F, SF.K, SF.avgh, SF.FScaleMode, 44100, SX, FX, FY, 0, 12, FShiftdBEdit->Text.ToInt());
	}

  ForceUpdate=false;
}

void __fastcall TSFDemoForm::WaveView1OpMode(TObject* Sender, TShiftState Shift, int& OpMode)
{
  if (Shift.Contains(ssShift)) OpMode=0; //drag mode
  else OpMode=1; //select mode
}

int __fastcall TSFDemoForm::WaveView1CustomInfo(TObject* Sender)
{
  TWaveView* WV=(TWaveView*)Sender;
  TStringList* List=new TStringList;
	double fs=WV->StartPos*1.0/WV->SamplesPerSec, fe=WV->EndPos*1.0/WV->SamplesPerSec;
	List->Add(AnsiString().sprintf(" Time (%.4gs): from %.4gs to %.4gs. ", fe-fs, fs, fe));
	List->Add(AnsiString().sprintf(" Frequency: from %.1fhz to %.1fhz. ", WV->StartDigiFreq*WV->SamplesPerSec, WV->EndDigiFreq*WV->SamplesPerSec));
	return int(List);
}

void __fastcall TSFDemoForm::Image1MouseMove(TObject *Sender,
			TShiftState Shift, int X, int Y)
{
	if (!HS) return;
	TImage* Image=(TImage*)Sender;
	Image->Parent->SetFocus();

	double t, sps=WaveView1->SamplesPerSec, f; AnsiString S, as; TFont* F=Image->Canvas->Font;

	if (Sender==Image0)
	{
		if (Shift.Contains(ssLeft))
		{
			if (CurrentB<0 || CurrentB>=SF.P) return;
			if (X<0 || X>=Image->Width || CurrentB>0 && X<=BX2[CurrentB-1] || CurrentB<SF.P-1 && X>=BX1[CurrentB+1]){}
			else
			{
				double dlp=1.0*(X-StartDragX)*SF.L/Image->Width;
				SF.lp[CurrentB]+=dlp;
				StartDragX=X;
			}
		}
		else
		{
			CurrentB=-1;
			for (int i=0; i<SF.P; i++) if (X>=BX1[i] && X<=BX2[i]) {CurrentB=i; break;}
		}
	}
	else if (Sender==Image1)
	{
		if (Shift.Contains(ssLeft))
		{
			if (CurrentC<0 || CurrentC>=cyclefrcount) return;
			double df=-(Y-StartDrag)*f_ex/Image->Height;
			cyclefs[CurrentC]+=df;
			StartDrag=Y;
			Smooth_Interpolate(SF.F0C, SF.L, cyclefrcount+1, cyclefs, cyclefrs);
			SynthesizeSF(HS, &SF, WaveView1->SamplesPerSec);
			UpdateDisplay(true, true, true, false);

				t=HS->Partials[0][0].t+SF.offst*cyclefrs[CurrentC];
				f=cyclefs[CurrentC]*sps;
				as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(f)/C4));
				S.sprintf("%.3gs, %.4ghz (%s)     ", t/sps, f, as.c_str());
				F->Color=clRed;
				Image->Canvas->TextOut(1, Image->Canvas->TextHeight("0"), S);
		}
		else
		{
			CurrentC=-1;
			for (int i=0; i<cyclefrcount; i++)
				if (X>=CX1[i] && X<=CX2[i] && Y>=CY1[i] && Y<=CY2[i]) {CurrentC=i; break;}
		}
			if (CurrentC>=0 && CurrentC<cyclefrcount)
			{
				t=HS->Partials[0][0].t+SF.offst*cyclefrs[CurrentC];
				f=cyclefs[CurrentC]*sps;
				as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(f)/C4));
				S.sprintf("%.3gs, %.4ghz (%s)     ", t/sps, f, as.c_str());
				F->Color=clRed;
				Image->Canvas->TextOut(1, Image->Canvas->TextHeight("0"), S);
			}
			else
			{
				S.sprintf("                                        ", t/sps, f, as.c_str());
				Image->Canvas->TextOut(1, Image->Canvas->TextHeight("0"), S);
			}
	}

	t=WaveView2->Length*1.0*(X+0.5)/Image->Width;
	f=(f_c-(Y-0.5*Image->Height)*f_ex/Image->Height)*sps;
	as=SemitoneToPitch(12*Log2(WV2_LOG_FREQ(f)/C4));
	S.sprintf("%.3gs, %.4ghz (%s)       ", t/sps, f, as.c_str());
	F->Color=clAqua; F->Height=12; F->Name="Ariel";
	Image->Canvas->TextOut(1, 0, S);


	if (Sender==Image0)
	{
		int Wid=HS->Partials[0][1].s, hWid=Wid/2, Offst=HS->Partials[0][1].t-HS->Partials[0][0].t;
		CurrentFr=(t-hWid)/Offst+0.5;
		if (!Shift.Contains(ssLeft)) CurrentP=-1; while (CurrentP<SF.P-1 && CurrentFr>SF.lp[CurrentP+1]) CurrentP++; if (CurrentP>=SF.P-1) CurrentP=-1;

		if (CurrentP>=0)
    {
      PageControl1->ActivePage=AmpCycleSheet;
      SFCheck->Parent=Panel14;
      MACheck->Parent=Panel14;
    }
		else
    {
      PageControl1->ActivePage=AmpOverSheet;
      SFCheck->Parent=Panel12;
      MACheck->Parent=Panel12;
    }
		UpdateDisplay(true, 0, CurrentP>=0, 0);

			if (CurrentB>=0 && CurrentB<SF.P)
			{
				t=HS->Partials[0][0].t+SF.offst*SF.lp[CurrentB];
				S.sprintf("%.3gs", t/sps);
				F->Color=clYellow;
				Image->Canvas->TextOut(1, Image->Canvas->TextHeight("0"), S);
			}
			else
			{
				S.sprintf("          ");
				Image->Canvas->TextOut(1, Image->Canvas->TextHeight("0"), S);
			}
	}
}
//---------------------------------------------------------------------------

void TSFDemoForm::DrawModulator()
{
	if (PageControl1->ActivePage==AmpCycleSheet && CurrentP>=0 && CurrentP<SF.P-1)
	{
			double fst=0, fen=fmaxEdit->Text.ToDouble();
			int lp0=ceil(SF.lp[0]);
			double updb=updbEdit->Text.ToDouble(), downdb=downdbEdit->Text.ToDouble(), dbshift=SFPicShift->Text.ToDouble();
			dbrange=updb-downdb;
			if (SFCheck->Checked)
				DrawSFr(AImage2->Picture->Bitmap, AImage2->ClientRect, fst, fen, updb, downdb, SF.M, ceil(SF.lp[CurrentP]), ceil(SF.lp[CurrentP+1]), HS->Partials, useA0?SF.logA0:SF.logA0C, SF.b[CurrentFr-lp0], SF.F, SF.K, SF.h[CurrentFr-lp0], SF.FScaleMode, SF.Fs);
			else
				DrawAF(AImage2->Picture->Bitmap, AImage2->ClientRect, fst, fen, updb, downdb, SF.M, ceil(SF.lp[CurrentP]), ceil(SF.lp[CurrentP+1]), HS->Partials, MACheck->Checked, false, WaveView1->SamplesPerSec);
			dbrange=updb-downdb;
			SXcyclec=DrawSF(AImage4->Picture->Bitmap, AImage4->ClientRect, fst, fen, updb+dbshift, downdb+dbshift, SF.F0[CurrentFr], SF.M, SF.b[CurrentFr-lp0], SF.F, SF.K, SF.h[CurrentFr-lp0], SF.FScaleMode, SF.Fs, SXcycle, 0, 0, 0, 12, FShiftdBEdit->Text.ToInt());
	}
}


void __fastcall TSFDemoForm::SFCheckClick(TObject *Sender)
{
  if (Sender==SFCheck)
  {
    MACheck->Visible=!SFCheck->Checked;
  }
  UpdateDisplay();
}
//---------------------------------------------------------------------------


void __fastcall TSFDemoForm::ExternalInput()
{
  AnsiString FileName=ExtractFilePath(Application->ExeName)+"tsfin";
	FILE* file;
	if (file=fopen(FileName.c_str(), "rb"))
  {
    SF.LoadFromFileHandle(file);
    fclose(file);
    DeleteFile(FileName);
		SynthesizeSF(HS, &SF, WaveView1->SamplesPerSec);

		SaveSF();
    Synthesize();
    UpdateDisplay();
  }
}

//---------------------------------------------------------------------------
void __fastcall TSFDemoForm::SaveSF()
{
	FILE* file;
	if (file=fopen((ExtractFilePath(Application->ExeName)+"tsfout").c_str(), "wb"))
  {
    SF.SaveToFileHandle(file);
    fclose(file);
  }
}

void __fastcall TSFDemoForm::SaveSF0()
{
  SF0.Duplicate(SF);
}

void __fastcall TSFDemoForm::SaveWave()
{
  THS* HS=Form1->HS; int dst, den;
  double* xrec=SynthesisHSp(HS, dst, den);
  delete[] Wave1; Wave1=new __int16[WaveView2->Length];
  DoubleToInt(Wave1, 2, xrec, den-dst);
  free8(xrec);
}
//---------------------------------------------------------------------------
void __fastcall TSFDemoForm::MouseWheelHandler(TMessage& Msg)
{
  TWMMouseWheel* WMsg=(TWMMouseWheel*)&Msg;

  bool Handled=false;
  TMouseWheelEvent Event=NULL;
  TControl* Sender;
  TShiftState Shift; memcpy(&Shift, &WMsg->Keys, 1);

  if (ActiveControl==Image1->Parent) Event=Image1MouseWheel, Sender=Image1;
	else if (ActiveControl==AImage3->Parent) Event=AImage3MouseWheel, Sender=AImage3;
	else if (ActiveControl==AImage1->Parent) Event=AImage3MouseWheel, Sender=AImage1;

  if (Event) Event(Sender, Shift, WMsg->WheelDelta, Sender->ScreenToClient(TPoint(WMsg->Pos.x, WMsg->Pos.y)), Handled);
  WMsg->Result=Handled;

  if (!Handled) TCustomForm::MouseWheelHandler(Msg);
}

//---------------------------------------------------------------------------
void __fastcall TSFDemoForm::Image1MouseWheel(TObject* Sender, TShiftState Shift, int WheelDelta, const TPoint &MousePos, bool &Handled)
{
  double cent;
  if (Shift.Contains(ssShift)) cent=1;
  else if (Shift.Contains(ssCtrl)) cent=100;
  else cent=10;

  if (WheelDelta<0) cent=-cent;
  double rate=pow(2, cent/1200);

  for (int l=0; l<SF.L; l++) SF.F0C[l]*=rate, SF.F0D[l]*=rate; SF.F0Overall*=rate;
  for (int fr=0; fr<cyclefrcount; fr++) cyclefs[fr]*=rate;

	SynthesizeSF(HS, &SF, WaveView1->SamplesPerSec);
	SaveSF();
	S_U();
  Handled=true;
}

void __fastcall TSFDemoForm::AImage3MouseWheel(TObject* Sender, TShiftState Shift, int WheelDelta, const TPoint &MousePos, bool &Handled)
{
  double ddb;
  if (Shift.Contains(ssShift)) ddb=0.2;
  else if (Shift.Contains(ssCtrl)) ddb=5;
  else ddb=1;
  double db=20*log10e*SF.avgb[CurrentPartial];
  if (WheelDelta<0) ddb=-ddb;

  db+=ddb;
  AnsiString S; S.sprintf("Partial %d: %.2fdB", CurrentPartial+1, db);
	Label7->Caption=S;
	  ddb=ddb/(20*log10e);

	int lp0=ceil(SF.lp[0]), lpp=ceil(SF.lp[SF.P-1]);
	for (int l=0; l<lpp-lp0; l++) SF.b[l][CurrentPartial]+=ddb;
	SF.avgb[CurrentPartial]+=ddb;

	SynthesizeSF(HS, &SF, WaveView1->SamplesPerSec);
  SaveSF();
  S_U(true);
  Handled=true;
}


//---------------------------------------------------------------------------
void __fastcall TSFDemoForm::AImage3MouseMove(TObject *Sender,
			TShiftState Shift, int X, int Y)
{
	if (!HS) return;
	TImage* Image=(TImage*)Sender;
	int lp0=ceil(SF.lp[0]), lpp=ceil(SF.lp[SF.P-1]);

	if (Shift.Contains(ssLeft))
	{
		if (CurrentFB<0 && CurrentPartial<0) return;
		int fh=Image->Canvas->TextHeight("0");
		int YY=Image->Height-Bdw-Bdwc-fh;
		double dbperpixel=dbrange/YY;
		int dY=Y-StartDrag;
    if (dY==0) return;
		StartDrag=Y;
		double ddb=-dY*dbperpixel;

		if (Sender==AImage3 && CurrentFB>=0)
		{
			double db=20*log10e*SF.avgh[CurrentPartial];
			db+=ddb;
			AnsiString S; S.sprintf("Filter control %d: %.2fdB", CurrentFB, db);
			Label7->Caption=S;
			ddb=ddb/(20*log10e);

			for (int l=0; l<lpp-lp0; l++) SF.h[l][CurrentFB]+=ddb;
			SF.avgh[CurrentFB]+=ddb;

			if (Shift.Contains(ssShift))
				for (int k=CurrentFB+1; k<=SF.K+1; k++) {for (int l=0; l<lpp-lp0; l++) SF.h[l][k]+=ddb; SF.avgh[k]+=ddb;}
			if (Shift.Contains(ssCtrl)) {S_F_b(SF, HS->Partials); SaveSF0();}
		}
		else if (Sender==AImage1 || Sender==AImage3)
		{
			double db=20*log10e*SF.avgb[CurrentPartial];
			db+=ddb;
			AnsiString S; S.sprintf("Partial %d: %.2fdB", CurrentPartial+1, db);
			Label7->Caption=S;
			ddb=ddb/(20*log10e);

			for (int l=0; l<lpp-lp0; l++) SF.b[l][CurrentPartial]+=ddb;
			SF.avgb[CurrentPartial]+=ddb;

      int dm=0; for (int i=0x31; i<0x3a; i++) if (GetKeyState(i)<0) dm+=(i-0x30);
      if (dm==0 && Shift.Contains(ssShift)) dm=CurrentPartial+1;
      if (dm>0) for (int m=CurrentPartial+dm; m<SF.M; m+=dm){for (int l=0; l<lpp-lp0; l++) SF.b[l][m]+=ddb; SF.avgb[m]+=ddb;}
		}
		else if (Sender==AImage2 || Sender==AImage4)
		{
			if (CurrentFr>=lp0 && CurrentFr<lpp)
			{
				double db=20*log10e*SF.b[CurrentFr-lp0][CurrentPartial];
				db+=ddb;
				AnsiString S; S.sprintf("Partial %d: %.2fdB", CurrentPartial+1, db);
				Label8->Caption=S;
				ddb=ddb/(20*log10e);

				SF.b[CurrentFr-lp0][CurrentPartial]+=ddb;
			}
		}

		SynthesizeSF(HS, &SF, WaveView1->SamplesPerSec);
		UpdateDisplay();
	}
	else
	{
		Image->Parent->SetFocus();

		int lSXc, *lSX;
		double* LogAS;
		if (Sender==AImage1 || Sender==AImage3) lSXc=SXc, lSX=SX, LogAS=SF.avgb, X=Bdw+(X-Bdw)*AImage3->Width/((TImage*)Sender)->Width;
		else if (Sender==AImage2 || Sender==AImage4)
		{
			if (CurrentFr<lp0 || CurrentFr>=lpp) return;
			lSXc=SXcyclec, lSX=SXcycle, LogAS=SF.b[CurrentFr-lp0], X=Bdw+(X-Bdw)*AImage4->Width/((TImage*)Sender)->Width;
		}

    if (lSXc<=0) return;

		CurrentPartial=0; while (CurrentPartial<lSXc && lSX[CurrentPartial]<X) CurrentPartial++;
		if (CurrentPartial>=lSXc) CurrentPartial=lSXc-1;
		if (CurrentPartial>=1 && X-lSX[CurrentPartial-1]<lSX[CurrentPartial]-X) CurrentPartial--;
		CurrentFB=-1; int lfb=0; while (lfb<=SF.K && X>FX[lfb]) lfb++;
		if (Sender==AImage3)
		{
			if (X>=FX[lfb]-ftrr && X<=FX[lfb]+ftrr && Y>=FY[lfb]-ftrr && Y<=FY[lfb]+ftrr) CurrentFB=lfb;
			if (X>=FX[lfb-1]-ftrr && X<=FX[lfb-1]+ftrr && Y>=FY[lfb-1]-ftrr && Y<=FY[lfb-1]+ftrr)
			{
				if (CurrentFB==-1) CurrentFB=lfb-1;
				else if (abs(X-FX[lfb-1])<abs(X-FX[lfb])) CurrentFB=lfb-1;
			}
		}
		
		TLabel* Label; AnsiString S;
		if (Sender==AImage1 || Sender==AImage3) Label=Label7;
		else if (Sender==AImage2 || Sender==AImage4) Label=Label8;

		if (CurrentFB>=0)
		{
			S.sprintf("Filter control %d: %.2fdB", CurrentFB, 20*log10e*SF.avgh[CurrentFB]);
			Label->Font->Color=clWhite;
		}
		else
		{
			double dB=20*log10e*LogAS[CurrentPartial];
			S.sprintf("Partial %d: %.2fdB", CurrentPartial+1, dB);
			Label->Font->Color=RotateColors[CurrentPartial%10];
		}
		Label->Caption=S;
	}
}
//---------------------------------------------------------------------------


void __fastcall TSFDemoForm::AImage1MouseDown(TObject *Sender,
			TMouseButton Button, TShiftState Shift, int X, int Y)
{
	FirstStartDrag=StartDrag=Y;
	FirstStartDragX=StartDragX=X;
	if (Sender==Image0 && Shift.Contains(ssShift))
	{
		if (CurrentB>=0) //delete current peakmark
		{
			for (int p=CurrentB; p<SF.P-1; p++) SF.lp[p]=SF.lp[p+1];
			SF.P--;
			CurrentB=-1;
			UpdateDisplay(1, 0, 0, 0);
			analyze2=true;
		}
		else
		{
			double newlpb=1.0*X/Image0->Width*SF.L, *newlp=new double[SF.P+3];
			int newb=0; while (newb<SF.P && X>BX1[newb]) newb++;

			if (newb>0) memcpy(newlp, SF.lp, sizeof(double)*newb);
			if (newb<SF.P) memcpy(&newlp[newb+1], &SF.lp[newb], sizeof(double)*(SF.P-newb));
			newlp[newb]=newlpb;
			CurrentB=newb;
			SF.P++;
			delete[] SF.lp; SF.lp=newlp;
			UpdateDisplay(1, 0, 0, 0);
			analyze2=true;
		}
	}
  if ((Sender==AImage2 || Sender==AImage4) && Shift.Contains(ssRight))
    PageControl1->ActivePage=AmpOverSheet;
}
//---------------------------------------------------------------------------

void __fastcall TSFDemoForm::AImage1MouseUp(TObject *Sender,
			TMouseButton Button, TShiftState Shift, int X, int Y)
{
	if (Sender==Image0 && (CurrentB>=0 && X!=FirstStartDragX || analyze2))
	{
    //do analysis using the new segmentation
		double sps=WaveView1->SamplesPerSec;
		double h=WaveView1->SpecOffst;
		int FSMode=MethodListBox->ItemIndex;
		double FSF=FEdit->Text.ToDouble();
		int FSFScale=FScaleCombo->ItemIndex;
		if (FSFScale==0) FSF/=sps;
		else FSF/=Amel*log(1+sps/700);
		double FStheta=ThetaEdit->Text.ToDouble();
		AnalyzeSF_2(*HS, SF, cyclefrs, cyclefs, sps, &cyclefrcount, FSMode, FSF, FSFScale, FStheta);
    SaveSF0();
		SaveSF();
		S_U(true);
		UpdateDisplay();
		analyze2=false;
	}
	else if (Y!=FirstStartDrag)
	{
		SaveSF();
		S_U();
		UpdateDisplay();
	}
}
//---------------------------------------------------------------------------


//---------------------------------------------------------------------------

void __fastcall TSFDemoForm::PageControl1Change(TObject *Sender)
{
  if (PageControl1->ActivePage==AmpOverSheet)
  {
    SFCheck->Parent=Panel12;
    MACheck->Parent=Panel12;
  }
  else
  {
    SFCheck->Parent=Panel14;
    MACheck->Parent=Panel14;
  }
  UpdateDisplay();
}
//---------------------------------------------------------------------------

void __fastcall TSFDemoForm::SaveButtonClick(TObject *Sender)
{
  int k=1; AnsiString FileName="1.sf";
  while (TmplListBox->Items->IndexOf(FileName)>=0) {k++; FileName=AnsiString(k)+".sf";}
  if (GetKeyState(VK_SHIFT)<0)
  {
    AnsiString oldpath=SaveDialog1->InitialDir;
    SaveDialog1->InitialDir=ExtractFileDir(Application->ExeName);
    if (SaveDialog1->Execute()) FileName=SaveDialog1->FileName;
    SaveDialog1->InitialDir=oldpath;
  }
  else FileName=ExtractFilePath(Application->ExeName)+FileName;
  SF.SaveToFile(FileName.c_str());
  TmplListBox->Items->Add(ExtractFileName(FileName));
}
//---------------------------------------------------------------------------


  //calculates increment dh on sf1.avgh so that it becomes sf2.avgh
  // dh should have the same dimenstions as sf1.avgh
  void* CompareFilters(double* dh, TSF& sf1, TSF& sf2)
  {
    int K1=sf1.K, FScale1=sf1.FScaleMode;
    double F1=sf1.F, Fs=sf1.Fs;
    for (int k=0; k<K1+2; k++)
    {
      double h1=sf1.avgh[k], f=F1*k;
      if (FScale1)
      {
  			double fmel=Amel*log(1+Fs/700)*f;
	  		double fhz=700*(exp(fmel/Amel)-1);
		  	f=fhz/Fs;
      }
      double h2=sf2.LogAF(f);
      dh[k]=h2-h1;
    }
  }
  double* CompareFilters(TSF& sf1, TSF& sf2)
  {
    double* dh=new double[sf1.K+2];
    CompareFilters(dh, sf1, sf2);
    return dh;
  }

void __fastcall TSFDemoForm::UseButtonClick(TObject *Sender)
{
  if (TmplListBox->ItemIndex<0) return;
	AnsiString FileName=ExtractFilePath(Application->ExeName)+TmplListBox->Items->Strings[TmplListBox->ItemIndex];
  if (!FileExists(FileName)) return;
  double dbamp=TrackBar1->Position/10.0;
  double dbgain=TrackBar2->Position/5.0-10;

  TSF* sf=new TSF;
	FILE* file;
	if (file=fopen(FileName.c_str(), "rb"))
	{
    sf->LoadFromFileHandle(file);
    fclose(file);
  }
  if (dbgain!=0) sf->ShiftFilterByDB(dbgain);
  int L=ceil(SF.lp[SF.P-1])-ceil(SF.lp[0]), K=SF.K, M=SF.M;
  double avgbend=SF0.avgb[sf->M-1];
  switch (ComboBox1->ItemIndex)
  {
		case 0: //use source
    {
      for (int l=0; l<L; l++) memcpy(SF.h[l], SF0.h[l], sizeof(double)*(K+2));
      memcpy (SF.avgh, SF0.avgh, sizeof(double)*(K+2));
      for (int m=0; m<M; m++)
      {
        double sfb=(m<sf->M)?sf->avgb[m]:(sf->avgb[sf->M-1]-avgbend+SF0.avgb[m]);
        double db=(sfb-SF0.avgb[m])*dbamp;
        for (int l=0; l<L; l++) SF.b[l][m]=SF0.b[l][m]+db;
        SF.avgb[m]=SF0.avgb[m]+db;
      }
      break;
    }
    case 1: //use filter
    {
      for (int l=0; l<L; l++) memcpy(SF.b[l], SF0.b[l], sizeof(double)*M);
      memcpy(SF.avgb, SF0.avgb, sizeof(double)*M); 
      double* dh=CompareFilters(SF0, *sf);

      for (int k=0; k<K+2; k++){double dhk=dh[k]*dbamp;
        for (int l=0; l<L; l++) SF.h[l][k]=SF0.h[l][k]+dhk; SF.avgh[k]=SF0.avgh[k]+dhk;}
      delete[] dh;
      break;
    }
    case 2: //source and filter
    {
      for (int m=0; m<M; m++)
      {
        double sfb=(m<sf->M)?sf->avgb[m]:(sf->avgb[sf->M-1]-avgbend+SF0.avgb[m]);
        double db=(sfb-SF0.avgb[m])*dbamp;
        for (int l=0; l<L; l++) SF.b[l][m]=SF0.b[l][m]+db; SF.avgb[m]=SF0.avgb[m]+db;
      }
      double* dh=CompareFilters(SF0, *sf);

      for (int k=0; k<K+2; k++){double dhk=dh[k]*dbamp;
        for (int l=0; l<L; l++) SF.h[l][k]=SF0.h[l][k]+dhk; SF.avgh[k]=SF0.avgh[k]+dhk;}
      delete[] dh;
      break;
    }
    case 3: //towards
    {
      for (int m=0; m<M; m++)
      {
        double sfb=(m<sf->M)?sf->avgb[m]:(sf->avgb[sf->M-1]-avgbend+SF0.avgb[m]);
        double db=(sfb-SF0.avgb[m])*dbamp, avgb=0;
        for (int l=0; l<L; l++){SF.b[l][m]=SF0.b[l][m]+db*l/(L-1); avgb+=SF.b[l][m];}
        SF.avgb[m]=avgb/L;
      }
      double *dh=CompareFilters(SF0, *sf);

      for (int k=0; k<K+2; k++){double dhk=dh[k]*dbamp, avgh=0;
        for (int l=0; l<L; l++){SF.h[l][k]=SF0.h[l][k]+dhk*l/(L-1); avgh+=SF.h[l][k];} SF.avgh[k]=avgh/L;}
      delete[] dh;
      break;
    }
    case 4: //towards and back
    {
      for (int m=0; m<M; m++)
      {
        double sfb=(m<sf->M)?sf->avgb[m]:(sf->avgb[sf->M-1]-avgbend+SF0.avgb[m]);
        double db=(sfb-SF0.avgb[m])*dbamp, avgb=0;
        for (int l=0; l<L; l++){SF.b[l][m]=SF0.b[l][m]+db*(1-fabs(l*2.0/(L-1)-1)); avgb+=SF.b[l][m];}
        SF.avgb[m]=avgb/L;
      }
      double *dh=CompareFilters(SF0, *sf);

      for (int k=0; k<K+2; k++){double dhk=dh[k]*dbamp, avgh=0;
        for (int l=0; l<L; l++){SF.h[l][k]=SF0.h[l][k]+dhk*(1-fabs(l*2.0/(L-1)-1)); avgh+=SF.h[l][k];} SF.avgh[k]=avgh/L;}
      delete[] dh;
      break;
    }
      break;
  }

  delete sf;

	SynthesizeSF(HS, &SF, WaveView1->SamplesPerSec);
  SaveSF();
  S_U(true);
  UpdateDisplay();
}
//---------------------------------------------------------------------------
//S_U implements delayed updating of lengthy synthesizing to speed up operation
//  in case the HS is updated frequently, s.a. during mouse wheeling.
//This allows a later synthesis to start before an earlier synthesis finishes,
//  at which the latter is terminated and its incomplete result is discarded.
void __fastcall TSFDemoForm::S_U(bool sf)
{
  if (SUThread)
  {
    SUThread->Terminate();
  }
	SUThread=new TSUThread(true);
	SUThread->sf=sf;
	SUThread->OnTerminate=SUTerminate;
	SUThread->HS=new THS(HS);
	ThreadList[pThread++]=SUThread;
	if (pThread==ThreadCaps/2 && ThreadList[pThread])
		for (int i=pThread; i<ThreadCaps; i++) {delete ThreadList[i]; ThreadList[i]=0;}
	else if (pThread==ThreadCaps)
	{
		for (int i=0; i<ThreadCaps/2; i++) {delete ThreadList[i]; ThreadList[i]=0;}
		pThread=0;
	}
	SUThread->Resume();
}

void __fastcall TSFDemoForm::SUTerminate(TObject* Sender)
{
	TSUThread* CurrentThread=(TSUThread*)Sender;
  if (CurrentThread==SUThread)
  {
    double* xrec=CurrentThread->xrec;
		if (xrec)
		{
			int dst=CurrentThread->dst, den=CurrentThread->den, bps=WaveAudio2->BitsPerSample;
			if (dst<0) dst=0;
			double max=1<<(bps-1); for (int i=dst; i<den; i++) if (xrec[i]>max) max=xrec[i]; else if (xrec[i]<-max) max=-xrec[i];
			if (max>1<<(bps-1))
			{
				max=(1<<(bps-1))/max;
				for (int i=dst; i<den; i++) xrec[i]*=max;
			}
			void* data=new char[(den-dst)*bps];
			DoubleToInt(data, bps>>3, xrec, den-dst);
      WaveAudio2->Clear(NULL);
      WaveAudio2->WriteSamples(data, den-dst);
      delete[] data;
      UpdateDisplay();
    }
    SUThread=0;
  }
  else
  {
    UpdateDisplay(true, true, true, CurrentThread->sf);
  }
  free8(CurrentThread->xrec); CurrentThread->xrec=0;
}
//---------------------------------------------------------------------------

void __fastcall TSFDemoForm::Panel4Resize(TObject *Sender)
{
	ForceUpdate=true;
}
//---------------------------------------------------------------------------


void __fastcall TSFDemoForm::Panel1Resize(TObject *Sender)
{
  Panel2->Height=(Panel1->Height-Splitter1->Height)/2;  
}
//---------------------------------------------------------------------------



void __fastcall TSFDemoForm::AmpCycleSheetResize(TObject *Sender)
{
	UpdateDisplay();
}
//---------------------------------------------------------------------------


void __fastcall TSFDemoForm::MethodListBoxClick(TObject *Sender)
{
	Reset();	
}
//---------------------------------------------------------------------------

void __fastcall TSFDemoForm::AmpOverSheetResize(TObject *Sender)
{
	UpdateDisplay();
}
//---------------------------------------------------------------------------

void __fastcall TSFDemoForm::Splitter7Moved(TObject *Sender)
{
	if (Sender==Splitter7) AImage4->Parent->Parent->Height=AImage3->Parent->Parent->Height;
	else if (Sender==Splitter8) AImage3->Parent->Parent->Height=AImage4->Parent->Parent->Height;
	UpdateDisplay();
}
//---------------------------------------------------------------------------

void __fastcall TSFDemoForm::AImage3DblClick(TObject *Sender)
{
	Graphics::TBitmap* bmp=new Graphics::TBitmap;
	bmp->Width=800; bmp->Height=600; TRect Rect=TRect(0, 0, 800, 600);


			double fst=0, fen=fmaxEdit->Text.ToDouble();

			double updb=updbEdit->Text.ToDouble(), downdb=downdbEdit->Text.ToDouble(), dbshift=SFPicShift->Text.ToDouble();
			dbrange=updb-downdb;

			if (Sender==AImage3) DrawSF(bmp, Rect, fst, fen, updb+dbshift, downdb=dbshift, SF.F0Overall, SF.M, SF.avgb, SF.F, SF.K, SF.avgh, SF.FScaleMode, 44100, SX, 0, 0, 0, 36, FShiftdBEdit->Text.ToInt());
			else if (Sender==AImage1)
			{
				int lp0=ceil(SF.lp[0]), lpp=ceil(SF.lp[SF.P-1]);
				if (!SFCheck->Checked)
					DrawAF(bmp, Rect, fst, fen, updb, downdb, SF.M, lp0, lpp, HS->Partials, MACheck->Checked, false, WaveView1->SamplesPerSec, 36);
				else
					DrawSFr(bmp, Rect, fst, fen, updb, downdb, SF.M, lp0, lpp, HS->Partials, useA0?SF.logA0:SF.logA0C, SF.avgb, SF.F, SF.K, SF.avgh, SF.FScaleMode, SF.Fs, 36);
			}
  Form1->SaveDialog1->FileName="waveview.waveview1.bmp";
  Form1->SaveDialog1->FilterIndex=2;
  if (Form1->SaveDialog1->Execute()) bmp->SaveToFile(Form1->SaveDialog1->FileName);
	delete bmp;
}
//---------------------------------------------------------------------------

void __fastcall TSFDemoForm::Panel7Resize(TObject *Sender)
{
	UpdateDisplay(1, 0, 0, 0);
}
//---------------------------------------------------------------------------

void __fastcall TSFDemoForm::Panel8Resize(TObject *Sender)
{
	UpdateDisplay(0, 1, 0, 0);
}
//---------------------------------------------------------------------------


void __fastcall TSFDemoForm::CheckBox1Click(TObject *Sender)
{
	WaveView2->LoopPlay=CheckBox1->Checked;		
}
//---------------------------------------------------------------------------

void __fastcall TSFDemoForm::fmaxEditKeyPress(TObject *Sender, char &Key)
{
  if (Key==VK_RETURN) UpdateDisplay(0, 0, 0, 1);
}
//---------------------------------------------------------------------------

void __fastcall TSFDemoForm::Button1Click(TObject *Sender)
{
  int N=WaveView1->Length, Channel=Form1->HS->Channel;

  if (GetKeyState(VK_SHIFT)<0)
  {
    __int16* d2=WaveView2->Data16[0];
    Form1->PostWaveViewData(d2, Channel, StartPos, StartPos+N, Form1->FadeInCheck->Checked, Form1->FadeInCombo->Text.ToInt());
  }
  else
  {
    __int16 *data=new __int16[N], *data1=WaveView1->Data16[0], *data2=WaveView2->Data16[0];
    for (int n=0; n<N; n++) data[n]=data1[n]-datain[n]+data2[n];
    Form1->PostWaveViewData(data, Channel, StartPos, StartPos+N, Form1->FadeInCheck->Checked, Form1->FadeInCombo->Text.ToInt());
    delete[] data;
  }
}
//---------------------------------------------------------------------------