physbo 2.1.0__tar.gz → 2.2.0__tar.gz

{physbo-2.1.0 → physbo-2.2.0}/PKG-INFO

@@ -1,6 +1,6 @@
-Metadata-Version: 2.2
+Metadata-Version: 2.4
 Name: physbo
-Version: 2.1.0
+Version: 2.2.0
 Summary: optimization tool for PHYSics based on Bayesian Optimization
 Home-page: https://github.com/issp-center-dev/PHYSBO
 Author: PHYSBO developers

{physbo-2.1.0 → physbo-2.2.0}/physbo/__init__.py

@@ -14,4 +14,4 @@ from . import predictor
 from .predictor import base_predictor
 from .variable import variable
 
-__version__ = "2.1.0"
+__version__ = "2.2.0"

{physbo-2.1.0 → physbo-2.2.0}/physbo/gp/__init__.py

@@ -11,5 +11,6 @@ from . import lik
 
 from .core import prior
 from .core import model
+from .core import sfs
 from .core import learning
 from .predictor import predictor

{physbo-2.1.0 → physbo-2.2.0}/physbo/gp/core/__init__.py

@@ -7,5 +7,6 @@
 
 from .prior import prior
 from .model import model
+from .model import sfs
 
 # from predictor import predictor

{physbo-2.1.0 → physbo-2.2.0}/physbo/gp/core/model.py

@@ -11,7 +11,7 @@ from physbo import blm
 from physbo.gp import inf
 from physbo.gp.core import learning
 from physbo.gp.core.prior import prior
-
+from physbo.misc import set_config
 
 class model:
     def __init__(self, lik, mean, cov, inf="exact"):
@@ -310,7 +310,7 @@ class model:
 
         fmean = self.get_post_fmean(X, Z, params=None)
         fcov = self.get_post_fcov(X, Z, params=None, diag=False)
-        return np.random.multivariate_normal(fmean, fcov * alpha**2, N)
+        return np.random.multivariate_normal(fmean, fcov * alpha ** 2, N)
 
     def predict_sampling(self, X, Z, params=None, N=1):
         """
@@ -420,3 +420,81 @@ class model:
             params = bfgs.run(X, t)
 
         self.set_params(params)
+
+class sfs(model):
+
+    def __init__(self, lik, mean, cov, inf="exact",config=None):
+        super().__init__(lik, mean, cov, inf)
+
+        self.config = config
+
+    def prepare(self, X, t, params=None):
+        return super().prepare(X, t, params)
+
+    def fit(self, X, t):
+        """
+        Fitting function (update parameters)
+
+        Parameters
+        ----------
+        X: numpy.ndarray
+            N x d dimensional matrix. Each row of X denotes the d-dimensional feature vector of search candidate.
+
+        t:  numpy.ndarray
+            N dimensional array.
+            The negative energy of each search candidate (value of the objective function to be optimized).
+        config: physbo.misc.set_config object
+
+        """
+        #config = set_config()
+        method = self.config.learning.method
+
+        if method == "adam":
+            adam = learning.adam(self, self.config)
+            params = adam.run(X, t)
+
+        if method in ("bfgs", "batch"):
+            bfgs = learning.batch(self, self.config)
+            params = bfgs.run(X, t)
+
+        self.set_params(params)
+
+        self.prepare(X, t, params=None)
+
+        self.xtrain = X
+    
+    def get_post_fmean(self, X, Z, params=None):
+        return super().get_post_fmean(X, Z, params)
+
+    def predict(self, Z, params=None):
+        """
+        Calculating posterior mean of model (function)
+
+        Parameters
+        ==========
+        X: numpy.ndarray
+            inputs
+        Z: numpy.ndarray
+            feature maps
+        params: numpy.ndarray
+            Parameters
+        See also
+        ========
+        physbo.gp.inf.exact.get_post_fmean
+        """
+
+        if params is None:
+            params = np.copy(self.params)
+
+        if self.inf == "exact":
+            post_fmu = inf.exact.get_post_fmean(self, self.xtrain, Z, params)
+
+        return post_fmu
+    
+    def get_params(self,deep=True):
+
+        mean = self.prior.mean
+        cov = self.prior.cov
+        config = self.config
+
+        return {"lik":self.lik,"mean":mean,"cov":cov,"config":config}

{physbo-2.1.0 → physbo-2.2.0}/physbo/search/pareto.py

@@ -6,7 +6,7 @@
 # file, You can obtain one at https://mozilla.org/MPL/2.0/.
 
 import numpy as np
-
+import time
 
 class Rectangles(object):
     def __init__(self, n_dim, dtype):
@@ -156,7 +156,7 @@ class Pareto(object):
         self.cells.add(lb_idx, ub_idx)
 
     def __included_in_non_dom_region(self, p):
-        return np.all([np.any(pf <= p) for pf in self.front])
+        return np.all(np.any(self.front <= p, axis=1)) # revised on 2025/03/25, vectorized with np.any
 
     def __divide_using_binary_search(self):
         front = np.r_[
@@ -165,17 +165,23 @@ class Pareto(object):
             np.full((1, self.num_objectives), np.inf),
         ]
 
-        # Pareto front indices when sorted on each dimension's front value in ascending order.
-        # (indices start from 1)
-        # Index 0 means anti-ideal value, index `self.front.shape[0] + 1` means ideal point.
+        # Pareto front indices sorted in each dimension
         front_idx = np.r_[
             np.zeros((1, self.num_objectives), dtype=int),
             np.argsort(self.front, axis=0) + 1,
             np.full((1, self.num_objectives), self.front.shape[0] + 1, dtype=int),
         ]
 
+        print("front_idx shape", front_idx.shape)
+
         rect_candidates = [[np.copy(front_idx[0]), np.copy(front_idx[-1])]]
 
+        call_time = 0
+
+        start_time = time.time()
+
+        # modified on 2025/03/25, optimized checking of non-dominated region
+
         while rect_candidates:
             rect = rect_candidates.pop()
 
@@ -184,23 +190,38 @@ class Pareto(object):
             lb = [front[lb_idx[d], d] for d in range(self.num_objectives)]
             ub = [front[ub_idx[d], d] for d in range(self.num_objectives)]
 
+            # first check if lb is in the non-dominated region, if so, add it to the cells
             if self.__included_in_non_dom_region(lb):
                 self.cells.add([lb_idx], [ub_idx])
+                call_time += 1
+                continue
+
+            # when lb is not in the non-dominated region, check if ub is in the non-dominated region
+            # if ub is in the non-dominated region, it means that the whole rectangle is in the dominated region
+            if not self.__included_in_non_dom_region(ub):
+                call_time += 1
+                continue
+
+            # only when lb is not in the non-dominated region and ub is in the non-dominated region,
+            # we need to divide the rectangle into two sub-rectangles
+
+            rect_sizes = rect[1] - rect[0]
+            if np.any(rect_sizes > 1):
+                div_dim = np.argmax(rect_sizes)
+                div_point = rect[0][div_dim] + int(round(rect_sizes[div_dim] / 2.0))
+
+                # Left sub-rectangle
+                left_ub_idx = np.copy(rect[1])
+                left_ub_idx[div_dim] = div_point
+                rect_candidates.append([np.copy(rect[0]), left_ub_idx])
 
-            elif self.__included_in_non_dom_region(ub):
-                rect_sizes = rect[1] - rect[0]
+                # Right sub-rectangle
+                right_lb_idx = np.copy(rect[0])
+                right_lb_idx[div_dim] = div_point
+                rect_candidates.append([right_lb_idx, np.copy(rect[1])])
 
-                # divide rectangle by the dimension with largest size
-                if np.any(rect_sizes > 1):
-                    div_dim = np.argmax(rect_sizes)
-                    div_point = rect[0][div_dim] + int(round(rect_sizes[div_dim] / 2.0))
+        end_time = time.time()
 
-                    # add divided left rectangle
-                    left_ub_idx = np.copy(rect[1])
-                    left_ub_idx[div_dim] = div_point
-                    rect_candidates.append([np.copy(rect[0]), left_ub_idx])
+        print("Execution time:", end_time - start_time)
 
-                    # add divided right rectangle
-                    right_lb_idx = np.copy(rect[0])
-                    right_lb_idx[div_dim] = div_point
-                    rect_candidates.append([right_lb_idx, np.copy(rect[1])])
+        print("call_time", call_time)

{physbo-2.1.0 → physbo-2.2.0}/physbo.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
-Metadata-Version: 2.2
+Metadata-Version: 2.4
 Name: physbo
-Version: 2.1.0
+Version: 2.2.0
 Summary: optimization tool for PHYSics based on Bayesian Optimization
 Home-page: https://github.com/issp-center-dev/PHYSBO
 Author: PHYSBO developers

{physbo-2.1.0 → physbo-2.2.0}/setup.cfg

@@ -1,6 +1,6 @@
 [metadata]
 name = physbo
-version = 2.1.0
+version = 2.2.0
 description = optimization tool for PHYSics based on Bayesian Optimization
 long_description = file: README.md
 long_description_content_type = text/markdown