【C语言发布】非齐次线性方程组的解……（上）

RyanBern

前排。坐等更新。

另外这代码还得注意程序规范的问题，比如声明的时候通常不加变量名字，只说类型。
还有这个写法到底是啥意思

C 代码复制

scanf("%.0f\t", &[i]);

scanf("%.0f\t", &[i]);

RyanBern

寒冷魔王 发表于 2016-9-17 01:17
哇，小叶子也做这个了呢~
前排支持！

首先，是因为是叶子所以回复这么认真么/w\

如果要用C++(兼容C89但不兼容C99)编译呢，最好写成...

其实我早就跟叶子说了我不建议这样，但是谁让它是一个 C 语言程序并且自动打开 gcc 的 c99 标准呢……似乎我好像永远活在了 ANSI C 标准中无法自拔。

另外“使用C++编译”似乎是个不准确的说法呢……

用struct Matrix代替double[]的好处是不会和其他的double[]弄混，而且能够直接复制struct Matrix哦~

你确定？除非你显式写出复制的函数，否则是绝对不能这样“直接”复制的。

如果单纯考虑不弄混的话，可以直接定义一个类型，这样就不容易混了。

C 代码复制

typedef double* MatElePr;

typedef double* MatElePr;

另外从你的回帖里我能看出你尽量把程序做得面向对象一点，于是想方设法去写 Constructor Destructor 等东西。其实我在群里是这样教叶子的：不许用 struct，不许有面向对象的思维，给我用纯指针弄。

看来我的反面向对象思维越来越严重了。不过，有关的高性能的数学库确实不会考虑这些封装呢。

RyanBern

寒冷魔王 发表于 2016-9-18 00:50
/w\感谢您的回复，看到小叶子这么认真所以激动的有些忘乎所以了呢~

关于

C 代码复制

struct A1{
  double data[5];
};

struct A1{


  double data[5];


};

C 代码复制

struct A2{
  double *data;
};

struct A2{


  double *data;


};

下面两个有本质的区别啊，不要把它们当成一样的了喂！

C 代码复制

printf("%d %d\n", sizeof(struct A1), sizeof(struct A2));

printf("%d %d\n", sizeof(struct A1), sizeof(struct A2));

这样会看到很明显的不同。

概况来说，如果是动态分配的内存，是绝对不可以这样玩的。

C 代码复制

struct A{
  double *data;
};
 
int main(int argc, char **argv){
  struct A a1, a2;
  /* suppose we call 'constructors' for a1 and a2 */
  init(&a1, n); init(&a2, n);
  a2 = a1; /* the problem occurs ! */
  a1.data[0] = 1; /* then what is a2.data[0]? */
  /* suppose we call 'destructors' for a1 and a2 */
  finalize(&a1); finalize(&a2); /* core dumped !*/
  return 0;
}

struct A{


  double *data;


};


 


int main(int argc, char **argv){


  struct A a1, a2;


  /* suppose we call 'constructors' for a1 and a2 */


  init(&a1, n); init(&a2, n);


  a2 = a1; /* the problem occurs ! */


  a1.data[0] = 1; /* then what is a2.data[0]? */


  /* suppose we call 'destructors' for a1 and a2 */


  finalize(&a1); finalize(&a2); /* core dumped !*/


  return 0;


}

写得有点仓促，但是我觉得意思应该表示清楚了。

RyanBern

糊了一个 PLU，里面只有 PLU 分解，没有解方程组的部分。

叶子看完了可以把解方程组的部分补上。不知道有没有 BUG，没仔细看过。

源码：
https://github.com/RyanBernX/math/tree/master/lu
说明书：
https://ryanbernx.github.io/math/