optimization
parent
2e919a82d1
commit
b843de7707
|
@ -1,9 +1,8 @@
|
||||||
! Copyright (C) 2005 Slava Pestov.
|
! Copyright (C) 2005, 2006 Slava Pestov.
|
||||||
! See http://factor.sf.net/license.txt for BSD license.
|
! See http://factorcode.org/license.txt for BSD license.
|
||||||
IN: math
|
IN: math
|
||||||
USING: arrays generic kernel sequences ;
|
USING: arrays generic kernel sequences ;
|
||||||
|
|
||||||
! Vectors
|
|
||||||
: vneg ( v -- v ) [ neg ] map ;
|
: vneg ( v -- v ) [ neg ] map ;
|
||||||
|
|
||||||
: n*v ( n v -- v ) [ * ] map-with ;
|
: n*v ( n v -- v ) [ * ] map-with ;
|
||||||
|
@ -18,22 +17,13 @@ USING: arrays generic kernel sequences ;
|
||||||
: vmax ( v v -- v ) [ max ] 2map ;
|
: vmax ( v v -- v ) [ max ] 2map ;
|
||||||
: vmin ( v v -- v ) [ min ] 2map ;
|
: vmin ( v v -- v ) [ min ] 2map ;
|
||||||
|
|
||||||
|
: v. ( v v -- x ) 0 [ * + ] 2reduce ;
|
||||||
|
: norm-sq ( v -- n ) 0 [ absq + ] reduce ;
|
||||||
|
: norm ( vec -- n ) norm-sq sqrt ;
|
||||||
|
: normalize ( vec -- uvec ) dup norm v/n ;
|
||||||
|
|
||||||
: sum ( v -- n ) 0 [ + ] reduce ;
|
: sum ( v -- n ) 0 [ + ] reduce ;
|
||||||
: product ( v -- n ) 1 [ * ] reduce ;
|
: product ( v -- n ) 1 [ * ] reduce ;
|
||||||
|
|
||||||
: set-axis ( x y axis -- v )
|
: set-axis ( x y axis -- v )
|
||||||
2dup v* >r >r drop dup r> v* v- r> v+ ;
|
dup length [ >r 0 = pick pick ? r> swap nth ] 2map 2nip ;
|
||||||
|
|
||||||
: v. ( v v -- x )
|
|
||||||
#! Dot product.
|
|
||||||
0 [ * + ] 2reduce ;
|
|
||||||
|
|
||||||
: norm-sq ( v -- n ) 0 [ absq + ] reduce ;
|
|
||||||
|
|
||||||
: norm ( vec -- n )
|
|
||||||
#! Length of a vector.
|
|
||||||
norm-sq sqrt ;
|
|
||||||
|
|
||||||
: normalize ( vec -- uvec )
|
|
||||||
#! Unit vector with same direction as vec.
|
|
||||||
dup norm v/n ;
|
|
||||||
|
|
Loading…
Reference in New Issue