File tree Expand file tree Collapse file tree 2 files changed +0
-37
lines changed Expand file tree Collapse file tree 2 files changed +0
-37
lines changed Original file line number Diff line number Diff line change @@ -35,7 +35,6 @@ namespace deform {
3535 void estimateRotations ();
3636 void estimatePositions ();
3737 void prepareLinearSystem ();
38- void computeB ();
3938
4039
4140
Original file line number Diff line number Diff line change @@ -289,40 +289,4 @@ namespace deform {
289289 }
290290 }
291291 }
292-
293-
294- void AsRigidAsPossibleDeformation::computeB ()
295- {
296- /*
297- Mesh &m = *_data->mesh;
298-
299- _data->b.setZero();
300-
301- for (auto viter = m.vertices_begin(); viter != m.vertices_end(); ++viter) {
302- int i = viter->idx();
303-
304- if (m.property(_data->isConstrained, *viter)) {
305- // If this vertex is constraint we enforce lhs = rhs
306- _data->b.row(i) = m.property(_data->constraintLocations, *viter);
307- }
308- else {
309- // else w/2 * (Ri+Rj) * (pi - pj)
310- Eigen::Vector3d sum = Eigen::Vector3d::Zero();
311-
312- for (auto vh = m.voh_begin(*viter); vh != m.voh_end(*viter); ++vh) {
313-
314- Mesh::VertexHandle vj = m.to_vertex_handle(*vh);
315- int j = vj.idx();
316-
317- Eigen::Matrix3d r = (_data->rotations[i] + _data->rotations[j]);
318- Eigen::Vector3d v = _data->points.col(i) - _data->points.col(j);
319- Eigen::Vector3d t = r * v;
320- const double w = m.property(_data->cotanWeights, m.edge_handle(*vh));
321- sum += w * t * 0.5;
322- }
323- _data->b.row(i) = sum.transpose();
324- }
325- }
326- */
327- }
328292}
You can’t perform that action at this time.
0 commit comments