1 /* -*- c -*- */ 2 3 /* 4 * vector.h 5 * 6 * metapixel 7 * 8 * Copyright (C) 1997-1999 Mark Probst 9 * 10 * This program is free software; you can redistribute it and/or 11 * modify it under the terms of the GNU General Public License 12 * as published by the Free Software Foundation; either version 2 13 * of the License, or (at your option) any later version. 14 * 15 * This program is distributed in the hope that it will be useful, 16 * but WITHOUT ANY WARRANTY; without even the implied warranty of 17 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 18 * GNU General Public License for more details. 19 * 20 * You should have received a copy of the GNU General Public License 21 * along with this program; if not, write to the Free Software 22 * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. 23 */ 24 25 typedef struct _Vector2D 26 { 27 double x; 28 double y; 29 } Vector2D; 30 31 typedef struct _Vector3D 32 { 33 double x; 34 double y; 35 double z; 36 } Vector3D; 37 38 typedef double Vector4D[4]; 39 40 typedef double Matrix3D[3][3]; 41 typedef double Matrix4D[4][4]; 42 typedef double Matrix4x3[4][3]; 43 44 Vector2D MakeVector2D (double x, double y); 45 Vector2D* InitVector2D (Vector2D *vec, double x, double y); 46 Vector2D* CopyVector2D (Vector2D *dest, const Vector2D *src); 47 double Abs2D (const Vector2D *vec); 48 Vector2D* AddVectors2D (Vector2D *dest, const Vector2D *vec1, 49 const Vector2D *vec2); 50 Vector2D* SubVectors2D (Vector2D *dest, const Vector2D *vec1, 51 const Vector2D *vec2); 52 Vector2D* MultScalar2D (Vector2D *dest, double d, const Vector2D *vec); 53 Vector2D* MultVectors2D (Vector2D *vec, const Vector2D *vec1, 54 const Vector2D *vec2); 55 Vector2D* Unity2D (Vector2D *dest, Vector2D *vec); 56 Vector2D* Rectangular2DToPolar (Vector2D *dest, const Vector2D *src); 57 Vector2D* Polar2DToRectangular (Vector2D *dest, const Vector2D *src); 58 59 Vector3D MakeVector3D (double x, double y, double z); 60 Vector3D* InitVector3D (Vector3D *vec, double x, double y, double z); 61 Vector3D* CopyVector3D (Vector3D *dest, const Vector3D *src); 62 double Abs3D (const Vector3D *vec); 63 Vector3D* AddVectors3D (Vector3D *dest, const Vector3D *vec1, 64 const Vector3D *vec2); 65 Vector3D* SubVectors3D (Vector3D *dest, const Vector3D *vec1, 66 const Vector3D *vec2); 67 Vector3D* MultScalar3D (Vector3D *dest, double d, const Vector3D *vec); 68 Vector3D* MultVectors3D (Vector3D *vec, const Vector3D *vec1, 69 const Vector3D *vec2); 70 double DotProduct3D (const Vector3D *vec1, const Vector3D *vec2); 71 Vector3D* CrossProduct3D (Vector3D *vec, const Vector3D *vec1, 72 const Vector3D *vec2); 73 Vector3D* Unity3D (Vector3D *dest, const Vector3D *vec); 74 75 Matrix3D* InitMatrix3D (Matrix3D *mat); 76 Matrix3D* CopyMatrix3D (Matrix3D *dest, const Matrix3D *src); 77 Vector3D* MultMatrixVector3D (Vector3D *dest, const Matrix3D *mat, const Vector3D *vec); 78 79 Matrix4D* InitMatrix4D (Matrix4D *mat); 80 Matrix4D* CopyMatrix4D (Matrix4D *dest, const Matrix4D *src); 81 Matrix4D* MultMatrix4D (Matrix4D *dest, 82 const Matrix4D *mat1, const Matrix4D *mat2); 83 Vector3D* ApplyTransformation (Vector3D *dest, 84 const Matrix4D *mat, const Vector3D *vec); 85 Matrix4D* InitTranslationMatrix (Matrix4D *mat, double x, double y, double z); 86 Matrix4D* InitXRotationMatrix (Matrix4D *mat, double theta); 87 Matrix4D* InitYRotationMatrix (Matrix4D *mat, double theta); 88 Matrix4D* InitZRotationMatrix (Matrix4D *mat, double theta); 89 90 Matrix4x3* InitMatrix4x3 (Matrix4x3 *mat); 91 Matrix4x3* MakeMatrix4x3 (Matrix4x3 *dest, 92 const Vector3D *vec1, const Vector3D *vec2, 93 const Vector3D *vec3, const Vector3D *vec4); 94 Matrix4x3* CopyMatrix4x3 (Matrix4x3 *dest, const Matrix4x3 *src); 95 Matrix4x3* MultMatrix4x3 (Matrix4x3 *dest, 96 const Matrix4D *mat1, const Matrix4x3 *mat2); 97 Vector3D* MultVector4DMatrix4x3 (Vector3D *dest, 98 const Vector4D *vec, const Matrix4x3 *mat); 99 100 Vector3D* CatmullRom (Vector3D *dest, 101 const Vector3D *vec1, const Vector3D *vec2, 102 const Vector3D *vec3, const Vector3D *vec4, 103 double t); 104 Matrix4x3* PrecompileCatmullRom (Matrix4x3 *dest, 105 const Vector3D *vec1, const Vector3D *vec2, 106 const Vector3D *vec3, const Vector3D *vec4); 107 Vector3D* QuickCatmullRom (Vector3D *dest, Matrix4x3 *matC, double t); 108