This call obviously needs the chaco package to run in
parallel since otherwise I will get a warning to
configure with --download-chaco.
So I ran it with chaco installed. The code works but
the domain decomposition is valid but far from good.
In the attached image, I am showing the domain
decomposition for a 2-dimensional 1x1 box with two
processes, using the following statements:
MeshCreatePCICE(comm, 3,
"bratu_2d.nodes","bratu_2d.lcon",PETSC_FALSE,PETSC_FALSE,
&serialMesh);
MeshDistribute(serialMesh, PETSC_NULL, ¶llelMesh);
I am wondering if it is possible to improve the domain
decomposition by passing some command arguments. I
realized that there are Chaco options such as
-mat_partitioning_chaco_global (found at
http://www-unix.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-current/docs/manualpages/MatOrderings/MAT_PARTITIONING_Chaco.html)
. However, when I pass these options to my code, it is
not recognized. I wonder how chaco is used in Petsc
and how I can change its behavior.
In addition, is it possible to use other graph
decomposition packages such as ParMetis to implement
MeshDistribute()? I tried to enable ParMetis without
chaco but the code didn't run and warned me to install
chaco.
Any advice is valuable. Thank you very much.
Shi
____________________________________________________________________________________
Finding fabulous fares is fun.
Let Yahoo! FareChase search your favorite travel sites to find flight and hotel bargains.
http://farechase.yahoo.com/promo-generic-14795097