extern shapeop_err shapeop_editConstraint(ShapeOpSolver *op, const char *constraintType, int constraint_id, const ShapeOpScalar *scalars, int nb_scl) { if (strcmp(constraintType, "EdgeStrain") == 0) { auto c = std::dynamic_pointer_cast<ShapeOp::EdgeStrainConstraint>(op->s->getConstraint(constraint_id)); if (!c) { return SO_UNMATCHING_CONSTRAINT_ID; } if (nb_scl != 3) { return SO_INVALID_ARGUMENT_LENGTH; } c->setEdgeLength(scalars[0]); c->setRangeMin(scalars[1]); c->setRangeMax(scalars[2]); return SO_SUCCESS; } if (strcmp(constraintType, "TriangleStrain") == 0) { auto c = std::dynamic_pointer_cast<ShapeOp::TriangleStrainConstraint>(op->s->getConstraint(constraint_id)); if (!c) { return SO_UNMATCHING_CONSTRAINT_ID; } if (nb_scl != 2) { return SO_INVALID_ARGUMENT_LENGTH; } c->setRangeMin(scalars[0]); c->setRangeMax(scalars[1]); return SO_SUCCESS; } if (strcmp(constraintType, "TetrahedronStrain") == 0) { auto c = std::dynamic_pointer_cast<ShapeOp::TetrahedronStrainConstraint>(op->s->getConstraint(constraint_id)); if (!c) { return SO_UNMATCHING_CONSTRAINT_ID; } if (nb_scl != 2) { return SO_INVALID_ARGUMENT_LENGTH; } c->setRangeMin(scalars[0]); c->setRangeMax(scalars[1]); return SO_SUCCESS; } if (strcmp(constraintType, "Area") == 0) { auto c = std::dynamic_pointer_cast<ShapeOp::AreaConstraint>(op->s->getConstraint(constraint_id)); if (!c) { return SO_UNMATCHING_CONSTRAINT_ID; } if (nb_scl != 2) { return SO_INVALID_ARGUMENT_LENGTH; } c->setRangeMin(scalars[0]); c->setRangeMax(scalars[1]); return SO_SUCCESS; } if (strcmp(constraintType, "Volume") == 0) { auto c = std::dynamic_pointer_cast<ShapeOp::VolumeConstraint>(op->s->getConstraint(constraint_id)); if (!c) { return SO_UNMATCHING_CONSTRAINT_ID; } if (nb_scl != 2) { return SO_INVALID_ARGUMENT_LENGTH; } c->setRangeMin(scalars[0]); c->setRangeMax(scalars[1]); return SO_SUCCESS; } if (strcmp(constraintType, "Bending") == 0) { auto c = std::dynamic_pointer_cast<ShapeOp::BendingConstraint>(op->s->getConstraint(constraint_id)); if (!c) { return SO_UNMATCHING_CONSTRAINT_ID; } if (nb_scl != 2) { return SO_INVALID_ARGUMENT_LENGTH; } c->setRangeMin(scalars[0]); c->setRangeMax(scalars[1]); return SO_SUCCESS; } if (strcmp(constraintType, "Closeness") == 0) { auto c = std::dynamic_pointer_cast<ShapeOp::ClosenessConstraint>(op->s->getConstraint(constraint_id)); if (!c) { return SO_UNMATCHING_CONSTRAINT_ID; } if (nb_scl != 3) { return SO_INVALID_ARGUMENT_LENGTH; } Eigen::Map<const ShapeOp::Vector3> p(scalars, 3, 1); c->setPosition(p); return SO_SUCCESS; } if (strcmp(constraintType, "Similarity") == 0 || strcmp(constraintType, "Rigid") == 0) { auto c = std::dynamic_pointer_cast<ShapeOp::SimilarityConstraint>(op->s->getConstraint(constraint_id)); if (!c) { return SO_UNMATCHING_CONSTRAINT_ID; } int nI = static_cast<int>(c->nIndices()); if ((nb_scl % (nI * 3)) != 0) { return SO_INVALID_ARGUMENT_LENGTH; } std::vector<ShapeOp::Matrix3X> shapes; int nShapes = nb_scl / (nI * 3); for (int i = 0; i < nShapes; ++i) { Eigen::Map<const ShapeOp::Matrix3X> s(scalars + i * nI * 3, 3, nI); shapes.push_back(s); } c->setShapes(shapes); return SO_SUCCESS; } if (strcmp(constraintType, "Angle") == 0) { auto c = std::dynamic_pointer_cast<ShapeOp::AngleConstraint>(op->s->getConstraint(constraint_id)); if (!c) { return SO_UNMATCHING_CONSTRAINT_ID; } if (nb_scl != 2) { return SO_INVALID_ARGUMENT_LENGTH; } c->setMinAngle(scalars[0]); c->setMaxAngle(scalars[1]); return SO_SUCCESS; } return SO_INVALID_CONSTRAINT_TYPE; }
extern void shapeop_editEdgeStrainConstraint(ShapeOpSolver *op, int constraint_id, ShapeOpScalar length) { auto c = std::dynamic_pointer_cast<ShapeOp::EdgeStrainConstraint>(op->s->getConstraint(constraint_id)); //TODO: this will need to be robustify c->setEdgeLength(length); }