You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
 
 
 
 
 
 

26 lines
910 B

#include <test_common.h>
#include <igl/boundary_conditions.h>
#include <igl/readMESH.h>
#include <igl/writeDMAT.h>
#include <igl/readTGF.h>
#include <igl/bbw.h>
TEST_CASE("bbw: decimated_knight", "[igl]" "[slow]")
{
Eigen::MatrixXd V,C;
Eigen::MatrixXi T,F,E;
igl::readMESH(test_common::data_path("decimated-knight.mesh"),V,T,F);
igl::readTGF(test_common::data_path("decimated-knight.tgf"),C,E);
Eigen::MatrixXd W_groundtruth, Was, Wmo;
igl::readDMAT(
test_common::data_path("decimated-knight-matlab-active-set.dmat"),W_groundtruth);
Eigen::VectorXi b;
Eigen::MatrixXd bc;
igl::boundary_conditions(V,T,C,Eigen::VectorXi(),E,Eigen::MatrixXi(),Eigen::MatrixXi(),b,bc);
igl::BBWData params;
params.active_set_params.max_iter = 100;
igl::bbw(V,T,b,bc,params,Was);
// igl::writeDMAT("decimated-knight-as.dmat",Was);
REQUIRE (1e-4 > (Was-W_groundtruth).array().abs().maxCoeff());
}