File: mtlb_diff.sci

package info (click to toggle)
scilab 5.2.2-9
  • links: PTS, VCS
  • area: main
  • in suites: squeeze
  • size: 334,832 kB
  • ctags: 52,586
  • sloc: xml: 526,945; ansic: 223,590; fortran: 163,080; java: 56,934; cpp: 33,840; tcl: 27,936; sh: 20,397; makefile: 9,908; ml: 9,451; perl: 1,323; cs: 614; lisp: 30
file content (46 lines) | stat: -rw-r--r-- 945 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) 2002-2004 - INRIA - Vincent COUVERT 
// 
// This file must be used under the terms of the CeCILL.
// This source file is licensed as described in the file COPYING, which
// you should have received as part of this distribution.  The terms
// are also available at    
// http://www.cecill.info/licences/Licence_CeCILL_V2-en.txt

function y=mtlb_diff(x,N,dim)
// Emulation function for Matlab diff()

[lhs,rhs]=argn()
if rhs==1 then N=1,end
if rhs<=2 then
  if length(size(x))>2 then
    dim=firstnonsingleton(x)
  elseif size(x,1)==1 | size(x,2)==1 then
    dim=0
  else
    dim='r'
  end
    
end

if dim==0 then
  y=diff(x,N)
else
  if N>=size(x,"*") then
    y=[]
    return
  else
    for k=1:N
      x=diff(x,1,dim)
      if size(x,dim)==1 then
	dim=firstnonsingleton(x)
      end
    end
  end
  if x==[] then
    y=0
  else
    y=x
  end
end
endfunction