From dc0b823c3ae508895ff43552ad11e81aed51a32a Mon Sep 17 00:00:00 2001
From: Jiri Pittner <jiri@pittnerovi.com>
Date: Thu, 11 Aug 2022 15:52:58 +0200
Subject: [PATCH] realpart() for sparsemat

---
 nonclass.cc | 23 +++++++++++++++++++++++
 1 file changed, 23 insertions(+)

diff --git a/nonclass.cc b/nonclass.cc
index 3a18a08..14690ca 100644
--- a/nonclass.cc
+++ b/nonclass.cc
@@ -1317,6 +1317,29 @@ for(int i=0; i<w.size(); ++i)
 	}
 }
 
+template<>
+const SparseMat<double> realpart<SparseMat< std::complex<double> > >(const SparseMat< std::complex<double> > &a)
+{
+SparseMat<double> result(a.nrows(), a.ncols());
+for(SparseMat< std::complex<double> >::iterator p=a.begin(); p!=a.end(); ++p)
+	{
+	result.add(p->row,p->col,p->elem.real());
+	}
+if(a.issymmetric()) result.setsymmetric();
+return result;
+}
+
+template<>
+const SparseMat<double> imagpart<SparseMat< std::complex<double> > >(const SparseMat< std::complex<double> > &a)
+{
+SparseMat<double> result(a.nrows(), a.ncols());
+for(SparseMat< std::complex<double> >::iterator p=a.begin(); p!=a.end(); ++p)
+	{
+	result.add(p->row,p->col,p->elem.imag());
+	}
+if(a.issymmetric()) result.setsymmetric();
+return result;
+}
 
 template<>
 const NRMat<double> realpart<NRMat< std::complex<double> > >(const NRMat< std::complex<double> > &a)