Ejemplo n.º 1
0
void eval_At_times_u(int N, const double u[], double Au[])
{
    int i,j;
    for(i=0; i<N; i++)
    {
        Au[i]=0;
        for(j=0; j<N; j++) Au[i]+=eval_A(j,i)*u[j];
    }
}
Ejemplo n.º 2
0
int prepare_A(int N)
{
  int i,j;

  N_global=N;
  A_global=(double*)malloc(N*N*sizeof(double));

  if(A_global==NULL) return -1;

  for(i=0;i<N;i++)
    {
      for(j=0;j<N;j++)
	{
	  A_global[i*N+j]=eval_A(i,j);
	}
    }
  return 0;
}