code: def loadDolfin(filename, c="gold", alpha=0.5, wire=None, bc=None):
      
      if not os.path.exists(filename):
          colors.printc("~noentry Error in loadDolfin: Cannot find", filename, c=1)
          return None
      import xml.etree.ElementTree as et
  
      if filename.endswith(".gz"):
          import gzip
  
          inF = gzip.open(filename, "rb")
          outF = open("/tmp/filename.xml", "wb")
          outF.write(inF.read())
          outF.close()
          inF.close()
          tree = et.parse("/tmp/filename.xml")
      else:
          tree = et.parse(filename)
      coords, connectivity = [], []
      for mesh in tree.getroot():
          for elem in mesh:
              for e in elem.findall("vertex"):
                  x = float(e.get("x"))
                  y = float(e.get("y"))
                  ez = e.get("z")
                  if ez is None:
                      coords.append([x, y])
                  else:
                      z = float(ez)
                      coords.append([x, y, z])
  
              tets = elem.findall("tetrahedron")
              if not len(tets):
                  tris = elem.findall("triangle")
                  for e in tris:
                      v0 = int(e.get("v0"))
                      v1 = int(e.get("v1"))
                      v2 = int(e.get("v2"))
                      connectivity.append([v0, v1, v2])
              else:
                  for e in tets:
                      v0 = int(e.get("v0"))
                      v1 = int(e.get("v1"))
                      v2 = int(e.get("v2"))
                      v3 = int(e.get("v3"))
                      connectivity.append([v0, v1, v2, v3])
  
      poly = buildPolyData(coords, connectivity)
      return Actor(poly, c, alpha, True, bc)
docs: Reads a `Fenics/Dolfin` file format. Return an ``Actor(vtkActor)`` object.