|
utility
|
|
cmd_head.c,
count_tokens.c,
cp_file.c,
errmsg.c,
file_opened.c,
find_line.c,
find_string.c,
fn_base.c,
fn_dir.c,
fn_ext.c,
fn_root.c,
gchoice.c,
gdata_line.c,
gdouble.c,
gfield.c,
gint.c,
gntoken.c,
graphics_ok.c,
gstring.c,
lock.c,
lower.c,
mat_multiply.c,
mat_transform.c,
mat_transpose.c,
mv_file.c,
new_extension.c,
nolead.c,
notrail.c,
ordinal.c,
pfield.c,
squish_blanks.c,
test.c,
unlock.c,
upper.c,
utility.h,
utility.make,
valchars.c,
vec_angle.c,
vec_crossprod.c,
vec_diff.c,
vec_dotprod.c,
vec_length.c,
vec_norm.c,
|
|
|
#include "utility.h"
void mat_multiply(matrix1,matrix2,matrix3)
/*============================================================================*/
/* PURPOSE: MULTIPLY TWO THREE BY THREE MATRIXES.
/* INPUTS:
/* MATRIX1 FIRST THREE BY THREE MATRIX.
/* MATRIX2 SECOND THREE BY THREE MATRIX.
/* OUTPUTS:
/* MATRIX3 PRODUCT OF MATRIX1 AND MATRIX2.
/* WRITTEN: M.V.GRIESHABER
/* LAST MODIFICATION: 11 JUNE 1991 MVG
*/
double matrix1[3][3];
double matrix2[3][3];
double matrix3[3][3];
{
int i;
int j;
for (i=0; i<3; i++)
{
for (j=0; j<3; j++)
{
matrix3[i][j]=matrix1[i][0]*matrix2[0][j] + matrix1[i][1]*matrix2[1][j] + matrix1[i][2]*matrix2[2][j];
}
}
return;
}
|
