annotate src/samer/mds/CorrelationTask.java @ 8:5e3cbbf173aa tip

Reorganise some more
author samer
date Fri, 05 Apr 2019 22:41:58 +0100
parents bf79fb79ee13
children
rev   line source
samer@0 1 package samer.mds;
samer@0 2
samer@0 3 // import samer.core.*;
samer@0 4 import samer.maths.*;
samer@0 5 import samer.tools.*;
samer@0 6
samer@0 7 /**
samer@0 8 Transfer correlation matrix to MDS distances. Assumes 1s
samer@0 9 down the main diagonal and a symmetric matrix.
samer@0 10 */
samer@0 11 public class CorrelationTask extends AnonymousTask
samer@0 12 {
samer@0 13 int N;
samer@0 14 double [] d; // linear array of distances
samer@0 15 double [][] _R; // matrix of correlation coefficients
samer@0 16
samer@0 17 /** link each object to all the others using distances in matrix, returns a task
samer@0 18 that can be used to refresh distances from original matrix */
samer@0 19 public CorrelationTask(MDSBase mds, Matrix R) {
samer@0 20 N = R.getRowDimension();
samer@0 21
samer@0 22 d=new double[N*(N - 1)/2];
samer@0 23 _R=R.getArray();
samer@0 24
samer@0 25 mds.clearLinks(d);
samer@0 26 for (int k=0, i=0; i<N; i++)
samer@0 27 for (int j=0; j<i; j++) mds.setLink(k++,i,j);
samer@0 28
samer@0 29 run();
samer@0 30 }
samer@0 31
samer@0 32 public Vec getDistances() { return new Vec.ForArray(d); }
samer@0 33 public void run() {
samer@0 34 for (int k=0, i=0; i<N; i++) {
samer@0 35 double [] Ri=_R[i];
samer@0 36 for (int j=0; j<i; j++)
samer@0 37 d[k++]=Math.sqrt(-Math.log(Math.abs(Ri[j])));
samer@0 38 }
samer@0 39 }
samer@0 40 }