svoid audioRestorationDemo_v1() {
  S text =
    //"please tell me something nice";
    "you can almost understand this now";
  bool sine = false;
  
  double seconds = 1.0;
  new ContinuousOscillators_TestAudioRestoration test;
  test.showImage(true);
  var co = test.co;
  S waveName;
  test.separateFrequencies(true);
  test.separateRealAndImag(false);
  test.cumulative(true);
  co.renderUsingSine = false;
  
  if (sine) {
    test.setAudio(seconds, new SineSoundSource(440.0, test.inputSampleRate));
    //if (true) ret;
    waveName = "sine";
  } else {
    test.loadMP3(new GSay(text).audioFile());
    waveName = text;
  }
  
  //test.co.setFrequencies(pianoFrequencies88());
  test.co.setFrequencies(map toFrequency(geometricIterator(55.0, test.inputSampleRate, 2.0)));
  test.silence = 0.01;
  //test.minFreq = 440;
  //test.maxFreq = 440*4;
  test.wavOut = javaxDataDir("Restored Audio/" + waveName + ".restored." + (test.co.renderUsingSine ? "sine" : "square") + ".wav");
  test.run();
  
  /*
  double interval = co.sampleRate/880;
  print(+interval);
  int n = 3;
  L<Double> period = map(x -> x/32767.0, asList(takeFirst(iceil(interval*n), test.resampled)));
  L<Average> partSum = AutoCreateList(Average);
  for i over period: {
    var sample = period.get(i);
    int part = ifloor(i/interval*4);
    var sum = partSum.get(part);
    sum.add(sample);
    printVars(+i, +part, sample := formatDouble2(sample), sum := formatDouble2(sum!));
  }
  pnl(partSum);
  
  for (int i = 0; i+3 < l(partSum); i += 4) {
    double[] vals = new[2];
    for phase to 2: {
      double a = partSum.get(i+((4-phase)&3))!;
      double b = partSum.get(i+((5-phase)&3))!;
      double c = partSum.get(i+((6-phase)&3))!;
      double d = partSum.get(i+((7-phase)&3))!;
    
      printVars(+a, +b, +c, +d);
      double hi = avg(a, b), lo = avg(c, d);
      double val = hi-lo;
      vals[phase] = val;
      print("i=" + i + ", phase " + phase + ": " + formatDouble2(hi) + " - " + formatDouble2(lo) + " = " + formatDouble2(val));
    }
    print("c=" + renderComplexWithAngle(complex(vals)));
  }
  */
  
  // test linear interpolation
  /*for (t : countIterator(0.0, 2.0, 0.1))
    print(t + "=" + co.audio.sampleSum(0, 0, t));*/
}