corecomplex.cpp

图像分割算法

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



#include "./CoreComplex.h"


namespace MathCore
{





	Matrix<ComplexFloat>& Rand(Matrix<ComplexFloat>& m, float max)
	{
		if(!RandomGen::Initialized())
		{
			RandomGen::Initialize();
		}
		for(int i=0;i<m.Length();i++)
		{
			m.ElemNC(i) = ComplexFloat( (float)(rand()/(double)RAND_MAX*max), (float)(rand()/(double)RAND_MAX*max) );
		}
		return m;
	}
	
	Vector<ComplexFloat>& Rand(Vector<ComplexFloat>& m, float max)
	{
		if(!RandomGen::Initialized())
		{
			RandomGen::Initialize();
		}
		for(int i=0;i<m.Length();i++)
		{
			m.ElemNC(i) = ComplexFloat( (float)(rand()/(double)RAND_MAX*max), (float)(rand()/(double)RAND_MAX*max) );
		}
		return m;
	}


	Matrix<ComplexDouble>& Rand(Matrix<ComplexDouble>& m, double max)
	{
		if(!RandomGen::Initialized())
		{
			RandomGen::Initialize();
		}
		for(int i=0;i<m.Length();i++)
		{
			m.ElemNC(i) = ComplexDouble( (double)(rand()/(double)RAND_MAX*max), (double)(rand()/(double)RAND_MAX*max) );
		}
		return m;
	}

	Vector<ComplexDouble>& Rand(Vector<ComplexDouble>& m, double max)
	{
		if(!RandomGen::Initialized())
		{
			RandomGen::Initialize();
		}
		for(int i=0;i<m.Length();i++)
		{
			m.ElemNC(i) = ComplexDouble( (double)(rand()/(double)RAND_MAX*max), (double)(rand()/(double)RAND_MAX*max) );
		}
		return m;
	}





	Vector<float> Abs(Vector<ComplexFloat>& m)
	{
		Vector<float> dummy(m.Numel());
		for(int i=0; i<m.Numel(); i++)
		{
			dummy.ElemNC(i) = abs(m.ElemNC(i));
		}
		return dummy;
	}

	Matrix<float> Abs(Matrix<ComplexFloat>& m)
	{
		Matrix<float> dummy(m.Rows(),m.Columns());
		for(int i=0; i<m.Numel(); i++)
		{
			dummy.ElemNC(i) = abs(m.ElemNC(i));
		}
		return dummy;
	}

	Vector<double> Abs(Vector<ComplexDouble>& m)
	{
		Vector<double> dummy(m.Numel());
		for(int i=0; i<m.Numel(); i++)
		{
			dummy.ElemNC(i) = abs(m.ElemNC(i));
		}
		return dummy;
	}

	Matrix<double> Abs(Matrix<ComplexDouble>& m)
	{
		Matrix<double> dummy(m.Rows(),m.Columns());
		for(int i=0; i<m.Numel(); i++)
		{
			dummy.ElemNC(i) = abs(m.ElemNC(i));
		}
		return dummy;
	}



	Vector<float> Angle(Vector<ComplexFloat>& m)
	{
		Vector<float> dummy(m.Numel());
		for(int i=0; i<m.Numel(); i++)
		{
			dummy.ElemNC(i) = arg(m.ElemNC(i));
		}
		return dummy;
	}

	Matrix<float> Angle(Matrix<ComplexFloat>& m)
	{
		Matrix<float> dummy(m.Rows(),m.Columns());
		for(int i=0; i<m.Numel(); i++)
		{
			dummy.ElemNC(i) = arg(m.ElemNC(i));
		}
		return dummy;
	}

	Vector<double> Angle(Vector<ComplexDouble>& m)
	{
		Vector<double> dummy(m.Numel());
		for(int i=0; i<m.Numel(); i++)
		{
			dummy.ElemNC(i) = arg(m.ElemNC(i));
		}
		return dummy;
	}

	Matrix<double> Angle(Matrix<ComplexDouble>& m)
	{
		Matrix<double> dummy(m.Rows(),m.Columns());
		for(int i=0; i<m.Numel(); i++)
		{
			dummy.ElemNC(i) = arg(m.ElemNC(i));
		}
		return dummy;
	}






	Vector<ComplexFloat> Conj(Vector<ComplexFloat>& m)
	{
		Vector<ComplexFloat> dummy(m.Numel());
		for(int i=0; i<m.Numel(); i++)
		{
			dummy.ElemNC(i) = conj(m.ElemNC(i));
		}
		return dummy;
	}

	Matrix<ComplexFloat> Conj(Matrix<ComplexFloat>& m)
	{
		Matrix<ComplexFloat> dummy(m.Rows(),m.Columns());
		for(int i=0; i<m.Numel(); i++)
		{
			dummy.ElemNC(i) = conj(m.ElemNC(i));
		}
		return dummy;
	}

	Vector<ComplexDouble> Conj(Vector<ComplexDouble>& m)
	{
		Vector<ComplexDouble> dummy(m.Numel());
		for(int i=0; i<m.Numel(); i++)
		{
			dummy.ElemNC(i) = conj(m.ElemNC(i));
		}
		return dummy;
	}

	Matrix<ComplexDouble> Conj(Matrix<ComplexDouble>& m)
	{
		Matrix<ComplexDouble> dummy(m.Rows(),m.Columns());
		for(int i=0; i<m.Numel(); i++)
		{
			dummy.ElemNC(i) = conj(m.ElemNC(i));
		}
		return dummy;
	}



	Vector<ComplexFloat>& ConjI(Vector<ComplexFloat>& m)
	{
		for(int i=0; i<m.Numel(); i++)
		{
			m.ElemNC(i) = conj(m.ElemNC(i));
		}
		return m;
	}

	Matrix<ComplexFloat>& ConjI(Matrix<ComplexFloat>& m)
	{
		for(int i=0; i<m.Numel(); i++)
		{
			m.ElemNC(i) = conj(m.ElemNC(i));
		}
		return m;
	}

	Vector<ComplexDouble>& ConjI(Vector<ComplexDouble>& m)
	{
		for(int i=0; i<m.Numel(); i++)
		{
			m.ElemNC(i) = conj(m.ElemNC(i));
		}
		return m;
	}

	Matrix<ComplexDouble>& ConjI(Matrix<ComplexDouble>& m)
	{
		for(int i=0; i<m.Numel(); i++)
		{
			m.ElemNC(i) = conj(m.ElemNC(i));
		}
		return m;
	}




	Vector<float> Real(Vector<ComplexFloat> &m)
	{
		int l = m.Numel();
		Vector<float> temp(l);
		for(int i=0; i<l; i++)
		{
			temp[i] = real(m[i]);
		}
		return temp;
	}

	Vector<float> Imag(Vector<ComplexFloat> &m)
	{
		int l = m.Numel();
		Vector<float> temp(l);
		for(int i=0; i<l; i++)
		{
			temp[i] = imag(m[i]);
		}
		return temp;
	}

	Matrix<float> Real(Matrix<ComplexFloat> &m)
	{
		int r = m.Rows();
		int c = m.Columns();
		int l = m.Numel();
		Matrix<float> temp(r,c);
		for(int i=0; i<l; i++)
		{
			temp.ElemNC(i) = real(m.ElemNC(i));
		}
		return temp;
	}

	Matrix<float> Imag(Matrix<ComplexFloat> &m)
	{
		int r = m.Rows();
		int c = m.Columns();
		int l = m.Numel();
		Matrix<float> temp(r,c);
		for(int i=0; i<l; i++)
		{
			temp.ElemNC(i) = imag(m.ElemNC(i));
		}
		return temp;
	}





	Vector<double> Real(Vector<ComplexDouble> &m)
	{
		int l = m.Numel();
		Vector<double> temp(l);
		for(int i=0; i<l; i++)
		{
			temp[i] = real(m[i]);
		}
		return temp;
	}

	Vector<double> Imag(Vector<ComplexDouble> &m)
	{
		int l = m.Numel();
		Vector<double> temp(l);
		for(int i=0; i<l; i++)
		{
			temp[i] = imag(m[i]);
		}
		return temp;
	}

	Matrix<double> Real(Matrix<ComplexDouble> &m)
	{
		int r = m.Rows();
		int c = m.Columns();
		int l = m.Numel();
		Matrix<double> temp(r,c);
		for(int i=0; i<l; i++)
		{
			temp.ElemNC(i) = real(m.ElemNC(i));
		}
		return temp;
	}

	Matrix<double> Imag(Matrix<ComplexDouble> &m)
	{
		int r = m.Rows();
		int c = m.Columns();
		int l = m.Numel();
		Matrix<double> temp(r,c);
		for(int i=0; i<l; i++)
		{
			temp.ElemNC(i) = imag(m.ElemNC(i));
		}
		return temp;
	}










};