public TriMesh Recon(TriMesh mesh, int reconNum)
{
Eigen eigen = EigenManager.Instance.ComputeEigen(mesh,
EnumLaplaceMatrix.LapalceCot, reconNum + 10);
List <double> X = TriMeshUtil.GetX(mesh);
List <double> Y = TriMeshUtil.GetY(mesh);
List <double> Z = TriMeshUtil.GetZ(mesh);
List <double> factorX = new List <double>(reconNum);
List <double> factorY = new List <double>(reconNum);
List <double> factorZ = new List <double>(reconNum);
for (int i = 0; i < reconNum; i++)
{
List <double> eigenVector = eigen.GetEigenVector(i);
double valueX = TriMeshUtil.Multiply(eigenVector, X);
double valueY = TriMeshUtil.Multiply(eigenVector, Y);
double valueZ = TriMeshUtil.Multiply(eigenVector, Z);
factorX.Add(valueX);
factorY.Add(valueY);
factorZ.Add(valueZ);
}
double[] reconX = new double[mesh.Vertices.Count];
double[] reconY = new double[mesh.Vertices.Count];
double[] reconZ = new double[mesh.Vertices.Count];
for (int i = 0; i < reconNum; i++)
{
List <double> eigenVector = eigen.GetEigenVector(i);
List <double> tempX = TriMeshUtil.Multiply(eigenVector, factorX[i]);
reconX = TriMeshUtil.Add(reconX, tempX);
List <double> tempY = TriMeshUtil.Multiply(eigenVector, factorY[i]);
reconY = TriMeshUtil.Add(reconY, tempY);
List <double> tempZ = TriMeshUtil.Multiply(eigenVector, factorZ[i]);
reconZ = TriMeshUtil.Add(reconZ, tempZ);
}
for (int i = 0; i < mesh.Vertices.Count; i++)
{
mesh.Vertices[i].Traits.Position.x = reconX[i];
mesh.Vertices[i].Traits.Position.y = reconY[i];
mesh.Vertices[i].Traits.Position.z = reconZ[i];
}
return(mesh);
}
public double[][] ComputeLaplacianBasic(SparseMatrix L, TriMesh mesh)
{
if (L == null)
{
throw new Exception("Laplacian matrix is null");
}
int n = mesh.Vertices.Count;
double[] coordinate = new double[n];
double[][] lap = new double[3][];
coordinate = TriMeshUtil.GetX(mesh).ToArray();
lap[0] = L.Multiply(coordinate);
coordinate = TriMeshUtil.GetY(mesh).ToArray();
lap[1] = L.Multiply(coordinate);
coordinate = TriMeshUtil.GetZ(mesh).ToArray();
lap[2] = L.Multiply(coordinate);
return(lap);
}
public double[] ComputeY(TriMesh mesh)
{
return(TriMeshUtil.GetY(mesh).ToArray());
}