Computing quasiconformal maps using an auxiliary metric and discrete curvature flow