Make it compile
Implement fred2fcart()
Implement direct_sum()
Add forces