Implement a formula interface

[evangorstein]

• Refactor hdmm() so it does argument checks and construction for both input of matrices and input of formula
• Move actual fitting to _hdmm() workhorse

[evangorstein]

@ajinkya-k