content: /**
   * Copyright 2007 DFKI GmbH.
   * All Rights Reserved.  Use is subject to license terms.
   * <p/>
   * Permission is hereby granted, free of charge, to use and distribute
   * this software and its documentation without restriction, including
   * without limitation the rights to use, copy, modify, merge, publish,
   * distribute, sublicense, and/or sell copies of this work, and to
   * permit persons to whom this work is furnished to do so, subject to
   * the following conditions:
   * <p/>
   * 1. The code must retain the above copyright notice, this list of
   * conditions and the following disclaimer.
   * 2. Any modifications must be clearly marked as such.
   * 3. Original authors' names are not deleted.
   * 4. The authors' names are not used to endorse or promote products
   * derived from this software without specific prior written
   * permission.
   * <p/>
   * DFKI GMBH AND THE CONTRIBUTORS TO THIS WORK DISCLAIM ALL WARRANTIES WITH
   * REGARD TO THIS SOFTWARE, INCLUDING ALL IMPLIED WARRANTIES OF
   * MERCHANTABILITY AND FITNESS, IN NO EVENT SHALL DFKI GMBH NOR THE
   * CONTRIBUTORS BE LIABLE FOR ANY SPECIAL, INDIRECT OR CONSEQUENTIAL
   * DAMAGES OR ANY DAMAGES WHATSOEVER RESULTING FROM LOSS OF USE, DATA OR
   * PROFITS, WHETHER IN AN ACTION OF CONTRACT, NEGLIGENCE OR OTHER TORTIOUS
   * ACTION, ARISING OUT OF OR IN CONNECTION WITH THE USE OR PERFORMANCE OF
   * THIS SOFTWARE.
   */
  
  package marytts.signalproc.sinusoidal.hntm.synthesis;
  
  import java.util.Arrays;
  
  import marytts.signalproc.sinusoidal.hntm.analysis.HntmAnalyzerParams;
  import marytts.signalproc.sinusoidal.hntm.analysis.HntmPlusTransientsSpeechSignal;
  import marytts.signalproc.window.HammingWindow;
  import marytts.signalproc.window.Window;
  import marytts.util.signal.SignalProcUtils;
  
  /**
   * Syntehsizer for the transient part waveform segments.
   *
   * @author oytun.turk
   */
  public class TransientPartSynthesizer {
  
      public static double[] synthesize(HntmPlusTransientsSpeechSignal hnmSignal, HntmAnalyzerParams analysisParams) {
          int outputLen = SignalProcUtils.time2sample(hnmSignal.originalDurationInSeconds, hnmSignal.samplingRateInHz);
  
          double[] transientPart = new double[outputLen];
          Arrays.fill(transientPart, 0.0);
  
          if (hnmSignal.transients != null) {
              int i, j;
              int startInd, endInd;
              int windowLeftEndInd, windowRightStartInd;
              int ws = SignalProcUtils.time2sample(2 * analysisParams.overlapBetweenTransientAndNontransientSectionsInSeconds, hnmSignal.samplingRateInHz);
              if (ws % 2 == 0)
                  ws++;
              Window win = new HammingWindow(ws);
              win.normalizePeakValue(1.0f);
              int winMidInd = (ws - 1) / 2;
              for (i = 0; i < hnmSignal.transients.segments.length; i++) {
                  if (hnmSignal.transients.segments[i] != null && hnmSignal.transients.segments[i].waveform != null && hnmSignal.transients.segments[i].waveform.length > 0 && hnmSignal.transients.segments[i].startTime >= 0.0f) {
                      startInd = Math.min(SignalProcUtils.time2sample(hnmSignal.transients.segments[i].startTime, hnmSignal.samplingRateInHz), outputLen - 1);
                      windowLeftEndInd = Math.min(startInd + winMidInd, outputLen - 1);
                      endInd = Math.min(SignalProcUtils.time2sample(hnmSignal.transients.segments[i].startTime, hnmSignal.samplingRateInHz) + hnmSignal.transients.segments[i].waveform.length - 1, outputLen - 1);
                      windowRightStartInd = endInd - winMidInd;
  
                      for (j = startInd; j <= windowLeftEndInd; j++)
                          transientPart[j] = hnmSignal.transients.segments[i].waveform[j - startInd] * win.value(j - startInd);
                      for (j = windowLeftEndInd + 1; j < windowRightStartInd; j++)
                          transientPart[j] = hnmSignal.transients.segments[i].waveform[j - startInd];
                      for (j = windowRightStartInd; j <= endInd; j++)
                          transientPart[j] = hnmSignal.transients.segments[i].waveform[j - startInd] * win.value((j - endInd) + ws - 1);
                  }
              }
          }
  
          return transientPart;
      }
  }
  