P_link1 and P_link2 do not exist on the same rigid body so this is incorrect. It should be:

```python
P_link2.v2pt_theory(P1, I, B) 
```