Commit 94525844 authored by Alvarez, Gonzalo's avatar Alvarez, Gonzalo
Browse files

ContinuedFraction with complex weight

parent b8ac1889

src/ContinuedFraction.h

+5 −3
Original line number Diff line number Diff line
@@ -180,7 +180,8 @@ public: 
	ComplexType iOfOmega(const ComplexType& z,RealType offset,int isign) const



	{

		if (weight_==0) return ComplexType(0,0);

		if (PsimagLite::real(weight_) == 0 && PsimagLite::imag(weight_) == 0)

			return ComplexType(0,0);



		ComplexType sum = 0;

		for (SizeType l=0;l<intensity_.size();l++)
@@ -197,7 +198,8 @@ private: 


	void diagonalize()

	{

		if (weight_==0) return;

		if (PsimagLite::real(weight_) == 0 && PsimagLite::imag(weight_) == 0) return;



		MatrixType T;

		ab_.buildDenseMatrix(T);

		eigs_.resize(T.rows());
@@ -217,7 +219,7 @@ private: 
	FreqEnum freqEnum_;

	TridiagonalMatrixType ab_;

	RealType Eg_;

	RealType weight_;

	std::complex<RealType> weight_;

	int isign_;

	typename Vector<RealType>::Type eigs_;

	typename Vector<RealType>::Type intensity_;

src/Io/IoSimple.h

+2 −1
Original line number Diff line number Diff line
@@ -95,7 +95,8 @@ class IoSimple { 


	template<typename T>

	struct PrintWithEqualSign {

		enum { True = Loki::TypeTraits<T>::isArith || std::is_enum<T>::value };

		enum { True = Loki::TypeTraits<T>::isArith || std::is_enum<T>::value ||

		     IsComplexNumber<T>::True};

	};



public: