/*列主元消元法
% \5 r3 F% Y# F. S5 c 列主元消元法是在Gauss消元法的基础上的一个改进。, L2 z1 f) x0 H( E' E9 r
在Gauss消元法中,有些方程组虽然能够求解出结果,但是很难保证计算结果的可靠性& `2 Y$ J ^ P F( D) O' p# W4 ?
列主元消元法对此进行了改进,减少了在消元过程中的误差。, C6 s9 N& ~) H5 b
*/' s- I$ \6 @( w' O/ m6 Y- f2 J
/*函数名称:row_element_elimination_calculate;列主元消元法, I8 V( s A- T* G0 e( G
函数参数:int (*p)[3];线性方程组的系数行列式( O- @2 `* _' ~) p F6 l
int *B;线性方程组的右边常数向量
2 w, K0 R& p- B A" `1 Q) @ int size;线性方程组的阶数
3 S* H2 e, r% S. l9 ~, K 函数返回值:无,函数过程中输入线性方程组的求解结果 i. Q8 [1 F' ^8 o7 `+ ^, n6 ?
*/# G$ G( C2 q+ m
void row_element_elimination_calculate(double (*p)[3], double* B, int size) }: `; V% u( ^) ]7 e
{
! S7 ` Q! m, [( U+ H5 f0 @ double* X = new double[size];6 J7 b& c: w* H9 [
int i = 0, j = 0; V# f) M: w2 L* e# r
double largest = -1e10; //最大列主元
% z2 x L3 P( i6 t% Z& v4 Z int swap_index = -1;
( n' X3 s% N; |; Y/ Z' l7 |& \$ w for(i = 0; i < size; i++)
- P$ G6 w3 M! Y& g7 Z2 ~. Z {
0 z; \: l# y9 c1 Y; F* ]0 ~! w* T largest = -1e10; |