package edu.stanford.rsl.conrad.volume3d.operations;
public class MultiplySlabs extends ParallelVolumeOperation {
@Override
public void performOperation() {
switch (vol1.getInternalDimension()) {
case 1:
for (int indexX=beginIndexX; indexX<endIndexX; indexX++) {
for (int indexY=0; indexY<vol1.size[1]; indexY++) {
for (int indexZ=0; indexZ<vol1.size[2]; indexZ++) {
vol1.data[indexX][indexY][indexZ] *=
vol2.data[indexX][indexY][indexZ];
}
}
}
break;
case 2:
for (int indexX=beginIndexX; indexX<endIndexX; indexX++) {
for (int indexY=0; indexY<vol1.size[1]; indexY++) {
for (int indexZ=0; indexZ<vol1.size[2]; indexZ++) {
float tmp_re1 = vol1.data[indexX][indexY][indexZ*2];
float tmp_re2 = vol2.data[indexX][indexY][indexZ*2];
vol1.data[indexX][indexY][indexZ*2] =
vol1.data[indexX][indexY][indexZ*2] * vol2.data[indexX][indexY][indexZ*2]-
vol1.data[indexX][indexY][indexZ*2+1] * vol2.data[indexX][indexY][indexZ*2+1];
vol1.data[indexX][indexY][indexZ*2+1] =
vol1.data[indexX][indexY][indexZ*2+1]*tmp_re2+
tmp_re1*vol2.data[indexX][indexY][indexZ*2+1];
}
}
}
break;
default:
System.out.println("vol_mult: Invalid dimension\n");
}
}
@Override
public ParallelVolumeOperation clone() {
return new MultiplySlabs();
}
}
/*
* Copyright (C) 2010-2014 Andreas Maier
* CONRAD is developed as an Open Source project under the GNU General Public License (GPL).
*/