int main(int argc, char* argv[])
{
if (argc != 2) {
std::cout << "usage: " << argv[0] << " configfile" << std::endl;
return 1;
}
#ifdef _OPENMP
std::cout << "# of OpenMP threads = " << omp_get_max_threads() << std::endl;
#endif
Config conf;
std::cout << std::endl << "Read config file: " << argv[1] << std::endl;
conf.load(argv[1]);
conf.print();
std::cout << std::endl << "Polylib setting: " << conf.polylibConf << std::endl;
Polylib* pl = Polylib::get_instance();
if (pl->load(conf.polylibConf)) return 1;
pl->show_group_hierarchy();
RepairPolygonData(pl);
double org[3] = { 0.0, 0.0, 0.0 };
double pitch[3] = { 1.0/conf.ndim[0], 1.0/conf.ndim[1], 1.0/conf.ndim[2] };
GridAccessor* grid = new Node(org, pitch);
size_t nnode[3] = { conf.ndim[0]+1, conf.ndim[1]+1, conf.ndim[2]+1 };
CutPosArray *cutPosArray;
if (conf.cutPosType == "CutPos32") cutPosArray = new CutPos32Array(nnode);
if (conf.cutPosType == "CutPos8") cutPosArray = new CutPos8Array(nnode);
CutBidArray *cutBidArray;
if (conf.cutBidType == "CutBid8") cutBidArray = new CutBid8Array(nnode);
if (conf.cutBidType == "CutBid5") cutBidArray = new CutBid5Array(nnode);
CutNormalArray* cutNormalArray = new CutNormalArray(nnode);
std::cout << std::endl << "CalcCutInfo: " << std::endl;
int ret = CalcCutInfo(conf.ista, conf.nlen,
grid, pl, cutPosArray, cutBidArray, cutNormalArray);
std::cout << "return code = " << ret << std::endl;
if (conf.output != "") {
outputVtk(conf.output, grid, cutPosArray, cutBidArray,
cutNormalArray, conf.reverseNormal);
}
delete grid;
delete cutPosArray;
delete cutBidArray;
delete cutNormalArray;
return 0;
}
int main(int argc, char* argv[])
{
if (argc != 2) {
std::cout << "usage: " << argv[0] << " configfile" << std::endl;
return 1;
}
Config conf;
std::cout << std::endl << "Read config file: " << argv[1] << std::endl;
conf.load(argv[1]);
conf.print();
std::cout << std::endl << "Polylib setting: " << conf.polylibConf << std::endl;
Polylib* pl = Polylib::get_instance();
pl->load(conf.polylibConf);
pl->show_group_hierarchy();
RepairPolygonData(pl);
CutPosOctree *cutPos;
if (conf.cutPosType == "CutPos32") cutPos= new CutPos32Octree(conf.dIndex);
if (conf.cutPosType == "CutPos8") cutPos= new CutPos8Octree(conf.dIndex);
CutBidOctree *cutBid;
if (conf.cutBidType == "CutBid8") cutBid= new CutBid8Octree(conf.dIndex+cutPos->getSizeInFloat());
if (conf.cutBidType == "CutBid5") cutBid= new CutBid5Octree(conf.dIndex+cutPos->getSizeInFloat());
unsigned dLen = conf.dIndex + cutPos->getSizeInFloat() + cutBid->getSizeInFloat();
SklTree* tree = new SklTree;
tree->CreateTree(conf.nRootCell, conf.nLevel, dLen);
float org[3] = { 0.0, 0.0, 0.0 };
float d[3] = {
1.0 / conf.nRootCell[0],
1.0 / conf.nRootCell[1],
1.0 / conf.nRootCell[2],
};
tree->SetPitch(org, d);
std::cout << std::endl;
int ret;
if (conf.leafCellOnly) {
std::cout << "CalcCutInfoOctreeLeafCell: " << std::endl;
ret = CalcCutInfoOctreeLeafCell(tree, pl, cutPos, cutBid);
} else {
std::cout << "CalcCutInfoOctreeAllCell: " << std::endl;
ret = CalcCutInfoOctreeAllCell(tree, pl, cutPos, cutBid);
// ret = CalcCutInfoOctreeAllCell0(tree, pl, cutPos, cutBid);
}
std::cout << "return code = " << ret << std::endl;
if (ret != 0) return ret;
if (conf.output != "") {
if (conf.leafCellOnly) {
outputVtkLeafCell(conf.output, tree, cutPos, cutBid);
} else {
outputVtkAllCell(conf.output, tree, conf.nLevel, cutPos, cutBid);
}
}
delete cutPos;
delete cutBid;
delete tree;
return 0;
}