Go to the documentation of this file.00001
00002
00003
00004 #include "arith.h"
00005
00006 srot(n,sx,isx,jsx,incx,sy,isy,jsy,incy,c,s)
00007 MATRIX sx,sy;
00008 REAL c,s;
00009 int n,incx,incy,isx,jsx,isy,jsy;
00010 {
00011 REAL stemp;
00012 int i,ix,iy;
00013
00014 if(n <= 0)
00015 return;
00016
00017 if((incx != 1) || (incy != 1)){
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036 return;
00037 }
00038
00039
00040
00041 for(i=0; i<n; i++){
00042 stemp= c*sx[isx+i][jsx]+s*sy[isy+i][jsy];
00043 sy[isy+i][jsy]= c*sy[isy+i][jsy]-s*sx[isx+i][jsx];
00044 sx[isx+i][jsx]=stemp;
00045 }
00046
00047 return;
00048 }
00049
00050
00051