本文整理汇总了C++中TPZCompMesh::ComputeFillIn方法的典型用法代码示例。如果您正苦于以下问题:C++ TPZCompMesh::ComputeFillIn方法的具体用法?C++ TPZCompMesh::ComputeFillIn怎么用?C++ TPZCompMesh::ComputeFillIn使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类TPZCompMesh
的用法示例。
在下文中一共展示了TPZCompMesh::ComputeFillIn方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: main
int main(int argc, char *argv[])
{
std::string dirname = PZSOURCEDIR;
#ifdef LOG4CXX
InitializePZLOG();
#endif
gRefDBase.InitializeRefPatterns();
TPZReadGIDGrid readGid;
TPZGeoMesh *gmesh = readGid.GeometricGIDMesh("../rail.dump");
{
std::ofstream out("../GmeshOrig.vtk");
TPZVTKGeoMesh::PrintGMeshVTK(gmesh, out);
}
InsertBoundaryElements(gmesh);
// CASO NAO TEM TRILHO 2
SwitchBCTrilho2(gmesh);
AddZZeroFaces(gmesh);
RefineGMesh(gmesh);
{
std::ofstream out("../GmeshRef.vtk");
TPZVTKGeoMesh::PrintGMeshVTK(gmesh, out);
}
const int porder = 2;
TPZCompMesh * cmesh = ComputationalElasticityMesh3D(gmesh,porder);
AdaptPOrders(cmesh, -400., 2);
// toto
// {
// std::set<int> mats;
// mats.insert(bcbottom);
// TPZManVector<STATE> force(1,0.);
// force = cmesh->Integrate("StressZ", mats);
// std::cout << "force = " << force << std::endl;
// }
{
std::ofstream out("../CmeshRef.vtk");
TPZVTKGeoMesh::PrintCMeshVTK(cmesh, out);
}
{
TPZFMatrix<STATE> visualf;
cmesh->ComputeFillIn(150, visualf);
VisualMatrix(visualf,"../VisualMatrixBefore.vtk");
}
TPZAnalysis an;
TPZAutoPointer<TPZRenumbering> renumber;
// renumber = new TPZSloan;
// renumber = new TPZCutHillMcKee;
renumber = new TPZMetis;
an.SetRenumber(renumber);
an.SetCompMesh(cmesh, true);
#ifdef PZDEBUG
{
std::ofstream out("../gmesh.txt");
gmesh->Print(out);
}
{
std::ofstream out("../cmesh.txt");
cmesh->Print(out);
}
#endif
{
TPZFMatrix<STATE> visualf;
cmesh->ComputeFillIn(150, visualf);
VisualMatrix(visualf,"../MatrixMetis.vtk");
}
std::cout << "NEquations " << an.Solution().Rows() << std::endl;
SolveSist(&an,cmesh);
std::set<int> mats;
mats.insert(bcbottom);
TPZManVector<STATE> Force(1,0.);
Force = cmesh->Integrate("StressZ", mats);
std::cout << "Integrated sigma_Z " << Force << std::endl;
std::cout << "Post processing" << std::endl;
PostProcessElasticity(an, "../postProc.vtk");
std::cout << "Finished\n";
// delete cmesh;
// delete gmesh;
return EXIT_SUCCESS;
//.........这里部分代码省略.........