add residuals

Author: tomvanmele

plugin/FF_solve_fd.py

@@ -1,6 +1,7 @@
 #! python3
 import rhinoscriptsyntax as rs  # type: ignore  # noqa: F401
 
+from compas.geometry import Vector
 from compas.scene import Scene
 from compas_fd.solvers import fd_constrained_numpy
 from compas_fofin.datastructures import CableMesh
@@ -51,7 +52,7 @@ def RunCommand(is_interactive):
         attr["x"] = result.vertices[index, 0]
         attr["y"] = result.vertices[index, 1]
         attr["z"] = result.vertices[index, 2]
-        attr["residual"] = result.residuals[index]
+        attr["residual"] = Vector(*result.residuals[index])
 
     for index, (edge, attr) in enumerate(mesh.edges(data=True)):
         attr["_f"] = result.forces[index]

plugin/FoFin.rhproj

@@ -2,7 +2,7 @@
   "id": "5df3b821-822e-49e0-b24e-aebbe671c3d1",
   "identity": {
     "name": "COMPAS-FormFinder",
-    "version": "0.1.34-beta",
+    "version": "0.1.36-beta",
     "publisher": {
       "email": "tom.v.mele@gmail.com",
       "name": "Tom Van Mele",

src/compas_fofin/datastructures/cablemesh.py

@@ -43,13 +43,13 @@ def __init__(self, constraints=None, **kwargs):
 
     def vertex_residual(self, vertex):
         residual = self.vertex_attribute(vertex, "residual")
-        if not residual:
+        if residual is None:
             residual = Vector(0, 0, 0)
         return residual
 
     def vertex_load(self, vertex):
         load = self.vertex_attribute(vertex, "load")
-        if not load:
+        if load is None:
             load = Vector(0, 0, 0)
         return load