add the dual KL method

just out of curiosity. are there any update on the linalg stuff for solving linear systems? @lambday