/**
* Combine bias data from small files to a integrated file
*
* @param conf
* @param pathName
* @param outName
* @throws IOException
*/
public static void extractBias(Configuration conf, String pathName, String outName, int n)
throws IOException {
List<String> resList = DataOperators.readTextFromHDFS(conf, pathName);
DenseDoubleMatrix1D vBias = new DenseDoubleMatrix1D(n);
for (String line : resList) {
String[] items = line.split("\t");
vBias.set(Integer.parseInt(items[0]), Double.parseDouble(items[1]));
}
MatrixIO.saveDenseMatrix1D2HDFS(conf, outName, vBias, true);
}
public static void main(String[] args) {
double[][] arg1 = {{1, 1, 1}, {2, 3, 2}, {4, 7, 3}};
DenseDoubleMatrix2D obj1 = new DenseDoubleMatrix2D(arg1); // rot
double[][] arg2 = {{5, 5, 5}, {8, 8, 8}, {2, 2, 2}};
DenseDoubleMatrix2D obj2 = new DenseDoubleMatrix2D(arg2); // grŸn
double[][] arg3 = {{4, 5, 6}, {3, 2, 1}, {1, 2, 1}};
DenseDoubleMatrix2D obj3 = new DenseDoubleMatrix2D(arg3); // schwarz
System.out.println("obj1\n" + obj1);
System.out.println("obj2\n" + obj2);
System.out.println("obj3\n" + obj3);
double[][][] pair21 = {arg1, arg2};
Transformation trans21 = Kabsch.calculateTransformation(pair21); // grŸn -> rot
DenseDoubleMatrix2D rot21 = trans21.peekRotation();
DenseDoubleMatrix1D trl21 = trans21.peekTranslation();
System.out.println("rot21\n" + rot21);
System.out.println("trl21\n" + trl21);
double[][][] pair12 = {arg2, arg1};
Transformation trans12 = Kabsch.calculateTransformation(pair12); // grŸn -> rot
DenseDoubleMatrix2D rot12 = trans12.peekRotation();
DenseDoubleMatrix1D trl12 = trans12.peekTranslation();
System.out.println("rot12\n" + rot12);
System.out.println("trl12\n" + trl12);
double[][][] pair32 = {arg2, arg3};
Transformation trans32 = Kabsch.calculateTransformation(pair32); // schwarz -> grŸn
DenseDoubleMatrix2D rot32 = trans32.peekRotation();
DenseDoubleMatrix1D trl32 = trans32.peekTranslation();
System.out.println("rot32\n" + rot32);
System.out.println("trl32\n" + trl32.toString());
DenseDoubleMatrix2D arg2onarg1 =
new DenseDoubleMatrix2D(trans21.transform(arg2)); // grŸn -> rot
System.out.println(arg2onarg1.toString());
DenseDoubleMatrix2D uni21 = createUniMatrix(rot21, trl21);
DenseDoubleMatrix2D uni32 = createUniMatrix(rot32, trl32);
System.out.println("uni21\n" + uni21.toString());
System.out.println("uni32\n" + uni32.toString());
Algebra alg = new Algebra();
DenseDoubleMatrix2D result =
(DenseDoubleMatrix2D)
alg.mult(
DoubleFactory2D.dense.identity(4),
computeTransitiveTransformation(
(DenseDoubleMatrix2D) DoubleFactory2D.dense.identity(4), uni32, uni21));
System.out.println("action results\n" + result);
System.out.println("moved points\n" + movePoints(result, obj3).toString());
}
public static DenseDoubleMatrix2D createUniMatrix(
DenseDoubleMatrix2D rot, DenseDoubleMatrix1D trl) {
DenseDoubleMatrix2D res = new DenseDoubleMatrix2D(4, 4);
for (int i = 0; i != 3; i++) {
for (int j = 0; j != 3; j++) {
res.set(i, j, rot.get(i, j));
}
}
for (int i = 0; i != 3; i++) {
res.set(i, 3, trl.get(i));
}
for (int i = 0; i != 3; i++) {
res.set(3, i, 0);
}
res.set(3, 3, 1);
return res;
}
