109 {
111 int jp1;
115
116
117 if( m1>m2 || n1>n2 )
118 {
119 return;
120 }
121
122
123
124
125 if( isforward )
126 {
127 if( n1!=n2 )
128 {
129
130
131
132
133 for(
j=m1;
j<=m2-1;
j++)
134 {
137 if( ctemp!=1 || stemp!=0 )
138 {
145 }
146 }
147 }
148 else
149 {
150
151
152
153
154 for(
j=m1;
j<=m2-1;
j++)
155 {
158 if( ctemp!=1 || stemp!=0 )
159 {
161 a(
j+1,n1) = ctemp*temp-stemp*a(
j,n1);
162 a(
j,n1) = stemp*temp+ctemp*a(
j,n1);
163 }
164 }
165 }
166 }
167 else
168 {
169 if( n1!=n2 )
170 {
171
172
173
174
175 for(
j=m2-1;
j>=m1;
j--)
176 {
179 if( ctemp!=1 || stemp!=0 )
180 {
187 }
188 }
189 }
190 else
191 {
192
193
194
195
196 for(
j=m2-1;
j>=m1;
j--)
197 {
200 if( ctemp!=1 || stemp!=0 )
201 {
203 a(
j+1,n1) = ctemp*temp-stemp*a(
j,n1);
204 a(
j,n1) = stemp*temp+ctemp*a(
j,n1);
205 }
206 }
207 }
208 }
209 }
raw_vector< T > getvector(int iStart, int iEnd)
raw_vector< T > getrow(int iRow, int iColumnStart, int iColumnEnd)
const CanonicalForm int s
void vmove(raw_vector< T > vdst, const_raw_vector< T > vsrc)
void vadd(raw_vector< T > vdst, const_raw_vector< T > vsrc)
void vmul(raw_vector< T > vdst, T2 alpha)
void vsub(raw_vector< T > vdst, const_raw_vector< T > vsrc)