Replace gaussj implementation