source: trunk/SRC/Obsolete/extrait.pro @ 134

;------------------------------------------------------------
;------------------------------------------------------------
;------------------------------------------------------------
;+
;
; @fie-comments 
; extractoin of subdomains of matrixes; Even if the subdomain is "pierced" (see the example)
; By default, IDL can make extractions of subdomain:
;
;      IDL> a=indgen(5,5) 
;      IDL> print, a 
;             0       1       2       3       4
;             5       6       7       8       9
;            10      11      12      13      14
;            15      16      17      18      19
;            20      21      22      23      24
;      IDL> print, a[[0,2],3] 
;            15      17
;      IDL> print, a[[0,2],*]  
;             0       2
;             5       7
;            10      12
;            15      17
;            20      22
; but
;      IDL> print, a[[0,2],[3,4]] 
;            15      22
; while
;      IDL> print, extrait(a,[0,2],[3,4])  
;            15      17
;            20      22
;
;
; you better use extrac2
;
; @obsolete
;
; @categories utilities
;
; @param tab {in}{required} a 1,2,3 or 4 dim table
;
; @param indicex {in}{required} can have 2 forms:
;              1)a vector containing indexes of lines we want to keep
;              2)the string '*' if we want to keep all lines.
;
; @param indicey {in}{required} the same thing that indicex but for dim 2.
;
; @param indicez {in}{required} the same thing that indicex but for dim 3.
; 
; @param indicet {in}{required} the same thing that indicex but for dim 4.
;
; @returns a matrix 1,2,3 or 4d extract from tab
;
; @restrictions res=-1 in case of mistake
;
; 
; @examples I have a dim 2 matrix named A. I want extract a small intersection 
;          matrix 2d of the line 2,3 and 7 and of the column 0 and 1:
;      
;      res=extrait(A,[2,3,7],[0,1])
;
;other ex:
;      IDL> print, a 
;      a b c
;      d e f
;      g h i
;      IDL> print, extrait(a,[0,2],[0,2])        
;      a c
;      g i
;
; @history Sebastien Masson (smasson@lodyc.jussieu.fr)
;                       12/1/1999
;                       29/4/1999: correction of a bug and complement of the heading
;
; @version $Id$
;
;-
;------------------------------------------------------------
;------------------------------------------------------------
;------------------------------------------------------------
FUNCTION extrait, tab, indicex, indicey, indicez, indicet
;------------------------------------------------------------
;
  compile_opt idl2, strictarrsubs
;
  case n_params() of
      0:return, extrac2()
      1:return, extrac2(tab)
      2:return, extrac2(tab, indicex)
      3:return, extrac2(tab, indicex, indicey)
      4:return, extrac2(tab, indicex, indicey, indicez)
      5:return, extrac2(tab, indicex, indicey, indicez, indicet)
  endcase
end