Example #1
0
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);
}
Example #2
0
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); }