void eval_AtA_times_u(int N, const double u[], double AtAu[]) { double *v = malloc(N * sizeof(double)); eval_A_times_u(N,u,v); eval_At_times_u(N,v,AtAu); free(v); }
void eval_AtA_times_u(int N, const double u[], double AtAu[]) { double v[N]; eval_A_times_u(N,u,v); eval_At_times_u(N,v,AtAu); }