/** @file */
#ifndef __RHS_H__
#define __RHS_H__
#include <petscksp.h>

PetscErrorCode calculate_rhs(Vec rhs, PetscInt rhs_length, PetscInt elems_x,
                             PetscInt elems_y, PetscScalar* u1_vec,
			     PetscScalar* u2_vec, PetscScalar* p_vec,
			     PetscScalar k, PetscScalar nu,
			     PetscInt p_vec_length, PetscInt vel_y_vec_length);


PetscScalar evaluate_rhs1(PetscInt glob_index,
                   PetscInt ele_index,
                   PetscScalar x1,
		   PetscScalar x2,
		   PetscInt elems_x,
		   PetscScalar k,
		   PetscScalar nu,
		   PetscScalar* u1_vec,
		   PetscScalar* u2_vec,
		   PetscScalar* p_vec);


PetscScalar evaluate_rhs2(PetscInt glob_index,
                   PetscInt ele_index,
                   PetscScalar x1,
		   PetscScalar x2,
		   PetscInt elems_x,
		   PetscScalar k,
		   PetscScalar nu,
		   PetscScalar* u1_vec,
		   PetscScalar* u2_vec,
		   PetscScalar* p_vec);

#endif