From: ageay Date: Wed, 27 Nov 2013 09:40:31 +0000 (+0000) Subject: For 0.11.0 scipy users. spdiags used instead of diags. X-Git-Tag: V7_3_0a1~2 X-Git-Url: http://git.salome-platform.org/gitweb/?a=commitdiff_plain;h=ffbbe9e7e9ec71b7e0416843fa1ed122ae6da97a;p=modules%2Fmed.git For 0.11.0 scipy users. spdiags used instead of diags. --- diff --git a/src/MEDCoupling_Swig/MEDCouplingRemapperTest.py b/src/MEDCoupling_Swig/MEDCouplingRemapperTest.py index 9c463ca08..9d4b26947 100644 --- a/src/MEDCoupling_Swig/MEDCouplingRemapperTest.py +++ b/src/MEDCoupling_Swig/MEDCouplingRemapperTest.py @@ -591,7 +591,7 @@ class MEDCouplingBasicsTest(unittest.TestCase): def testGetCrudeCSRMatrix1(self): """ testing CSR matrix output using numpy/scipy. """ - from scipy.sparse import diags + from scipy.sparse import spdiags #diags import scipy from numpy import array arr=DataArrayDouble(3) ; arr.iota() @@ -620,7 +620,8 @@ class MEDCouplingBasicsTest(unittest.TestCase): self.assertEqual(diff.getnnz(),0) # IntegralGlobConstraint (division by sum of cols) colSum=m.sum(axis=0) - m_0=m*diags(array(1/colSum),[0]) + # version 0.12.0 # m_0=m*diags(array(1/colSum),[0]) + m_0=m*spdiags(array(1/colSum),[0],colSum.shape[1],colSum.shape[1]) del colSum self.assertAlmostEqual(m_0[0,0],0.625,12) self.assertAlmostEqual(m_0[1,0],0.25,12) @@ -632,7 +633,8 @@ class MEDCouplingBasicsTest(unittest.TestCase): self.assertEqual(m_0.getnnz(),7) # ConservativeVolumic (division by sum of rows) rowSum=m.sum(axis=1) - m_1=diags(array(1/rowSum.transpose()),[0])*m + # version 0.12.0 # m_1=diags(array(1/rowSum.transpose()),[0])*m + m_1=spdiags(array(1/rowSum.transpose()),[0],rowSum.shape[0],rowSum.shape[0])*m del rowSum self.assertAlmostEqual(m_1[0,0],1.,12) self.assertAlmostEqual(m_1[1,0],0.4,12)