Changeset 2585 in Sophya for trunk/SophyaLib/TArray/tmatrix.cc

Timestamp:

Jul 30, 2004, 3:11:51 PM (21 years ago)

Author:

cmv

Message:

Correction bug Multiply segment gault cmv 30/07/04

File:

• trunk/SophyaLib/TArray/tmatrix.cc (modified) (9 diffs)

trunk/SophyaLib/TArray/tmatrix.cc

@@ -1 @@
-// $Id: tmatrix.cc,v 1.27 2004-07-30 10:24:12 cmv Exp $
+// $Id: tmatrix.cc,v 1.28 2004-07-30 13:11:51 cmv Exp $
 //                         C.Magneville          04/99
 #include "machdefs.h"
@@ -405 +405 @@
 }
 
-////////////////////////////////////////////////////////////////
-//****  Multiplication matricielle  *****
-////////////////////////////////////////////////////////////////
+//////////////////////////////////////////////////////////
+/////////////// Multiplication matricielle ///////////////
+//////////////////////////////////////////////////////////
 
 //! Return the matrix product C = (*this)*B
@@ -454 +454 @@
 
   // Commentaire: pas de difference de vitesse notable selon le mapping de la matrice produit "rm"
-  if (mm == SameMemoryMapping) mm = GetMemoryMapping();  
+  if (mm == SameMemoryMapping) mm = GetMemoryMapping();
   TMatrix<T> rm(NRows(), b.NCols(), mm);
 
@@ -515 +515 @@
 
   if(stepa==1 && stepb==1) {
-    //cout<<"A.col packed && B.row not packed "<<stepa<<" "<<stepb<<endl;
-     T * pea = new T[rm.NCols()];
+    //cout<<"A.col packed && B.row packed "<<stepa<<" "<<stepb<<endl;
+     T * pea = new T[NCols()];
     const T * peb;
     for(r=0; r<rm.NRows(); r++) {     // Boucle sur les lignes de A
@@ -540 +540 @@
     //cout<<"A.col packed && B.row not packed "<<stepa<<" "<<stepb<<endl;
     const T * pea;
-    T * peb = new T[rm.NCols()];
+    T * peb = new T[NCols()];
     for(c=0; c<rm.NCols(); c++) {     // Boucle sur les colonnes de B
       pe = &(b(0,c));
@@ -562 +562 @@
   if(stepa!=1 && stepb==1) {
     //cout<<"A.col not packed && B.row packed "<<stepa<<" "<<stepb<<endl;
-    T * pea = new T[rm.NCols()];
+    T * pea = new T[NCols()];
     const T * peb;
     for(r=0; r<rm.NRows(); r++) {     // Boucle sur les lignes de A
@@ -608 +608 @@
     //cout<<"A.col not packed && B.row not packed ==> optimize on B "<<stepa<<" "<<stepb<<endl;
     const T * pea;
-    T * peb = new T[rm.NCols()];
+    T * peb = new T[NCols()];
     for(c=0; c<rm.NCols(); c++) {     // Boucle sur les colonnes de B
       pe = &(b(0,c));
@@ -627 +627 @@
   if(stepa!=1 && stepb!=1) {
     //cout<<"A.col not packed && B.row not packed ==> optimize on A "<<stepa<<" "<<stepb<<endl;
-    T * pea = new T[rm.NCols()];
+    T * pea = new T[NCols()];
     const T * peb;
     for(r=0; r<rm.NRows(); r++) {     // Boucle sur les lignes de A
@@ -650 +650 @@
   throw(SzMismatchError("TMatrix<T>::Multiply(b) Optimize case not treated... Please report BUG !!! ") );
   return rm;
-}
+
+}
+
 
 ///////////////////////////////////////////////////////////////