@@ -502,7 +502,9 @@ def bfgs_sample_dense(
502
502
503
503
logdet = 2.0 * pt .sum (pt .log (pt .abs (pt .diagonal (Lchol , axis1 = - 2 , axis2 = - 1 ))), axis = - 1 )
504
504
505
- mu = x - pt .batched_dot (H_inv , g )
505
+ with _warnings .catch_warnings ():
506
+ _warnings .simplefilter ("ignore" , category = FutureWarning )
507
+ mu = x - pt .batched_dot (H_inv , g )
506
508
507
509
phi = pt .matrix_transpose (
508
510
# (L, N, 1)
@@ -572,14 +574,16 @@ def bfgs_sample_sparse(
572
574
logdet += pt .sum (pt .log (alpha ), axis = - 1 )
573
575
574
576
# NOTE: changed the sign from "x + " to "x -" of the expression to match Stan which differs from Zhang et al., (2022). same for dense version.
575
- mu = x - (
576
- # (L, N), (L, N) -> (L, N)
577
- pt .batched_dot (alpha_diag , g )
578
- # beta @ gamma @ beta.T
579
- # (L, N, 2J), (L, 2J, 2J), (L, 2J, N) -> (L, N, N)
580
- # (L, N, N), (L, N) -> (L, N)
581
- + pt .batched_dot ((beta @ gamma @ pt .matrix_transpose (beta )), g )
582
- )
577
+ with _warnings .catch_warnings ():
578
+ _warnings .simplefilter ("ignore" , category = FutureWarning )
579
+ mu = x - (
580
+ # (L, N), (L, N) -> (L, N)
581
+ pt .batched_dot (alpha_diag , g )
582
+ # beta @ gamma @ beta.T
583
+ # (L, N, 2J), (L, 2J, 2J), (L, 2J, N) -> (L, N, N)
584
+ # (L, N, N), (L, N) -> (L, N)
585
+ + pt .batched_dot ((beta @ gamma @ pt .matrix_transpose (beta )), g )
586
+ )
583
587
584
588
phi = pt .matrix_transpose (
585
589
# (L, N, 1)
0 commit comments