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]; } }
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; }