***********************************************
integration_loadBalanceModule_in_LibMultiScale
***********************************************


shared variables LibMultiScale
*******************************

GLOBAL VARIABLE DEFINITION ARE IN
 ~/Bureau/src_libmultiscale/simulpack/mechanics/data.h

extern double * coordinates;		/* nodal coordinates			(x)		*/ mesh->nodes  with adaptation
extern int * connectivity;		/* element connectivity			(ix)		*/
extern int * element_type;		/* element type				(ie)		*/ depends on mesh->elm_size
extern int nodes;				/* current number of nodes		(numnp)	*/ mesh->nb_nodes
extern int nodes_element;		/* number of nodes/element		(nen)	*/ mesh->elm_size
extern int elements;			/* current number of elements	(numel)	*/ mesh->nb_cells
extern int dof_node;			/* number of dof per node		(ndf)	*/
extern int spatial_dimension;	/* spatial dimension			(ndm)	*/


coordinates = malloc(max_node_label*spatial_dimension*sizeof(double));
connectivity =  malloc(elements*nodes_element*sizeof(int));
element_type = malloc(elements*sizeof(int));

DEBUG 

printf("\n***************\nmax_node_label = %d\n***************\n", max_node_label);

printf("\n*******************************************************************\n");
printf("adlib_myrank = %d\n", adlib_myrank);
printf("nodes = %d\n", nodes);
printf("nodes_element = %d\n", nodes_element);
printf("elements = %d\n", elements);
printf("dof_node = %d\n", dof_node);
printf("spatial_dimension = %d\n", spatial_dimension);
int i = 0;
printf("coordinates = [ ");  for(i = 0 ; i<nodes*spatial_dimension ; i++) printf("%.2f ", coordinates[i]); printf("]\n");
printf("connectivity = [ "); for(i = 0 ; i<elements*nodes_element ; i++) printf("%d ", connectivity[i]); printf("]\n");
printf("element_type = [ "); for(i = 0 ; i<elements ; i++) printf("%d ", element_type[i]); printf("]\n");
printf("\n*******************************************************************\n");


*******************

la variable globale (mais locale à partitionate.c) contien les informations sur les résultats de partitionnement

____________________________________************************************************____________________________________


integration LibMultiScale
**************************

/home/elafrit/Bureau/src_libmultiscale/src/amel/AMEL.cpp
/home/elafrit/Bureau/src_libmultiscale/src/parser/parser.h
/home/elafrit/Bureau/src_libmultiscale/src/continuum/simulpack/material_simulpack.cpp
/home/elafrit/Bureau/src_libmultiscale/src/continuum/simulpack/simulpack_cfunctions.c
/home/elafrit/Bureau/src_libmultiscale/simulpack/mesh_io/mesh_io.h
/home/elafrit/Bureau/src_libmultiscale/simulpack/mesh_io/unv_io.c
/home/elafrit/Bureau/src_libmultiscale/simulpack/pmechanics/partitionate.h
/home/elafrit/Bureau/src_libmultiscale/simulpack/pmechanics/partitionate.c
/home/elafrit/Bureau/src_libmultiscale/simulpack/mechanics/data.h

grep -nir "load_mesh_UNV" .

printf("\n*************\n******OK*****\n*************\n");

kdevelop /home/elafrit/Bureau/src_libmultiscale/src/amel/AMEL.cpp /home/elafrit/Bureau/src_libmultiscale/src/parser/parser.h /home/elafrit/Bureau/src_libmultiscale/src/continuum/simulpack/material_simulpack.cpp /home/elafrit/Bureau/src_libmultiscale/src/continuum/simulpack/simulpack_cfunctions.c /home/elafrit/Bureau/src_libmultiscale/simulpack/mesh_io/mesh_io.h /home/elafrit/Bureau/src_libmultiscale/simulpack/mesh_io/unv_io.c /home/elafrit/Bureau/src_libmultiscale/simulpack/mesh_io/msh_io.c /home/elafrit/Bureau/src_libmultiscale/simulpack/pmechanics/partitionate.h /home/elafrit/Bureau/src_libmultiscale/simulpack/pmechanics/partitionate.c /home/elafrit/Bureau/src_libmultiscale/simulpack/mechanics/data.h &

les fichier cmake modifié pour la compilation :

kdevelop /home/elafrit/Bureau/src_libmultiscale/simulpack/mesh_io/CMakeLists.txt /home/elafrit/Bureau/src_libmultiscale/simulpack/CMakeLists.txt /home/elafrit/Bureau/src_libmultiscale/src/amel/CMakeLists.txt /home/elafrit/Bureau/src_libmultiscale/src/cmake/FindHGGEN.cmake &

************************************************************

### mettre à jour le fichier /home/elafrit/Bureau/src_libmultiscale/simulpack/mesh_io/CMakeLists.txt pour prendre en compte le nouveau fichier msh_io.c

SET(libmesh_ioSources
unv_io.c
msh_io.c
)

### mettre à jour le fichier /home/elafrit/Bureau/src_libmultiscale/simulpack/CMakeLists.txt pour prendre en compte de hggen

OPTION(USING_HGGEN "use hggen" OFF)
IF (USING_HGGEN)
FIND_PATH(HGGEN_INCLUDE_DIR hggen.h
	/usr/include
	$ENV{HGGEN_TOP}/include
	${HGGEN_DIR}/include
	${HGGEN_DIR}/src
)

FIND_LIBRARY(HGGEN_LIBRARY hggen
	/usr/lib
	$ENV{HGGEN_TOP}/lib
	${HGGEN_DIR}/lib
	${HGGEN_DIR}/src
)

IF(HGGEN_LIBRARY AND HGGEN_INCLUDE_DIR)
	SET(HGGEN_FOUND TRUE)
ELSE(HGGEN_LIBRARY AND HGGEN_INCLUDE_DIR)
	SET(HGGEN_FOUND FALSE)
	MESSAGE("HGGEN not found, set HGGEN_DIR.")
	SET(HGGEN_DIR "" CACHE PATH "Root path to hggen library")
ENDIF(HGGEN_LIBRARY AND HGGEN_INCLUDE_DIR)
ENDIF (USING_HGGEN)


### mettre à jour le fichier /home/elafrit/Bureau/src_libmultiscale/src/amel/CMakeLists.txt  pour prendre en compte de hggen

IF(USING_HGGEN)
TARGET_LINK_LIBRARIES(AMEL ${HGGEN_LIBRARY})
TARGET_LINK_LIBRARIES(LMPOST ${HGGEN_LIBRARY})
ENDIF(USING_HGGEN)

### ajouter le fichier /home/elafrit/Bureau/src_libmultiscale/src/cmake/FindHGGEN.cmake pour prendre en compte de hggen  et qui contient 

OPTION(USING_HGGEN "use hggen" OFF)
IF (USING_HGGEN)
FIND_PATH(HGGEN_INCLUDE_DIR hggen.h
	/usr/include
	$ENV{HGGEN_TOP}/include
	${HGGEN_DIR}/include
	${HGGEN_DIR}/src
)

FIND_LIBRARY(HGGEN_LIBRARY hggen
	/usr/lib
	$ENV{HGGEN_TOP}/lib
	${HGGEN_DIR}/lib
	${HGGEN_DIR}/src
)

IF(HGGEN_LIBRARY AND HGGEN_INCLUDE_DIR)
	SET(HGGEN_FOUND TRUE)
ELSE(HGGEN_LIBRARY AND HGGEN_INCLUDE_DIR)
	SET(HGGEN_FOUND FALSE)
	MESSAGE("HGGEN not found, set HGGEN_DIR.")
	SET(HGGEN_DIR "" CACHE PATH "Root path to hggen library")
ENDIF(HGGEN_LIBRARY AND HGGEN_INCLUDE_DIR)
ENDIF (USING_HGGEN)

### ajout de struct.h dans /home/elafrit/Bureau/src_libmultiscale/simulpack/

qui contient la définitions des structures 

### dans  /home/elafrit/Bureau/src_libmultiscale/src/continuum/simulpack/simulpack_cfunctions.c

mettre à jour la fonction :
void simulpack_init(double ts)
avec :
...

int NEW_STRATEGY = 1;
int OLD_STRATEGY = 0;
int strategy = NEW_STRATEGY;
if (strategy == NEW_STRATEGY) {
  load_mesh_MSH(fname);
  partition_mesh_hg(adlib_world_size);
}
else if (strategy == OLD_STRATEGY){
  load_mesh_UNV(fname);
  partition_mesh(adlib_world_size);
}

...

### dans /home/elafrit/Bureau/src_libmultiscale/simulpack/mesh_io/mesh_io.h

ajouter la declaration:
extern void load_mesh_MSH(const char * inputfile);


### dans /home/elafrit/Bureau/src_libmultiscale/simulpack/mesh_io/unv_io.c

ajouter la fonction:
void load_mesh_MSH(const char * inputfile){

	  /* load mesh */
	Mesh mesh;
	loadMSHMesh(inputfile, CELL_FIELD_MATERIAL_ID, &mesh);  // add CELL_MATERIAL variable

	  /* fil in global variabals */
	nodes_element = mesh.elm_size;
	elements = mesh.nb_cells;
	nodes = mesh.nb_nodes;
	coordinates = malloc(nodes*spatial_dimension*sizeof(double));
	connectivity =  malloc(elements*nodes_element*sizeof(int));
	element_type = malloc(elements*sizeof(int));
	int i;
	/* filling in nodes cordinates array */
	for (i = 0 ; i<nodes ; i++){
		coordinates[spatial_dimension*i] = mesh.nodes[3 * nodeRank];
		if (spatial_dimension > 1)
			coordinates[spatial_dimension*i+1] = mesh.nodes[3 * nodeRank + 1];
		if (spatial_dimension > 2)
			coordinates[spatial_dimension*i+2] = mesh.nodes[3 * nodeRank + 2];
	}
	/* filling in connectivity array ( in case of triangle element_i --> [node1, node2, node3] ) */
	for (i = 0 ; i<elements*nodes_element ; i++){ 
		connectivity[i] = mesh.cells[i] + 1; // +1 because fortran style
	}
	/* if the element is a triangle fill in element_type array by VTK_TRIANGLE (= 5) */
	if (mesh.elm_size == 3) {
		for (i = 0 ; i<nodes*spatial_dimension ; i++){
			element_type[i] = 5;
		}
	}
}

### dans /home/elafrit/Bureau/src_libmultiscale/simulpack/pmechanics/partitionate.h

ajouter la declaration:
void partition_mesh_hg(int partnbr);

### dans /home/elafrit/Bureau/src_libmultiscale/simulpack/pmechanics/partitionate.c

ajouter la fonction:
void partition_mesh_hg(int partnbr){

  if (adlib_myrank != 0) return;
  fprintf(stderr,"(%d) partitioning mesh in %d pieces\n",adlib_myrank,partnbr);
  fprintf(stderr,"HYPERGRAPH PARTITIONING WITH FIX CELLS\n");
  
  /* generating bridging hypergraph */  
  Hypergraph brghg;
  generateBridgingHypergraph(&mesh, &graph, &atombox, &brghg, 
			     normal_vtx_wght, normal_he_cost, bridging_vtx_wght, bridging_he_cost, coupling_he_cost);
  // printHypergraph(&hg);
  // printHypergraph(&brghg);
  
  int ubfactor_brghg = UBFACTOR;
  int part_brghg[brghg.nvtxs];
  int edgecut_brghg = 0;
  
  partitionHypergraph(&brghg, nparts, ubfactor_brghg, part_brghg, &edgecut_brghg);
  
  /* partitioning with fix cells */
  int ubfactor_fixhg = UBFACTOR;
  int part_fixhg[hg.nvtxs];
  int edgecut_fixhg = 0;
  for(int i = 0 ; i < hg.nvtxs ; i++) part_fixhg[i] = -1; // free cells
  for(int i = 0 ; i < brghg.nvtxs ; i++) part_fixhg[i] = part_brghg[i]; // pre-assigned cells
  
  partitionHypergraphWithFixCells(&hg, nparts, ubfactor_fixhg, part_fixhg, &edgecut_fixhg);
  
  /* add hypergraph partitioning variable */
  double parttmp_fixhg[hg.nvtxs];
  for(int i = 0 ; i < hg.nvtxs ; i++) parttmp_fixhg[i] = part_fixhg[i];
  
  addMeshCellVariable(&mesh, CELL_FIELD_FIXHGPART_ID, 1, parttmp_fixhg); 

  /*use parttab for the following functions*/
  partab = parttmp_fixhg;
  fprintf(stderr,"(%d) partitioning of mesh in %d pieces : done\n",adlib_myrank,adlib_world_size);
}