sig
  type t = private { data : Lacaml.D.mat array; n : int; }
  val create : Lacaml.D.mat array -> Block_diag.t
  val copy : Block_diag.t -> Block_diag.t
  val potrf : ?jitter:float -> Block_diag.t -> unit
  val potri : ?jitter:float -> ?factorize:bool -> Block_diag.t -> unit
end