@@ -418,8 +418,8 @@ def solve_semi_dual_entropic(a, b, M, reg, method, numItermax=10000, lr=None,
418418 return None
419419
420420 opt_alpha = c_transform_entropic (b , M , reg , opt_beta )
421- pi = (np .exp ((opt_alpha [:, None ] + opt_beta [None , :] - M [:, :]) / reg )
422- * a [:, None ] * b [None , :])
421+ pi = (np .exp ((opt_alpha [:, None ] + opt_beta [None , :] - M [:, :]) / reg ) *
422+ a [:, None ] * b [None , :])
423423
424424 if log :
425425 log = {}
@@ -520,8 +520,8 @@ def batch_grad_dual(a, b, M, reg, alpha, beta, batch_size, batch_alpha,
520520 arXiv preprint arxiv:1711.02283.
521521 '''
522522
523- G = - (np .exp ((alpha [batch_alpha , None ] + beta [None , batch_beta ]
524- - M [batch_alpha , :][:, batch_beta ]) / reg ) *
523+ G = - (np .exp ((alpha [batch_alpha , None ] + beta [None , batch_beta ] -
524+ M [batch_alpha , :][:, batch_beta ]) / reg ) *
525525 a [batch_alpha , None ] * b [None , batch_beta ])
526526 grad_beta = np .zeros (np .shape (M )[1 ])
527527 grad_alpha = np .zeros (np .shape (M )[0 ])
@@ -702,8 +702,8 @@ def solve_dual_entropic(a, b, M, reg, batch_size, numItermax=10000, lr=1,
702702
703703 opt_alpha , opt_beta = sgd_entropic_regularization (a , b , M , reg , batch_size ,
704704 numItermax , lr )
705- pi = (np .exp ((opt_alpha [:, None ] + opt_beta [None , :] - M [:, :]) / reg )
706- * a [:, None ] * b [None , :])
705+ pi = (np .exp ((opt_alpha [:, None ] + opt_beta [None , :] - M [:, :]) / reg ) *
706+ a [:, None ] * b [None , :])
707707 if log :
708708 log = {}
709709 log ['alpha' ] = opt_alpha
0 commit comments