let store_matrix f dm =
    let oc = open_out_bin f in
    output_value oc dm;
    close_out oc