diff src/rubberband-1.8.1/main/main.cpp @ 10:37bf6b4a2645

Add FFTW3

--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/src/rubberband-1.8.1/main/main.cpp	Wed Mar 20 15:35:50 2013 +0000
@@ -0,0 +1,648 @@
+/* -*- c-basic-offset: 4 indent-tabs-mode: nil -*-  vi:set ts=8 sts=4 sw=4: */
+
+/*
+    Rubber Band Library
+    An audio time-stretching and pitch-shifting library.
+    Copyright 2007-2012 Particular Programs Ltd.
+
+    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.  See the file
+    COPYING included with this distribution for more information.
+
+    Alternatively, if you have a valid commercial licence for the
+    Rubber Band Library obtained by agreement with the copyright
+    holders, you may redistribute and/or modify it under the terms
+    described in that licence.
+
+    If you wish to distribute code using the Rubber Band Library
+    under terms other than those of the GNU General Public License,
+    you must obtain a valid commercial licence before doing so.
+*/
+
+#include "rubberband/RubberBandStretcher.h"
+
+#include <iostream>
+#include <sndfile.h>
+#include <cmath>
+#include <time.h>
+#include <cstdlib>
+#include <cstring>
+#include <string>
+
+#include <fstream>
+
+#include "system/sysutils.h"
+
+#ifdef __MSVC__
+#include "getopt/getopt.h"
+#else
+#include <getopt.h>
+#include <unistd.h>
+#include <sys/time.h>
+#endif
+
+#include "base/Profiler.h"
+
+using namespace std;
+using namespace RubberBand;
+
+#ifdef _WIN32
+using RubberBand::gettimeofday;
+#endif
+
+#ifdef __MSVC__
+using RubberBand::usleep;
+#endif
+
+double tempo_convert(const char *str)
+{
+    char *d = strchr((char *)str, ':');
+
+    if (!d || !*d) {
+        double m = atof(str);
+        if (m != 0.0) return 1.0 / m;
+        else return 1.0;
+    }
+
+    char *a = strdup(str);
+    char *b = strdup(d+1);
+    a[d-str] = '\0';
+    double m = atof(a);
+    double n = atof(b);
+    free(a);
+    free(b);
+    if (n != 0.0 && m != 0.0) return m / n;
+    else return 1.0;
+}
+
+int main(int argc, char **argv)
+{
+    int c;
+
+    double ratio = 1.0;
+    double duration = 0.0;
+    double pitchshift = 0.0;
+    double frequencyshift = 1.0;
+    int debug = 0;
+    bool realtime = false;
+    bool precise = true;
+    int threading = 0;
+    bool lamination = true;
+    bool longwin = false;
+    bool shortwin = false;
+    bool smoothing = false;
+    bool hqpitch = false;
+    bool formant = false;
+    bool together = false;
+    bool crispchanged = false;
+    int crispness = -1;
+    bool help = false;
+    bool version = false;
+    bool quiet = false;
+
+    bool haveRatio = false;
+
+    std::string mapfile;
+
+    enum {
+        NoTransients,
+        BandLimitedTransients,
+        Transients
+    } transients = Transients;
+
+    enum {
+        CompoundDetector,
+        PercussiveDetector,
+        SoftDetector
+    } detector = CompoundDetector;
+
+    while (1) {
+        int optionIndex = 0;
+
+        static struct option longOpts[] = {
+            { "help",          0, 0, 'h' },
+            { "version",       0, 0, 'V' },
+            { "time",          1, 0, 't' },
+            { "tempo",         1, 0, 'T' },
+            { "duration",      1, 0, 'D' },
+            { "pitch",         1, 0, 'p' },
+            { "frequency",     1, 0, 'f' },
+            { "crisp",         1, 0, 'c' },
+            { "crispness",     1, 0, 'c' },
+            { "debug",         1, 0, 'd' },
+            { "realtime",      0, 0, 'R' },
+            { "loose",         0, 0, 'L' },
+            { "precise",       0, 0, 'P' },
+            { "formant",       0, 0, 'F' },
+            { "no-threads",    0, 0, '0' },
+            { "no-transients", 0, 0, '1' },
+            { "no-lamination", 0, 0, '2' },
+            { "centre-focus",  0, 0, '7' },
+            { "window-long",   0, 0, '3' },
+            { "window-short",  0, 0, '4' },
+            { "bl-transients", 0, 0, '8' },
+            { "detector-perc", 0, 0, '5' },
+            { "detector-soft", 0, 0, '6' },
+            { "smoothing",     0, 0, '9' },
+            { "pitch-hq",      0, 0, '%' },
+            { "threads",       0, 0, '@' },
+            { "quiet",         0, 0, 'q' },
+            { "timemap",       1, 0, 'M' },
+            { 0, 0, 0, 0 }
+        };
+
+        c = getopt_long(argc, argv,
+                        "t:p:d:RLPFc:f:T:D:qhVM:",
+                        longOpts, &optionIndex);
+        if (c == -1) break;
+
+        switch (c) {
+        case 'h': help = true; break;
+        case 'V': version = true; break;
+        case 't': ratio *= atof(optarg); haveRatio = true; break;
+        case 'T': ratio *= tempo_convert(optarg); haveRatio = true; break;
+        case 'D': duration = atof(optarg); haveRatio = true; break;
+        case 'p': pitchshift = atof(optarg); haveRatio = true; break;
+        case 'f': frequencyshift = atof(optarg); haveRatio = true; break;
+        case 'd': debug = atoi(optarg); break;
+        case 'R': realtime = true; break;
+        case 'L': precise = false; break;
+        case 'P': precise = true; break;
+	case 'F': formant = true; break;
+        case '0': threading = 1; break;
+        case '@': threading = 2; break;
+        case '1': transients = NoTransients; crispchanged = true; break;
+        case '2': lamination = false; crispchanged = true; break;
+        case '3': longwin = true; crispchanged = true; break;
+        case '4': shortwin = true; crispchanged = true; break;
+        case '5': detector = PercussiveDetector; crispchanged = true; break;
+        case '6': detector = SoftDetector; crispchanged = true; break;
+        case '7': together = true; break;
+        case '8': transients = BandLimitedTransients; crispchanged = true; break;
+        case '9': smoothing = true; crispchanged = true; break;
+        case '%': hqpitch = true; break;
+        case 'c': crispness = atoi(optarg); break;
+        case 'q': quiet = true; break;
+        case 'M': mapfile = optarg; break;
+        default:  help = true; break;
+        }
+    }
+
+    if (version) {
+        cerr << RUBBERBAND_VERSION << endl;
+        return 0;
+    }
+
+    if (help || !haveRatio || optind + 2 != argc) {
+        cerr << endl;
+	cerr << "Rubber Band" << endl;
+        cerr << "An audio time-stretching and pitch-shifting library and utility program." << endl;
+	cerr << "Copyright 2007-2012 Particular Programs Ltd." << endl;
+        cerr << endl;
+	cerr << "   Usage: " << argv[0] << " [options] <infile.wav> <outfile.wav>" << endl;
+        cerr << endl;
+        cerr << "You must specify at least one of the following time and pitch ratio options." << endl;
+        cerr << endl;
+        cerr << "  -t<X>, --time <X>       Stretch to X times original duration, or" << endl;
+        cerr << "  -T<X>, --tempo <X>      Change tempo by multiple X (same as --time 1/X), or" << endl;
+        cerr << "  -T<X>, --tempo <X>:<Y>  Change tempo from X to Y (same as --time X/Y), or" << endl;
+        cerr << "  -D<X>, --duration <X>   Stretch or squash to make output file X seconds long" << endl;
+        cerr << endl;
+        cerr << "  -p<X>, --pitch <X>      Raise pitch by X semitones, or" << endl;
+        cerr << "  -f<X>, --frequency <X>  Change frequency by multiple X" << endl;
+        cerr << endl;
+        cerr << "  -M<F>, --timemap <F>    Use file F as the source for key frame map" << endl;
+        cerr << endl;
+        cerr << "A map file consists of a series of lines each having two numbers separated" << endl;
+        cerr << "by a single space.  These are source and target sample frame numbers for fixed" << endl;
+        cerr << "time points within the audio data, defining a varying stretch factor through" << endl;
+        cerr << "the audio.  You must specify an overall stretch factor using e.g. -t as well." << endl;
+        cerr << endl;
+        cerr << "The following options provide a simple way to adjust the sound.  See below" << endl;
+        cerr << "for more details." << endl;
+        cerr << endl;
+        cerr << "  -c<N>, --crisp <N>      Crispness (N = 0,1,2,3,4,5,6); default 5 (see below)" << endl;
+	cerr << "  -F,    --formant        Enable formant preservation when pitch shifting" << endl;
+        cerr << endl;
+        cerr << "The remaining options fine-tune the processing mode and stretch algorithm." << endl;
+        cerr << "These are mostly included for test purposes; the default settings and standard" << endl;
+        cerr << "crispness parameter are intended to provide the best sounding set of options" << endl;
+        cerr << "for most situations.  The default is to use none of these options." << endl;
+        cerr << endl;
+        cerr << "  -L,    --loose          Relax timing in hope of better transient preservation" << endl;
+        cerr << "  -P,    --precise        Ignored: The opposite of -L, this is default from 1.6" << endl;
+        cerr << "  -R,    --realtime       Select realtime mode (implies --no-threads)" << endl;
+        cerr << "         --no-threads     No extra threads regardless of CPU and channel count" << endl;
+        cerr << "         --threads        Assume multi-CPU even if only one CPU is identified" << endl;
+        cerr << "         --no-transients  Disable phase resynchronisation at transients" << endl;
+        cerr << "         --bl-transients  Band-limit phase resync to extreme frequencies" << endl;
+        cerr << "         --no-lamination  Disable phase lamination" << endl;
+        cerr << "         --window-long    Use longer processing window (actual size may vary)" << endl;
+        cerr << "         --window-short   Use shorter processing window" << endl;
+        cerr << "         --smoothing      Apply window presum and time-domain smoothing" << endl;
+        cerr << "         --detector-perc  Use percussive transient detector (as in pre-1.5)" << endl;
+        cerr << "         --detector-soft  Use soft transient detector" << endl;
+        cerr << "         --pitch-hq       In RT mode, use a slower, higher quality pitch shift" << endl;
+        cerr << "         --centre-focus   Preserve focus of centre material in stereo" << endl;
+        cerr << "                          (at a cost in width and individual channel quality)" << endl;
+        cerr << endl;
+        cerr << "  -d<N>, --debug <N>      Select debug level (N = 0,1,2,3); default 0, full 3" << endl;
+        cerr << "                          (N.B. debug level 3 includes audible ticks in output)" << endl;
+        cerr << "  -q,    --quiet          Suppress progress output" << endl;
+        cerr << endl;
+        cerr << "  -V,    --version        Show version number and exit" << endl;
+        cerr << "  -h,    --help           Show this help" << endl;
+        cerr << endl;
+        cerr << "\"Crispness\" levels:" << endl;
+        cerr << "  -c 0   equivalent to --no-transients --no-lamination --window-long" << endl;
+        cerr << "  -c 1   equivalent to --detector-soft --no-lamination --window-long (for piano)" << endl;
+        cerr << "  -c 2   equivalent to --no-transients --no-lamination" << endl;
+        cerr << "  -c 3   equivalent to --no-transients" << endl;
+        cerr << "  -c 4   equivalent to --bl-transients" << endl;
+        cerr << "  -c 5   default processing options" << endl;
+        cerr << "  -c 6   equivalent to --no-lamination --window-short (may be good for drums)" << endl;
+        cerr << endl;
+	return 2;
+    }
+
+    if (crispness >= 0 && crispchanged) {
+        cerr << "WARNING: Both crispness option and transients, lamination or window options" << endl;
+        cerr << "         provided -- crispness will override these other options" << endl;
+    }
+
+    switch (crispness) {
+    case -1: crispness = 5; break;
+    case 0: detector = CompoundDetector; transients = NoTransients; lamination = false; longwin = true; shortwin = false; break;
+    case 1: detector = SoftDetector; transients = Transients; lamination = false; longwin = true; shortwin = false; break;
+    case 2: detector = CompoundDetector; transients = NoTransients; lamination = false; longwin = false; shortwin = false; break;
+    case 3: detector = CompoundDetector; transients = NoTransients; lamination = true; longwin = false; shortwin = false; break;
+    case 4: detector = CompoundDetector; transients = BandLimitedTransients; lamination = true; longwin = false; shortwin = false; break;
+    case 5: detector = CompoundDetector; transients = Transients; lamination = true; longwin = false; shortwin = false; break;
+    case 6: detector = CompoundDetector; transients = Transients; lamination = false; longwin = false; shortwin = true; break;
+    };
+
+    if (!quiet) {
+        cerr << "Using crispness level: " << crispness << " (";
+        switch (crispness) {
+        case 0: cerr << "Mushy"; break;
+        case 1: cerr << "Piano"; break;
+        case 2: cerr << "Smooth"; break;
+        case 3: cerr << "Balanced multitimbral mixture"; break;
+        case 4: cerr << "Unpitched percussion with stable notes"; break;
+        case 5: cerr << "Crisp monophonic instrumental"; break;
+        case 6: cerr << "Unpitched solo percussion"; break;
+        }
+        cerr << ")" << endl;
+    }
+
+    std::map<size_t, size_t> mapping;
+    
+    if (mapfile != "") {
+        std::ifstream ifile(mapfile.c_str());
+        if (!ifile.is_open()) {
+            cerr << "ERROR: Failed to open time map file \"" << mapfile << "\""
+                 << endl;
+            return 1;
+        }
+        std::string line;
+        int lineno = 0;
+        while (!ifile.eof()) {
+            std::getline(ifile, line);
+            while (line.length() > 0 && line[0] == ' ') line = line.substr(1);
+            if (line == "") {
+                ++lineno;
+                continue;
+            }
+            std::string::size_type i = line.find_first_of(" ");
+            if (i == std::string::npos) {
+                cerr << "ERROR: Time map file \"" << mapfile
+                     << "\" is malformed at line " << lineno << endl;
+                return 1;
+            }
+            size_t source = atoi(line.substr(0, i).c_str());
+            while (i < line.length() && line[i] == ' ') ++i;
+            size_t target = atoi(line.substr(i).c_str());
+            mapping[source] = target;
+            if (debug > 0) {
+                cerr << "adding mapping from " << source << " to " << target << endl;
+            }
+            ++lineno;
+        }
+        ifile.close();
+
+        if (!quiet) {
+            cerr << "Read " << mapping.size() << " line(s) from map file" << endl;
+        }
+    }
+
+    char *fileName = strdup(argv[optind++]);
+    char *fileNameOut = strdup(argv[optind++]);
+
+    SNDFILE *sndfile;
+    SNDFILE *sndfileOut;
+    SF_INFO sfinfo;
+    SF_INFO sfinfoOut;
+    memset(&sfinfo, 0, sizeof(SF_INFO));
+
+    sndfile = sf_open(fileName, SFM_READ, &sfinfo);
+    if (!sndfile) {
+	cerr << "ERROR: Failed to open input file \"" << fileName << "\": "
+	     << sf_strerror(sndfile) << endl;
+	return 1;
+    }
+
+    if (duration != 0.0) {
+        if (sfinfo.frames == 0 || sfinfo.samplerate == 0) {
+            cerr << "ERROR: File lacks frame count or sample rate in header, cannot use --duration" << endl;
+            return 1;
+        }
+        double induration = double(sfinfo.frames) / double(sfinfo.samplerate);
+        if (induration != 0.0) ratio = duration / induration;
+    }
+
+    sfinfoOut.channels = sfinfo.channels;
+    sfinfoOut.format = sfinfo.format;
+    sfinfoOut.frames = int(sfinfo.frames * ratio + 0.1);
+    sfinfoOut.samplerate = sfinfo.samplerate;
+    sfinfoOut.sections = sfinfo.sections;
+    sfinfoOut.seekable = sfinfo.seekable;
+
+    sndfileOut = sf_open(fileNameOut, SFM_WRITE, &sfinfoOut) ;
+    if (!sndfileOut) {
+	cerr << "ERROR: Failed to open output file \"" << fileNameOut << "\" for writing: "
+	     << sf_strerror(sndfileOut) << endl;
+	return 1;
+    }
+    
+    int ibs = 1024;
+    size_t channels = sfinfo.channels;
+
+    RubberBandStretcher::Options options = 0;
+    if (realtime)    options |= RubberBandStretcher::OptionProcessRealTime;
+    if (precise)     options |= RubberBandStretcher::OptionStretchPrecise;
+    if (!lamination) options |= RubberBandStretcher::OptionPhaseIndependent;
+    if (longwin)     options |= RubberBandStretcher::OptionWindowLong;
+    if (shortwin)    options |= RubberBandStretcher::OptionWindowShort;
+    if (smoothing)   options |= RubberBandStretcher::OptionSmoothingOn;
+    if (formant)     options |= RubberBandStretcher::OptionFormantPreserved;
+    if (hqpitch)     options |= RubberBandStretcher::OptionPitchHighQuality;
+    if (together)    options |= RubberBandStretcher::OptionChannelsTogether;
+
+    switch (threading) {
+    case 0:
+        options |= RubberBandStretcher::OptionThreadingAuto;
+        break;
+    case 1:
+        options |= RubberBandStretcher::OptionThreadingNever;
+        break;
+    case 2:
+        options |= RubberBandStretcher::OptionThreadingAlways;
+        break;
+    }
+
+    switch (transients) {
+    case NoTransients:
+        options |= RubberBandStretcher::OptionTransientsSmooth;
+        break;
+    case BandLimitedTransients:
+        options |= RubberBandStretcher::OptionTransientsMixed;
+        break;
+    case Transients:
+        options |= RubberBandStretcher::OptionTransientsCrisp;
+        break;
+    }
+
+    switch (detector) {
+    case CompoundDetector:
+        options |= RubberBandStretcher::OptionDetectorCompound;
+        break;
+    case PercussiveDetector:
+        options |= RubberBandStretcher::OptionDetectorPercussive;
+        break;
+    case SoftDetector:
+        options |= RubberBandStretcher::OptionDetectorSoft;
+        break;
+    }
+
+    if (pitchshift != 0.0) {
+        frequencyshift *= pow(2.0, pitchshift / 12);
+    }
+
+    cerr << "Using time ratio " << ratio;
+    cerr << " and frequency ratio " << frequencyshift << endl;
+
+#ifdef _WIN32
+    RubberBand::
+#endif
+    timeval tv;
+    (void)gettimeofday(&tv, 0);
+
+    RubberBandStretcher::setDefaultDebugLevel(debug);
+
+    RubberBandStretcher ts(sfinfo.samplerate, channels, options,
+                           ratio, frequencyshift);
+
+    ts.setExpectedInputDuration(sfinfo.frames);
+
+    float *fbuf = new float[channels * ibs];
+    float **ibuf = new float *[channels];
+    for (size_t i = 0; i < channels; ++i) ibuf[i] = new float[ibs];
+
+    int frame = 0;
+    int percent = 0;
+
+    sf_seek(sndfile, 0, SEEK_SET);
+
+    if (!realtime) {
+
+        if (!quiet) {
+            cerr << "Pass 1: Studying..." << endl;
+        }
+
+        while (frame < sfinfo.frames) {
+
+            int count = -1;
+
+            if ((count = sf_readf_float(sndfile, fbuf, ibs)) <= 0) break;
+        
+            for (size_t c = 0; c < channels; ++c) {
+                for (int i = 0; i < count; ++i) {
+                    float value = fbuf[i * channels + c];
+                    ibuf[c][i] = value;
+                }
+            }
+
+            bool final = (frame + ibs >= sfinfo.frames);
+
+            ts.study(ibuf, count, final);
+
+            int p = int((double(frame) * 100.0) / sfinfo.frames);
+            if (p > percent || frame == 0) {
+                percent = p;
+                if (!quiet) {
+                    cerr << "\r" << percent << "% ";
+                }
+            }
+
+            frame += ibs;
+        }
+
+        if (!quiet) {
+            cerr << "\rCalculating profile..." << endl;
+        }
+
+        sf_seek(sndfile, 0, SEEK_SET);
+    }
+
+    frame = 0;
+    percent = 0;
+
+    if (!mapping.empty()) {
+        ts.setKeyFrameMap(mapping);
+    }
+    
+    size_t countIn = 0, countOut = 0;
+
+    while (frame < sfinfo.frames) {
+
+        int count = -1;
+
+	if ((count = sf_readf_float(sndfile, fbuf, ibs)) < 0) break;
+        
+        countIn += count;
+
+        for (size_t c = 0; c < channels; ++c) {
+            for (int i = 0; i < count; ++i) {
+                float value = fbuf[i * channels + c];
+                ibuf[c][i] = value;
+            }
+        }
+
+        bool final = (frame + ibs >= sfinfo.frames);
+
+        if (debug > 2) {
+            cerr << "count = " << count << ", ibs = " << ibs << ", frame = " << frame << ", frames = " << sfinfo.frames << ", final = " << final << endl;
+        }
+
+        ts.process(ibuf, count, final);
+
+        int avail = ts.available();
+        if (debug > 1) cerr << "available = " << avail << endl;
+
+        if (avail > 0) {
+            float **obf = new float *[channels];
+            for (size_t i = 0; i < channels; ++i) {
+                obf[i] = new float[avail];
+            }
+            ts.retrieve(obf, avail);
+            countOut += avail;
+            float *fobf = new float[channels * avail];
+            for (size_t c = 0; c < channels; ++c) {
+                for (int i = 0; i < avail; ++i) {
+                    float value = obf[c][i];
+                    if (value > 1.f) value = 1.f;
+                    if (value < -1.f) value = -1.f;
+                    fobf[i * channels + c] = value;
+                }
+            }
+//            cout << "fobf mean: ";
+//    double d = 0;
+//    for (int i = 0; i < avail; ++i) {
+//        d += fobf[i];
+//    }
+//    d /= avail;
+//    cout << d << endl;
+            sf_writef_float(sndfileOut, fobf, avail);
+            delete[] fobf;
+            for (size_t i = 0; i < channels; ++i) {
+                delete[] obf[i];
+            }
+            delete[] obf;
+        }
+
+        if (frame == 0 && !realtime && !quiet) {
+            cerr << "Pass 2: Processing..." << endl;
+        }
+
+	int p = int((double(frame) * 100.0) / sfinfo.frames);
+	if (p > percent || frame == 0) {
+	    percent = p;
+            if (!quiet) {
+                cerr << "\r" << percent << "% ";
+            }
+	}
+
+        frame += ibs;
+    }
+
+    if (!quiet) {
+        cerr << "\r    " << endl;
+    }
+    int avail;
+
+    while ((avail = ts.available()) >= 0) {
+
+        if (debug > 1) {
+            cerr << "(completing) available = " << avail << endl;
+        }
+
+        if (avail > 0) {
+            float **obf = new float *[channels];
+            for (size_t i = 0; i < channels; ++i) {
+                obf[i] = new float[avail];
+            }
+            ts.retrieve(obf, avail);
+            countOut += avail;
+            float *fobf = new float[channels * avail];
+            for (size_t c = 0; c < channels; ++c) {
+                for (int i = 0; i < avail; ++i) {
+                    float value = obf[c][i];
+                    if (value > 1.f) value = 1.f;
+                    if (value < -1.f) value = -1.f;
+                    fobf[i * channels + c] = value;
+                }
+            }
+
+            sf_writef_float(sndfileOut, fobf, avail);
+            delete[] fobf;
+            for (size_t i = 0; i < channels; ++i) {
+                delete[] obf[i];
+            }
+            delete[] obf;
+        } else {
+            usleep(10000);
+        }
+    }
+
+    sf_close(sndfile);
+    sf_close(sndfileOut);
+
+    if (!quiet) {
+
+        cerr << "in: " << countIn << ", out: " << countOut << ", ratio: " << float(countOut)/float(countIn) << ", ideal output: " << lrint(countIn * ratio) << ", error: " << abs(lrint(countIn * ratio) - int(countOut)) << endl;
+
+#ifdef _WIN32
+        RubberBand::
+#endif
+        timeval etv;
+        (void)gettimeofday(&etv, 0);
+        
+        etv.tv_sec -= tv.tv_sec;
+        if (etv.tv_usec < tv.tv_usec) {
+            etv.tv_usec += 1000000;
+            etv.tv_sec -= 1;
+        }
+        etv.tv_usec -= tv.tv_usec;
+        
+        double sec = double(etv.tv_sec) + (double(etv.tv_usec) / 1000000.0);
+        cerr << "elapsed time: " << sec << " sec, in frames/sec: " << countIn/sec << ", out frames/sec: " << countOut/sec << endl;
+    }
+
+    RubberBand::Profiler::dump();
+
+    return 0;
+}
+
+