[讨论]高斯消元法的程序对比

god

下面我收集的高斯消元法的程序。还是请大家看看，并提点意见。
3 z; \0 i  c& H5 K2 Q0 f

弘道

谢谢楼主……辛苦啦！………………

0.9清1.8清2.7清

wuyuwuyul 自己程序拿出来看看呗！干嘛还要体力A！

Rekcahpu

刚学高斯消元

deven1985

，不错，前几天的数值分析课上也是做的这个！

xiaozhangwh

谢谢分享，来试一下

hduscimc

#include<iostream>
5 l' O7 P2 I8 R- e#include<cmath>
using namespace std;
( C* d% I+ S3 A% T) xint main()
{
8 K* g: q; l: V% v9 p        int i;
6 ^" ?$ P* c. C8 w        float *x;
5 W3 {: z+ X6 L  n: h/ }5 i* @        float c[3][4] ={0.101,2.304,3.555,1.183,
7 u! b  m' J' {5 e0 @                                        -1.347,3.712,4.623,2.137,
6 P+ P& l: p& _( ?7 X                                        -2.835,1.072,5.643,3.035};
1 k9 ~5 {8 y; U: G        float *ColPivot(float *,int);
        x=ColPivot(c[0],3);
3 }- p. y5 b) I- P. b5 {: _* O2 ?2 ?* }        for( i=0;i<=2;i++ )
5 C5 v: F" x1 c1 |                cout<<"x("<<i<<")="<<x[i]<<endl;
}

float *ColPivot( float *c,int n )
2 a& t; S3 P% U9 I7 U{
( t/ y1 F4 s" n2 R& X" n% T! e  ^        int i,j,t,k;
        float *x,p;
        x=new float[n*sizeof(float)];
) r) u- w/ D0 a        for( i=0;i<=n-2;i++)
8 }' ^3 O/ o2 l3 l. W9 [        {
                k=i;
3 Q, I. I3 b- d" D) S. V2 b                for(j=i+1;j<=n-1;j++)
                        if(fabs(*(c+j*(n+1)+i))>(fabs(*(c+k*(n+1)+i))))
- t. L) E% `/ a                                k=j;
                        if(k!=i)
+ F1 }" @0 a, T; x# V                                for( j=i;j<=n;j++ )
                                {
0 v! `1 F2 z# _  [* g                                        p=*(c+i*(n+1)+j);
                                        *(c+i*(n+1)+j)=*(c+k*(n+1)+j);
* o( ^1 s5 |' y+ e" s  M' g# {                                        *(c+k*(n+1)+j)=p;
; A! U3 L; J! H! E$ Z6 l5 T                                }
3 J1 i6 o" d3 D- l                for( j=i+1;j<=n-1;j++ )
( v; P" H  ^; I# w* N0 K4 n$ m                {
                        p=(*(c+j*(n+1)+i))/(*(c+i*(n+1)+i));
                        for( t=i;t<=n;t++ )
9 q. l1 N7 C7 H: Y                                *(c+j*(n+1)+t)-=p*(*(c+i*(n+1)+t));
                }
! J% O0 _; s& c  l$ c5 B$ S        }
- |3 _4 g: Y+ ^3 C        for( i=n-1;i>=0;i--)
        {
$ s5 M8 N. K# e- k4 G7 N% e                for( j=n-1;j>=i+1;j--)
1 T, h6 a( V) v7 s; k                        (*(c+i*(n+1)+n))-=x[j]*(*(c+i*(n+1)+j));
                x[i]=*(c+i*(n+1)+n)/(*(c+i*(n+1)+i));
+ P0 q# b/ E. [2 X+ K        }
        return x;
+ Q% z+ Q. \) M! t4 y# m}

shuainaruto

感谢分享，下来研究一下