-
-
Save erikzenker/c4dc42c8d5a8c1cd3e5a to your computer and use it in GitHub Desktop.
#include <cstddef> /* NULL */ | |
#include <metis.h> | |
#include <iostream> | |
// Install metis from: | |
// http://glaros.dtc.umn.edu/gkhome/fetch/sw/metis/metis-5.1.0.tar.gz | |
// Build with | |
// g++ metis.cc -lmetis | |
int main(){ | |
idx_t nVertices = 6; | |
idx_t nEdges = 7; | |
idx_t nWeights = 1; | |
idx_t nParts = 2; | |
idx_t objval; | |
idx_t part[nVertices]; | |
// Indexes of starting points in adjacent array | |
idx_t xadj[nVertices+1] = {0,2,5,7,9,12,14}; | |
// Adjacent vertices in consecutive index order | |
idx_t adjncy[2 * nEdges] = {1,3,0,4,2,1,5,0,4,3,1,5,4,2}; | |
// Weights of vertices | |
// if all weights are equal then can be set to NULL | |
idx_t vwgt[nVertices * nWeights]; | |
// int ret = METIS_PartGraphRecursive(&nVertices,& nWeights, xadj, adjncy, | |
// NULL, NULL, NULL, &nParts, NULL, | |
// NULL, NULL, &objval, part); | |
int ret = METIS_PartGraphKway(&nVertices,& nWeights, xadj, adjncy, | |
NULL, NULL, NULL, &nParts, NULL, | |
NULL, NULL, &objval, part); | |
std::cout << ret << std::endl; | |
for(unsigned part_i = 0; part_i < nVertices; part_i++){ | |
std::cout << part_i << " " << part[part_i] << std::endl; | |
} | |
return 0; | |
} |
Works fine for me...
GCC-8
Linux
Why I also have this problem?How do you deal with it?
The result is
1
0 1
1 4294967296
2 4294967296
3 0
4 0
5 0
Thanks very well!
@Yingshiyu These high numbers in the partitioning array look like uninitialized memory. Best practice is to initialize memory before using it. Some platforms and compilers do it for you, this could explain the different result of you guys. I replaced the C arrays with proper initialized c++ vectors. Can you try again?
#include <cstddef> /* NULL */
#include <metis.h>
#include <iostream>
#include <vector>
// Install metis from:
// http://glaros.dtc.umn.edu/gkhome/fetch/sw/metis/metis-5.1.0.tar.gz
// Build with
// g++ metis.cc -lmetis
int main(){
idx_t nVertices = 6;
idx_t nEdges = 7;
idx_t nWeights = 1;
idx_t nParts = 2;
idx_t objval;
std::vector<idx_t> part(nVertices, 0);
// Indexes of starting points in adjacent array
std::vector<idx_t> xadj = {0,2,5,7,9,12,14};
// Adjacent vertices in consecutive index order
std::vector<idx_t> adjncy = {1,3,0,4,2,1,5,0,4,3,1,5,4,2};
// Weights of vertices
// if all weights are equal then can be set to NULL
std::vector<idx_t> vwgt(nVertices * nWeights, 0);
int ret = METIS_PartGraphKway(&nVertices,& nWeights, xadj.data(), adjncy.data(),
NULL, NULL, NULL, &nParts, NULL,
NULL, NULL, &objval, part.data());
std::cout << ret << std::endl;
for(unsigned part_i = 0; part_i < part.size(); part_i++){
std::cout << part_i << " " << part[part_i] << std::endl;
}
return 0;
}
@erikzenker
Thank you very much. But there are also have some problems when I use your code.
error: in C++98 'xadj' must be initialized by constructor, not by'{}'
error: could not convert{0,2,5,7,9,12,14} from to 'std::vector'
And I think this problem is related to the vector initialization(vector can not be initialized by {}). Therefore, I have modified it like:
idx_t a[nVertices+1]={0,2,5,7,9,12,14};
std::vector<idx_t> xadj(a,a+nVertices+1);
But the result is not change
1
0 1
1 4294967296
2 4294967296
3 0
4 0
5 0
Besides I cannot change nParts beyond 2.
My task is to use this API in my program, therefore, I want to try this API by a little example.
I really hope to get your reply again.
Thanks very much!
Which Compiler and oparating System are you using? Initialization with {} comes with c++11. I do not suggest to use c++98.
@erikzenker
OK, I have done it with -std=c+=11!
thanks very well!
@erikzenker
Sorry to disturb you again.
Your example can run on my computer and the answer is right. I want to ask are you suffer this problem?
If I am not declare vector like vector = {0,2,5,7,9,12,14}, because I don't know the size of this map or this map is too large.
Therefore I use std::vector<idx_t> xadj
xadj.push_back(i) to get this vector.
Finally I output the value of vector is right. But there are some problems
munmap_chunk():invalid pointer:0x00000000013fdc70
or
free(): invalid size:0x0000000000lacccc0
Whether should I change the source code? I think it perhaps about 'new' and 'delete'?
I want to include this library/code in a makefile project, so how can I add the tag "-lmetis" at the compilation? I will be building the project using makefile.
Running METIS is really difficult as I am a newbie here and C language is not my programming language. And even though I am using the METIS manual, I haven't been able to run the program. I want to use the road network dataset to partition it. And maybe use the output file as the input in the Apache spark to achieve distributed processing. Anyone, could you please share the tutorial or something helpful so that I can follow the pattern or understand how to use it? PLEASE! HELP!
Works fine for me...
GCC-8
Linux