Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025 #include "plplotP.h"
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052 void
00053 plP_affine_identity( PLFLT *affine_vector )
00054 {
00055 affine_vector[0] = 1.;
00056 affine_vector[1] = 0.;
00057 affine_vector[2] = 0.;
00058 affine_vector[3] = 1.;
00059 affine_vector[4] = 0.;
00060 affine_vector[5] = 0.;
00061 }
00062
00063
00064
00065 void
00066 plP_affine_translate( PLFLT *affine_vector, PLFLT xtranslate, PLFLT ytranslate )
00067 {
00068 affine_vector[0] = 1.;
00069 affine_vector[1] = 0.;
00070 affine_vector[2] = 0.;
00071 affine_vector[3] = 1.;
00072
00073
00074
00075 affine_vector[4] = -xtranslate;
00076 affine_vector[5] = -ytranslate;
00077 }
00078
00079
00080
00081 void
00082 plP_affine_scale( PLFLT *affine_vector, PLFLT xscale, PLFLT yscale )
00083 {
00084
00085
00086
00087 if ( xscale == 0. )
00088 {
00089 plwarn( "plP_affine_scale: attempt to scale X coordinates by zero." );
00090 xscale = 1.;
00091 }
00092 if ( yscale == 0. )
00093 {
00094 plwarn( "plP_affine_scale: attempt to scale Y coordinates by zero." );
00095 yscale = 1.;
00096 }
00097 affine_vector[0] = 1. / xscale;
00098 affine_vector[1] = 0.;
00099 affine_vector[2] = 0.;
00100 affine_vector[3] = 1. / yscale;
00101 affine_vector[4] = 0.;
00102 affine_vector[5] = 0.;
00103 }
00104
00105
00106
00107
00108 void
00109 plP_affine_rotate( PLFLT *affine_vector, PLFLT angle )
00110 {
00111 PLFLT cosangle = cos( PI * angle / 180. );
00112 PLFLT sinangle = sin( PI * angle / 180. );
00113 affine_vector[0] = cosangle;
00114 affine_vector[1] = -sinangle;
00115 affine_vector[2] = sinangle;
00116 affine_vector[3] = cosangle;
00117 affine_vector[4] = 0.;
00118 affine_vector[5] = 0.;
00119 }
00120
00121
00122
00123
00124 void
00125 plP_affine_xskew( PLFLT *affine_vector, PLFLT angle )
00126 {
00127 PLFLT tanangle = tan( PI * angle / 180. );
00128 affine_vector[0] = 1.;
00129 affine_vector[1] = 0.;
00130 affine_vector[2] = -tanangle;
00131 affine_vector[3] = 1.;
00132 affine_vector[4] = 0.;
00133 affine_vector[5] = 0.;
00134 }
00135
00136
00137
00138
00139 void
00140 plP_affine_yskew( PLFLT *affine_vector, PLFLT angle )
00141 {
00142 PLFLT tanangle = tan( PI * angle / 180. );
00143 affine_vector[0] = 1.;
00144 affine_vector[1] = -tanangle;
00145 affine_vector[2] = 0.;
00146 affine_vector[3] = 1.;
00147 affine_vector[4] = 0.;
00148 affine_vector[5] = 0.;
00149 }
00150
00151
00152
00153
00154
00155
00156
00157 void
00158 plP_affine_multiply(
00159 PLFLT *affine_vectorA,
00160 const PLFLT *affine_vectorB,
00161 const PLFLT *affine_vectorC )
00162 {
00163 int i;
00164 PLFLT result[NAFFINE];
00165
00166 result[0] = affine_vectorB[0] * affine_vectorC[0] +
00167 affine_vectorB[2] * affine_vectorC[1];
00168 result[2] = affine_vectorB[0] * affine_vectorC[2] +
00169 affine_vectorB[2] * affine_vectorC[3];
00170 result[4] = affine_vectorB[0] * affine_vectorC[4] +
00171 affine_vectorB[2] * affine_vectorC[5] +
00172 affine_vectorB[4];
00173
00174 result[1] = affine_vectorB[1] * affine_vectorC[0] +
00175 affine_vectorB[3] * affine_vectorC[1];
00176 result[3] = affine_vectorB[1] * affine_vectorC[2] +
00177 affine_vectorB[3] * affine_vectorC[3];
00178 result[5] = affine_vectorB[1] * affine_vectorC[4] +
00179 affine_vectorB[3] * affine_vectorC[5] +
00180 affine_vectorB[5];
00181
00182 for ( i = 0; i < NAFFINE; i++ )
00183 affine_vectorA[i] = result[i];
00184 }