Changeset 31 for trunk/ToBeReviewed/MATRICE/colle.pro

Timestamp:

05/02/06 15:54:11 (18 years ago)

Author:

pinsard

Message:

upgrade of MATRICE according to cerbere.lodyc.jussieu.fr: /usr/home/smasson/SAXO_RD/ : files

File:

• trunk/ToBeReviewed/MATRICE/colle.pro (copied) (copied from trunk/MATRICE/colle.pro) (1 diff)

trunk/ToBeReviewed/MATRICE/colle.pro

@@ -133 +133 @@
 ; pour cela on contient le vecteur permute qui donne la place que
 ; doivent prendre les dimensions ds la matrice transposee
+        siz = (size(*ptrtab[0]))[0]
+         if siz LT direc then $
+          *ptrtab[0] = reform(*ptrtab[0], [(size(*ptrtab[0]))[1:siz], replicate(1, direc-siz)], /over)
          permute = indgen((size(*ptrtab[0]))[0])
          permute[0] = direc-1
          permute[direc-1] = 0
-         if (size(*ptrtab[0]))[0] NE direc then $
-          *ptrtab[0] = reform(*ptrtab[0], [(size(*ptrtab[0]))[1:direc-1],  1], /over)
          res = transpose(*ptrtab[0], permute)
          if NOT keyword_set(sauve) then ptr_free, ptrtab[0]
          FOR n = 1,nbretab-1 DO BEGIN ; on colle suivant la dimension 1
-            if (size(*ptrtab[n]))[0] NE direc then $
-             *ptrtab[n] = reform(*ptrtab[n], [(size(*ptrtab[n]))[1:direc-1],  1])            
+            if (size(*ptrtab[n]))[0] LT direc then $
+             *ptrtab[n] = reform(*ptrtab[n], [(size(*ptrtab[n]))[1:siz], replicate(1, direc-siz)])            
             res = [temporary(res), transpose(*ptrtab[n], permute)]
             if NOT keyword_set(sauve) then ptr_free, ptrtab[n]