#ifndef __EVAL_H__
#define __EVAL_H__

#include "indices.h"
#include <petscksp.h>

void eval_a1(int ele_index, double x, double y);/*}}}*/

void eval_b1(int ele_index, double x, double y);/*}}}*/

void eval_c1(int ele_index, double x, double y);/*}}}*/

PetscScalar eval_FE_fct(PetscScalar* coeff_vec,
                 PetscScalar (**FE_space)(PetscScalar, PetscScalar),
         PetscInt ele_index,
         PetscInt elems_x,
         PetscScalar x1,
         PetscScalar x2,
         grid_type grid);/*}}}*/


PetscScalar eval_glob_basis(PetscScalar (**FE_space)(PetscScalar, PetscScalar),
                     PetscInt global_index,
             PetscInt ele_index,
             PetscInt elems_x,
             PetscScalar x1,
             PetscScalar x2,
             grid_type grid);/*}}}*/

 #endif