I defined a function in C as below:
int *face_2_cells(int boundary_field_id) { int b_cells_field[4]; if(boundary_field_id==1) { b_cells_field[0]=1; b_cell