convert ipiv to selection vector

Author: bilam

lu.ijs

@@ -18,5 +18,11 @@ for_i. i. d do.  NB. for each ring
   y =. col (<remvec;i)} y  NB. store the current column into the array as L; leave U unchanged
   y =. (((<remvec;remvec) { y) - col */ (<i;remvec) { y) (<remvec;remvec)} y  NB. subtract (current column */ (}.first row)) from lower-right block
 end.
-pv;y
-)
+NB. convert ipiv to selection vector
+z=. i.#pv
+for_i. i.-#pv do.
+ t=. i{z
+ z=. t (i{pv)} ((i{pv){z) i}z
+end.
+z;y
+)

manifest.ijs

@@ -6,7 +6,7 @@ DESCRIPTION=: 0 : 0
 J LU code, called from 128!:10
 )
 
-VERSION=: '1.0.3'
+VERSION=: '1.0.4'
 
 FILES=: 0 : 0
 lu.ijs