source: trunk/SRC/ToBeReviewed/CALCULS/determ2.pro @ 142

;+
;
; file_comments
; computes the determinant of n 2 by 2 arrays
; Z2DS is an 2*2*n array
;
; @categories 
; Without loop
; 
; @param A {in}{required}
; n element array
; Defined as z2ds[0, 0, *] = z1d00
;
; @param B {in}{required}
; n element array
; Defined as z2ds[0, 1, *] = z1d01
;
; @param C {in}{required}
; n element array
; Defined as z2ds[1, 0, *] = z1d10
;
; @param D {in}{required}
; n element array
; Defined as z2ds[1, 1, *] = z1d11
;
; @returns 
; n elements array, the determinent of each 2*2 arrrays
;
; @examples
; a=findgen(2,2,5)
; print, determ2(a)
; FOR i=0,4 DO print, determ(a[*,*,i]) ; IDL solution
;
; @history
;           S. Masson (smasson@lodyc.jussieu.fr)
;           July 11th, 2002
;
; @version
; $Id$
;
;-
FUNCTION determ2, a, b, c, d
;
  compile_opt idl2, strictarrsubs
;
  CASE n_params() OF
    1:res = a[0, 0, *]*a[1, 1, *]-a[0, 1, *]*a[1, 0, *]
    4:res = a[*]*d[*]-c[*]*b[*]
    ELSE:stop
  ENDCASE

  RETURN, res
END