Romain Prévost has kindly posted his reimplementation of our paper Bounded Biharmonic Weights for Real-time Deformation. He's calling it VoxelBBW. The main advantage/disadvantage of his implementation is that it uses a voxel grid fit to the input shape, rather than a tetrahedral mesh. This is an advantage because tetrahedral meshing of an input triangle mesh is elusive in practice, where as voxel grids always give you something. It's a disadvantage because voxel grids are dense, leading to slower computation. It all depends on your input, your application. This code might be perfect for you!