/*\r
* Jalview - A Sequence Alignment Editor and Viewer\r
-* Copyright (C) 2005 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle\r
+* Copyright (C) 2006 AM Waterhouse, J Procter, G Barton, M Clamp, S Searle\r
*\r
* This program is free software; you can redistribute it and/or\r
* modify it under the terms of the GNU General Public License\r
* along with this program; if not, write to the Free Software\r
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA\r
*/\r
-\r
package jalview.math;\r
\r
\r
-\r
-public class RotatableMatrix {\r
-\r
- float matrix[][];\r
-\r
- float[] temp;\r
-\r
- float[][] rot;\r
-\r
-\r
-\r
- public RotatableMatrix(int rows, int cols) {\r
-\r
- matrix = new float[rows][cols];\r
-\r
- temp = new float[3];\r
-\r
- rot = new float[3][3];\r
-\r
- }\r
-\r
-\r
-\r
- public void addElement(int i, int j, float value) {\r
-\r
- matrix[i][j] = value;\r
-\r
- }\r
-\r
-\r
-\r
- public void print() {\r
-\r
- System.out.println(matrix[0][0] + " " + matrix[0][1] + " " + matrix[0][2]);\r
-\r
- System.out.println(matrix[1][0] + " " + matrix[1][1] + " " + matrix[1][2]);\r
-\r
- System.out.println(matrix[2][0] + " " + matrix[2][1] + " " + matrix[2][2]);\r
-\r
- }\r
-\r
-\r
-\r
- public void rotate (float degrees, char axis) {\r
-\r
-\r
-\r
- float costheta = (float)Math.cos(degrees*Math.PI/(float)180.0);\r
-\r
- float sintheta = (float)Math.sin(degrees*Math.PI/(float)180.0);\r
-\r
-\r
-\r
- if (axis == 'z') {\r
-\r
-\r
-\r
- rot[0][0] = (float)costheta;\r
-\r
- rot[0][1] = (float)-sintheta;\r
-\r
- rot[0][2] = (float)0.0;\r
-\r
-\r
-\r
- rot[1][0] = (float)sintheta;\r
-\r
- rot[1][1] = (float)costheta;\r
-\r
- rot[1][2] = (float)0.0;\r
-\r
-\r
-\r
- rot[2][0] = (float)0.0;\r
-\r
- rot[2][1] = (float)0.0;\r
-\r
- rot[2][2] = (float)1.0;\r
-\r
-\r
-\r
- preMultiply(rot);\r
-\r
+/**\r
+ * DOCUMENT ME!\r
+ *\r
+ * @author $author$\r
+ * @version $Revision$\r
+ */\r
+public class RotatableMatrix\r
+{\r
+ float[][] matrix;\r
+ float[] temp;\r
+ float[][] rot;\r
+\r
+ /**\r
+ * Creates a new RotatableMatrix object.\r
+ *\r
+ * @param rows DOCUMENT ME!\r
+ * @param cols DOCUMENT ME!\r
+ */\r
+ public RotatableMatrix(int rows, int cols)\r
+ {\r
+ matrix = new float[rows][cols];\r
+\r
+ temp = new float[3];\r
+\r
+ rot = new float[3][3];\r
}\r
\r
- if (axis == 'x') {\r
-\r
- rot[0][0] = (float)1.0;\r
-\r
- rot[0][1] = (float)0.0;\r
-\r
- rot[0][2] = (float)0.0;\r
-\r
-\r
-\r
- rot[1][0] = (float)0.0;\r
-\r
- rot[1][1] = (float)costheta;\r
-\r
- rot[1][2] = (float)sintheta;\r
-\r
-\r
-\r
- rot[2][0] = (float)0.0;\r
-\r
- rot[2][1] = (float)-sintheta;\r
-\r
- rot[2][2] = (float)costheta;\r
-\r
-\r
-\r
- preMultiply(rot);\r
-\r
-\r
-\r
+ /**\r
+ * DOCUMENT ME!\r
+ *\r
+ * @param i DOCUMENT ME!\r
+ * @param j DOCUMENT ME!\r
+ * @param value DOCUMENT ME!\r
+ */\r
+ public void addElement(int i, int j, float value)\r
+ {\r
+ matrix[i][j] = value;\r
}\r
\r
- if (axis == 'y') {\r
-\r
- rot[0][0] = (float)costheta;\r
-\r
- rot[0][1] = (float)0.0;\r
-\r
- rot[0][2] = (float)-sintheta;\r
-\r
-\r
-\r
- rot[1][0] = (float)0.0;\r
-\r
- rot[1][1] = (float)1.0;\r
-\r
- rot[1][2] = (float)0.0;\r
-\r
-\r
-\r
- rot[2][0] = (float)sintheta;\r
-\r
- rot[2][1] = (float)0.0;\r
-\r
- rot[2][2] = (float)costheta;\r
-\r
-\r
-\r
- preMultiply(rot);\r
-\r
+ /**\r
+ * DOCUMENT ME!\r
+ */\r
+ public void print()\r
+ {\r
+ System.out.println(matrix[0][0] + " " + matrix[0][1] + " " +\r
+ matrix[0][2]);\r
\r
+ System.out.println(matrix[1][0] + " " + matrix[1][1] + " " +\r
+ matrix[1][2]);\r
\r
+ System.out.println(matrix[2][0] + " " + matrix[2][1] + " " +\r
+ matrix[2][2]);\r
}\r
\r
+ /**\r
+ * DOCUMENT ME!\r
+ *\r
+ * @param degrees DOCUMENT ME!\r
+ * @param axis DOCUMENT ME!\r
+ */\r
+ public void rotate(float degrees, char axis)\r
+ {\r
+ float costheta = (float) Math.cos((degrees * Math.PI) / (float) 180.0);\r
\r
+ float sintheta = (float) Math.sin((degrees * Math.PI) / (float) 180.0);\r
\r
- }\r
-\r
+ if (axis == 'z')\r
+ {\r
+ rot[0][0] = (float) costheta;\r
\r
+ rot[0][1] = (float) -sintheta;\r
\r
- public float[] vectorMultiply(float[] vect) {\r
+ rot[0][2] = (float) 0.0;\r
\r
- temp[0] = vect[0];\r
+ rot[1][0] = (float) sintheta;\r
\r
- temp[1] = vect[1];\r
+ rot[1][1] = (float) costheta;\r
\r
- temp[2] = vect[2];\r
+ rot[1][2] = (float) 0.0;\r
\r
+ rot[2][0] = (float) 0.0;\r
\r
+ rot[2][1] = (float) 0.0;\r
\r
- for (int i = 0; i < 3; i++) {\r
+ rot[2][2] = (float) 1.0;\r
\r
- temp[i] = matrix[i][0]*vect[0] + matrix[i][1]*vect[1] + matrix[i][2]*vect[2];\r
-\r
- }\r
+ preMultiply(rot);\r
+ }\r
\r
+ if (axis == 'x')\r
+ {\r
+ rot[0][0] = (float) 1.0;\r
\r
+ rot[0][1] = (float) 0.0;\r
\r
- vect[0] = temp[0];\r
+ rot[0][2] = (float) 0.0;\r
\r
- vect[1] = temp[1];\r
+ rot[1][0] = (float) 0.0;\r
\r
- vect[2] = temp[2];\r
+ rot[1][1] = (float) costheta;\r
\r
+ rot[1][2] = (float) sintheta;\r
\r
+ rot[2][0] = (float) 0.0;\r
\r
- return vect;\r
+ rot[2][1] = (float) -sintheta;\r
\r
- }\r
+ rot[2][2] = (float) costheta;\r
\r
+ preMultiply(rot);\r
+ }\r
\r
+ if (axis == 'y')\r
+ {\r
+ rot[0][0] = (float) costheta;\r
\r
- public void preMultiply(float mat[][]) {\r
+ rot[0][1] = (float) 0.0;\r
\r
- float tmp[][] = new float[3][3];\r
+ rot[0][2] = (float) -sintheta;\r
\r
+ rot[1][0] = (float) 0.0;\r
\r
+ rot[1][1] = (float) 1.0;\r
\r
- for (int i = 0; i < 3 ; i++) {\r
+ rot[1][2] = (float) 0.0;\r
\r
- for (int j = 0; j < 3; j++ ) {\r
+ rot[2][0] = (float) sintheta;\r
\r
- tmp[i][j] = mat[i][0]*matrix[0][j] +\r
+ rot[2][1] = (float) 0.0;\r
\r
- mat[i][1]*matrix[1][j] +\r
-\r
- mat[i][2]*matrix[2][j];\r
-\r
- }\r
-\r
- }\r
-\r
-\r
-\r
- for (int i = 0; i < 3 ; i++) {\r
-\r
- for (int j = 0; j < 3; j++ ) {\r
-\r
- matrix[i][j] = tmp[i][j];\r
-\r
- }\r
+ rot[2][2] = (float) costheta;\r
\r
+ preMultiply(rot);\r
+ }\r
}\r
\r
- }\r
-\r
-\r
-\r
- public void postMultiply(float mat[][]) {\r
-\r
- float tmp[][] = new float[3][3];\r
+ /**\r
+ * DOCUMENT ME!\r
+ *\r
+ * @param vect DOCUMENT ME!\r
+ *\r
+ * @return DOCUMENT ME!\r
+ */\r
+ public float[] vectorMultiply(float[] vect)\r
+ {\r
+ temp[0] = vect[0];\r
\r
+ temp[1] = vect[1];\r
\r
+ temp[2] = vect[2];\r
\r
- for (int i = 0; i < 3 ; i++) {\r
+ for (int i = 0; i < 3; i++)\r
+ {\r
+ temp[i] = (matrix[i][0] * vect[0]) + (matrix[i][1] * vect[1]) +\r
+ (matrix[i][2] * vect[2]);\r
+ }\r
\r
- for (int j = 0; j < 3; j++ ) {\r
+ vect[0] = temp[0];\r
\r
- tmp[i][j] = matrix[i][0]*mat[0][j] +\r
+ vect[1] = temp[1];\r
\r
- matrix[i][1]*mat[1][j] +\r
-\r
- matrix[i][2]*mat[2][j];\r
-\r
- }\r
+ vect[2] = temp[2];\r
\r
+ return vect;\r
}\r
\r
-\r
-\r
- for (int i = 0; i < 3 ; i++) {\r
-\r
- for (int j = 0; j < 3; j++ ) {\r
-\r
- matrix[i][j] = tmp[i][j];\r
-\r
- }\r
-\r
+ /**\r
+ * DOCUMENT ME!\r
+ *\r
+ * @param mat DOCUMENT ME!\r
+ */\r
+ public void preMultiply(float[][] mat)\r
+ {\r
+ float[][] tmp = new float[3][3];\r
+\r
+ for (int i = 0; i < 3; i++)\r
+ {\r
+ for (int j = 0; j < 3; j++)\r
+ {\r
+ tmp[i][j] = (mat[i][0] * matrix[0][j]) +\r
+ (mat[i][1] * matrix[1][j]) + (mat[i][2] * matrix[2][j]);\r
+ }\r
+ }\r
+\r
+ for (int i = 0; i < 3; i++)\r
+ {\r
+ for (int j = 0; j < 3; j++)\r
+ {\r
+ matrix[i][j] = tmp[i][j];\r
+ }\r
+ }\r
}\r
\r
- }\r
-\r
-\r
-\r
- public static void main(String[] args) {\r
-\r
-\r
-\r
- RotatableMatrix m = new RotatableMatrix(3,3);\r
-\r
- m.addElement(0,0,1);\r
-\r
- m.addElement(0,1,0);\r
-\r
- m.addElement(0,2,0);\r
-\r
- m.addElement(1,0,0);\r
-\r
- m.addElement(1,1,2);\r
-\r
- m.addElement(1,2,0);\r
-\r
- m.addElement(2,0,0);\r
-\r
- m.addElement(2,1,0);\r
-\r
- m.addElement(2,2,1);\r
-\r
-\r
-\r
- m.print();\r
-\r
-\r
-\r
- RotatableMatrix n = new RotatableMatrix(3,3);\r
-\r
- n.addElement(0,0,2);\r
-\r
- n.addElement(0,1,1);\r
-\r
- n.addElement(0,2,1);\r
-\r
- n.addElement(1,0,2);\r
-\r
- n.addElement(1,1,1);\r
+ /**\r
+ * DOCUMENT ME!\r
+ *\r
+ * @param mat DOCUMENT ME!\r
+ */\r
+ public void postMultiply(float[][] mat)\r
+ {\r
+ float[][] tmp = new float[3][3];\r
+\r
+ for (int i = 0; i < 3; i++)\r
+ {\r
+ for (int j = 0; j < 3; j++)\r
+ {\r
+ tmp[i][j] = (matrix[i][0] * mat[0][j]) +\r
+ (matrix[i][1] * mat[1][j]) + (matrix[i][2] * mat[2][j]);\r
+ }\r
+ }\r
+\r
+ for (int i = 0; i < 3; i++)\r
+ {\r
+ for (int j = 0; j < 3; j++)\r
+ {\r
+ matrix[i][j] = tmp[i][j];\r
+ }\r
+ }\r
+ }\r
\r
- n.addElement(1,2,1);\r
+ /**\r
+ * DOCUMENT ME!\r
+ *\r
+ * @param args DOCUMENT ME!\r
+ */\r
+ public static void main(String[] args)\r
+ {\r
+ RotatableMatrix m = new RotatableMatrix(3, 3);\r
\r
- n.addElement(2,0,2);\r
+ m.addElement(0, 0, 1);\r
\r
- n.addElement(2,1,1);\r
+ m.addElement(0, 1, 0);\r
\r
- n.addElement(2,2,1);\r
+ m.addElement(0, 2, 0);\r
\r
+ m.addElement(1, 0, 0);\r
\r
+ m.addElement(1, 1, 2);\r
\r
- n.print();\r
+ m.addElement(1, 2, 0);\r
\r
+ m.addElement(2, 0, 0);\r
\r
+ m.addElement(2, 1, 0);\r
\r
- //m.postMultiply(n.matrix);\r
+ m.addElement(2, 2, 1);\r
\r
- //m.print();\r
+ m.print();\r
\r
- // m.rotate(45,'z',new RotatableMatrix(3,3));\r
+ RotatableMatrix n = new RotatableMatrix(3, 3);\r
\r
+ n.addElement(0, 0, 2);\r
\r
+ n.addElement(0, 1, 1);\r
\r
- float vect[] = new float[3];\r
+ n.addElement(0, 2, 1);\r
\r
- vect[0] = 2;\r
+ n.addElement(1, 0, 2);\r
\r
- vect[1] = 4;\r
+ n.addElement(1, 1, 1);\r
\r
- vect[2] = 6;\r
+ n.addElement(1, 2, 1);\r
\r
+ n.addElement(2, 0, 2);\r
\r
+ n.addElement(2, 1, 1);\r
\r
- vect = m.vectorMultiply(vect);\r
+ n.addElement(2, 2, 1);\r
\r
- System.out.println(vect[0] + " " + vect[1] + " " + vect[2]);\r
+ n.print();\r
\r
+ //m.postMultiply(n.matrix);\r
+ //m.print();\r
+ // m.rotate(45,'z',new RotatableMatrix(3,3));\r
+ float[] vect = new float[3];\r
\r
+ vect[0] = 2;\r
\r
- }\r
+ vect[1] = 4;\r
\r
- public void setIdentity() {\r
+ vect[2] = 6;\r
\r
- matrix[0][0] = (float)1.0;\r
+ vect = m.vectorMultiply(vect);\r
\r
- matrix[1][1] = (float)1.0;\r
+ System.out.println(vect[0] + " " + vect[1] + " " + vect[2]);\r
+ }\r
\r
- matrix[2][2] = (float)1.0;\r
+ /**\r
+ * DOCUMENT ME!\r
+ */\r
+ public void setIdentity()\r
+ {\r
+ matrix[0][0] = (float) 1.0;\r
\r
- matrix[0][1] = (float)0.0;\r
+ matrix[1][1] = (float) 1.0;\r
\r
- matrix[0][2] = (float)0.0;\r
+ matrix[2][2] = (float) 1.0;\r
\r
- matrix[1][0] = (float)0.0;\r
+ matrix[0][1] = (float) 0.0;\r
\r
- matrix[1][2] = (float) 0.0;\r
+ matrix[0][2] = (float) 0.0;\r
\r
- matrix[2][0] = (float)0.0;\r
+ matrix[1][0] = (float) 0.0;\r
\r
- matrix[2][1] = (float)0.0;\r
+ matrix[1][2] = (float) 0.0;\r
\r
- }\r
+ matrix[2][0] = (float) 0.0;\r
\r
+ matrix[2][1] = (float) 0.0;\r
+ }\r
}\r
-\r
-\r
-\r