manolas
/
vcglib
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Added Reserve in sparse matrix init for speeding up poisson parametrization

This commit is contained in:
Paolo Cignoni 2017-09-18 11:26:40 +02:00
parent 288d86b160
commit d15e5fbb91
1 changed files with 8 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

17
vcg/complex/algorithms/parametrization/poisson_solver.h
View File

@ -92,15 +92,13 @@ class PoissonSolver
void FindFarthestVert(VertexType* &v0, VertexType* &v1) void FindFarthestVert(VertexType* &v0, VertexType* &v1)
{ {
v0=NULL; v0=v1=NULL;
v1=NULL; const int bestAxis = mesh.bbox.MaxDim();
int bestAxis = mesh.bbox.MaxDim(); for(VertexType &vv : mesh.vert) {
for(VertexType &vv : mesh.vert)
{
if(vv.P()[bestAxis] <= mesh.bbox.min[bestAxis]) v0 = &vv; if(vv.P()[bestAxis] <= mesh.bbox.min[bestAxis]) v0 = &vv;
if(vv.P()[bestAxis] >= mesh.bbox.max[bestAxis]) v1 = &vv; if(vv.P()[bestAxis] >= mesh.bbox.max[bestAxis]) v1 = &vv;
} }
assert(v0 && v1); assert( (v0!=v1) && v0 && v1);
} }
///set the value of b of the system Ax=b ///set the value of b of the system Ax=b
@ -398,6 +396,7 @@ class PoissonSolver
{ {
//--- Allocates the data for Ax=b //--- Allocates the data for Ax=b
A=Eigen::SparseMatrix<double>(total_size, total_size); // A A=Eigen::SparseMatrix<double>(total_size, total_size); // A
A.reserve(Eigen::VectorXi::Constant(total_size,32)); // This prealloaction trick greatly speed up the acc
b = Eigen::VectorXd::Zero(total_size); // x and b b = Eigen::VectorXd::Zero(total_size); // x and b
} }
@ -673,12 +672,12 @@ public:
PoissonSolver(MeshType &_mesh):mesh(_mesh) PoissonSolver(MeshType &_mesh):mesh(_mesh)
{ {
assert(mesh.vert.size()>3); assert(mesh.vert.size()>=3);
assert(mesh.face.size()>1); assert(mesh.face.size()>=1);
} }
}; // end class }; // end class
} //End Namespace Tri } // End Namespace tri
} // End Namespace vcg } // End Namespace vcg
#endif #endif