Initial implementation of Biloys solver