Analysis.h

Go to the documentation of this file.

//Copyright (c) 2004-2005, Baris Sumengen
//All rights reserved.
//
// CIMPL Matrix Performance Library
//
//Redistribution and use in source and binary
//forms, with or without modification, are
//permitted provided that the following
//conditions are met:
//
//    * No commercial use is allowed. 
//    This software can only be used
//    for non-commercial purposes. This 
//    distribution is mainly intended for
//    academic research and teaching.
//    * Redistributions of source code must
//    retain the above copyright notice, this
//    list of conditions and the following
//    disclaimer.
//    * Redistributions of binary form must
//    mention the above copyright notice, this
//    list of conditions and the following
//    disclaimer in a clearly visible part 
//    in associated product manual, 
//    readme, and web site of the redistributed 
//    software.
//    * Redistributions in binary form must
//    reproduce the above copyright notice,
//    this list of conditions and the
//    following disclaimer in the
//    documentation and/or other materials
//    provided with the distribution.
//    * The name of Baris Sumengen may not be
//    used to endorse or promote products
//    derived from this software without
//    specific prior written permission.
//
//THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT
//HOLDERS AND CONTRIBUTORS "AS IS" AND ANY
//EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT
//NOT LIMITED TO, THE IMPLIED WARRANTIES OF
//MERCHANTABILITY AND FITNESS FOR A PARTICULAR
//PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE
//CONTRIBUTORS BE LIABLE FOR ANY
//DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
//EXEMPLARY, OR CONSEQUENTIAL DAMAGES
//(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT
//OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
//DATA, OR PROFITS; OR BUSINESS INTERRUPTION)
//HOWEVER CAUSED AND ON ANY THEORY OF
//LIABILITY, WHETHER IN CONTRACT, STRICT
//LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR
//OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
//OF THIS SOFTWARE, EVEN IF ADVISED OF THE
//POSSIBILITY OF SUCH DAMAGE.



#pragma once
#ifndef ANALYSIS_H
#define ANALYSIS_H

#include "./cimpl.h"
using namespace CIMPL;

#include "./IO.h"


#include <time.h>

using namespace MathCore;
using namespace IO;


namespace Analysis
{
        void StatusDisplay(long status);

        
        enum Border
        {
                ZeroPad,
                Symmetric,
                Replicate,
                Circular
        };

        enum FilterArea
        {
                Full,
                Same,
                Valid
        };

        Vector<ComplexFloat> FFT(Vector<ComplexFloat>& x);
        Vector<ComplexFloat>& FFTI(Vector<ComplexFloat>& x);
        Vector<ComplexDouble> FFT(Vector<ComplexDouble>& x);
        Vector<ComplexDouble>& FFTI(Vector<ComplexDouble>& x);
        Vector<ComplexFloat> FFT(Vector<float>& x);
        Vector<ComplexDouble> FFT(Vector<double>& x);

        Vector<ComplexFloat> IFFT(Vector<ComplexFloat>& x);
        Vector<ComplexFloat>& IFFTI(Vector<ComplexFloat>& x);
        Vector<ComplexDouble> IFFT(Vector<ComplexDouble>& x);
        Vector<ComplexDouble>& IFFTI(Vector<ComplexDouble>& x);

        Matrix<ComplexFloat> FFT2(Matrix<ComplexFloat>& x);
        Matrix<ComplexFloat>& FFT2I(Matrix<ComplexFloat>& x);
        Matrix<ComplexDouble> FFT2(Matrix<ComplexDouble>& x);
        Matrix<ComplexDouble>& FFT2I(Matrix<ComplexDouble>& x);

        Matrix<ComplexFloat> IFFT2(Matrix<ComplexFloat>& x);
        Matrix<ComplexFloat>& IFFT2I(Matrix<ComplexFloat>& x);
        Matrix<ComplexDouble> IFFT2(Matrix<ComplexDouble>& x);
        Matrix<ComplexDouble>& IFFT2I(Matrix<ComplexDouble>& x);


        Vector<float> Conv(Vector<float>& input, Vector<float>& filter, FilterArea fa = Full, Border b = ZeroPad, bool PowerOfTwo = true);
        Vector<double> Conv(Vector<double>& input, Vector<double>& filter, FilterArea fa = Full, Border b = ZeroPad, bool PowerOfTwo = true);
        Vector<ComplexFloat> Conv(Vector<ComplexFloat>& input, Vector<ComplexFloat>& filter, FilterArea fa = Full, Border b = ZeroPad, bool PowerOfTwo = true);
        Vector<ComplexDouble> Conv(Vector<ComplexDouble>& input, Vector<ComplexDouble>& filter, FilterArea fa = Full, Border b = ZeroPad, bool PowerOfTwo = true);


        Matrix<float> Conv2(Matrix<float>& input, Matrix<float>& filter, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
        Matrix<float> Conv2(Vector<float>& filter1, Vector<float>& filter2, Matrix<float>& input, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
        Matrix<double> Conv2(Matrix<double>& input, Matrix<double>& filter, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
        Matrix<double> Conv2(Vector<double>& filter1, Vector<double>& filter2, Matrix<double>& input, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
        Matrix<ComplexFloat> Conv2(Matrix<ComplexFloat>& input, Matrix<ComplexFloat>& filter, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
        Matrix<ComplexFloat> Conv2(Vector<ComplexFloat>& filter1, Vector<ComplexFloat>& filter2, Matrix<ComplexFloat>& input, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
        Matrix<ComplexDouble> Conv2(Matrix<ComplexDouble>& input, Matrix<ComplexDouble>& filter, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);
        Matrix<ComplexDouble> Conv2(Vector<ComplexDouble>& filter1, Vector<ComplexDouble>& filter2, Matrix<ComplexDouble>& input, FilterArea fa = Same, Border b=Symmetric, bool PowerOfTwo = true);


        //FFTN();
        //IFFTN();
        //
        //ConvN();
        //Deconv();
        //Filter();
        //Filter2();
        //
        //Cov();
        //Corr2();
        //CorrCoef();
        //XCorr();
        //XCorr2();
        

};







#endif

---