samer@1: (import "samer.mds.*")
samer@1: 
samer@1: (put "colormap" (redgreen))
samer@1: (put "oz.background" (Color. 0 0 0))
samer@1: 
samer@1: (define N 16)
samer@1: (define D (Matrix. "D" N N))
samer@1: (define P (Matrix. "P" N 4))
samer@1: (define s (VVector. "s" N))
samer@1: (Mathx.set (.array s) (Constant. 0.5))
samer@1: 
samer@1: ;(MatrixBall.setBallSize 16 16)
samer@1: (define mds (MDS. P))
samer@1: (addtasks 
samer@1: 	(CovarianceTask. mds D)
samer@1: 	mds
samer@1: )
samer@1: (.setDimensionality mds 2)
samer@1: (expose)
samer@1: 
samer@1: (define gauss (samer.maths.random.NormalisedGaussian.))
samer@1: (define uniform (samer.maths.random.Uniform.))
samer@1: (define raleigh (samer.maths.random.Raleigh.))
samer@1: (define (rndP) (.set P gauss) (.changed P))
samer@1: (define (rndD) (.set D uniform) (.changed D))
samer@1: 
samer@1: ;(matexec D "load")
samer@1: (matexec D "image")
samer@1: ;(.setColumn (RowColumn. (Node. "Px") P) 0)
samer@1: ;(.setColumn (RowColumn. (Node. "Py") P) 1)
samer@1: (rndD)
samer@1: (rndP)
samer@1: 
samer@1: (define oz (node "oz" (MatrixPointViewer2. P s)))
samer@1: (expose oz "Oz")
samer@1: 
samer@1: