@@ -461,15 +461,15 @@ def fit(self, X: np.ndarray, Z: np.ndarray, d=None, g=None, dK: bool = True, aut
461461 # raise ValueError("derivative info not in GPsep; use newGPsep with dK=True")
462462
463463 # New: mleGPsep_optimize starts here:
464- print ( f"Starting MLE with d= { self . d } , g= { self . g } " )
464+
465465 # generate starting point p
466466 p = np .concatenate ([self .d , [self .g ]])
467- print (f"Starting point: { p } " )
468467 bounds = [(tmin [i ], tmax [i ]) for i in range (len (p ))]
469- print (f"bounds: { bounds } " )
470- print (f"p: { p } " )
471- # print(f"self.X: {self.X}")
472- # print(f"self.Z: {self.Z}")
468+ if self .verbosity > 0 :
469+ print (f"Starting MLE with d={ self .d } , g={ self .g } " )
470+ print (f"Starting point: { p } " )
471+ print (f"bounds: { bounds } " )
472+ print (f"p: { p } " )
473473 X = copy .deepcopy (self .X )
474474 Z = copy .deepcopy (self .Z )
475475
@@ -480,14 +480,17 @@ def gradient(par):
480480 return gradnlsep (par , X , Z , self .gradnlsep_method )
481481
482482 result = run_minimize_with_restarts (objective = objective , gradient = gradient , x0 = p , bounds = bounds , n_restarts_optimizer = self .n_restarts_optimizer , maxit = self .maxit , verb = self .verbosity )
483- print (f"result: { result } " )
484483
485484 d = result .x [:- 1 ]
486485 g = result .x [- 1 ]
487- print (f"Optimized d: { d } , g: { g } " )
486+
487+
488488 # set new parameters and build
489489 self .set_new_params (d , g )
490- print (f"Updated d: { self .d } , g: { self .g } " )
490+ if self .verbosity > 0 :
491+ print (f"result: { result } " )
492+ print (f"Optimized d: { d } , g: { g } " )
493+ print (f"Updated d: { self .d } , g: { self .g } " )
491494 self .build ()
492495 new_theta = np .concatenate ((self .get_d (), [self .get_g ()]))
493496 if np .sqrt (np .mean ((result .x - new_theta ) ** 2 )) > np .sqrt (np .finfo (float ).eps ):
@@ -531,7 +534,6 @@ def calc_ZtKiZ(self) -> None:
531534 KiZ = np .dot (self .Ki , Z )
532535 phi = np .dot (Z .T , KiZ )
533536 self .phi = phi [0 , 0 ]
534- print (f"phi: { self .phi } " )
535537 self .KiZ = KiZ
536538
537539 def build (self ) -> None :
0 commit comments