Function generateSDF

Tools/VolumeSampling/main.cpp:406–441

}

void generateSDF(SPH::TriangleMesh &mesh)
{
	//////////////////////////////////////////////////////////////////////////
	// Generate distance field of object using Discregrid
	//////////////////////////////////////////////////////////////////////////
#ifdef USE_DOUBLE
	Discregrid::TriangleMesh sdfMesh(&mesh.getVertices()[0][0], mesh.getFaces().data(), mesh.numVertices(), mesh.numFaces());
#else
	// if type is float, copy vector to double vector
	std::vector<double> doubleVec;
	doubleVec.resize(3 * mesh.numVertices());
	for (unsigned int i = 0; i < mesh.numVertices(); i++)
		for (unsigned int j = 0; j < 3; j++)
			doubleVec[3 * i + j] = mesh.getVertices()[i][j];
	Discregrid::TriangleMesh sdfMesh(&doubleVec[0], mesh.getFaces().data(), mesh.numVertices(), mesh.numFaces());
#endif

	Discregrid::TriangleMeshDistance md(sdfMesh);
	Eigen::AlignedBox3d domain;
	domain.extend(bbmin.cast<double>());
	domain.extend(bbmax.cast<double>());
	//domain.max() += 1.0e-3 * domain.diagonal().norm() * Eigen::Vector3d::Ones();
	//domain.min() -= 1.0e-3 * domain.diagonal().norm() * Eigen::Vector3d::Ones();
	domain.max() += Eigen::Vector3d::Ones() * radius * 4.0;
	domain.min() -= Eigen::Vector3d::Ones() * radius * 4.0;

	LOG_INFO << "Set SDF resolution: " << resolutionSDF[0] << ", " << resolutionSDF[1] << ", " << resolutionSDF[2];
	distanceField = std::make_shared<Discregrid::CubicLagrangeDiscreteGrid>(domain, std::array<unsigned int, 3>({ resolutionSDF[0], resolutionSDF[1], resolutionSDF[2] }));
	auto func = Discregrid::DiscreteGrid::ContinuousFunction{};

	// invert the distance field since the particles should stay inside
	func = [&md](Eigen::Vector3d const& xi) {return -md.signed_distance(xi).distance; };

	LOG_INFO << "Generate SDF";
	distanceField->addFunction(func, false);
}