Ejemplo n.º 1
0
int main(int argc, char **argv){
  static xyc *Z;
  static nde *N;
  FILE *fp;
  int i;

  initop(argc, argv);
  fp = argf(argc, argv);

  fp2msh(fp,&Z, &N);

  pltmsh(stdout,Z,N);
}
Ejemplo n.º 2
0
main(int argc, char **argv){
  int i, n = 512;
  static MX *A; static double *x, *b;
  initop(argc,argv);
  initmx(A,n+1,8); ary1(x, n+1); ary1(b, n+1);

  for ( i =1; i<=n; i++){
    mx(A,i,i) = 2.0; b[i] = 1.0;
  }
  for (i=1;i<n;i++) mx(A,i,i+1) = -1.0;
  for (i=1;i<n;i++) mx(A,i+1,i) = -1.0;

  solver(A,x,b);

  for (i=1;i<=n;i++) printf("%f\n", x[i] );
}
Ejemplo n.º 3
0
int
main(int argc, char **argv)
{
  static xyc *Z; static nde *N;
  static double **A, *u;

  initop(argc, argv);
  fp2mesh(stdfp(),Z,N);

  ary2(A,dim1(Z)+1, dim1(Z)+1); ary1(u,dim1(Z)+1);

  set_A(Z,N,A); set_u(Z,u);

  esolver(A,u);

  plt(NULL,NULL,Z,N,u); sleep(1000);
  return 0;
}