qupled 1.3.2__cp310-cp310-macosx_14_0_arm64.whl

qupled/hf.py

@@ -0,0 +1,263 @@
+from __future__ import annotations
+
+import numpy as np
+
+from . import database
+from . import mpi
+from . import native
+
+
+class HF:
+    """
+    Class used to solve the HF scheme.
+    """
+
+    def __init__(self):
+        self.inputs: Input = None
+        """The inputs used to solve the scheme. Default = ``None``"""
+        self.results: Result = Result()
+        """The results obtained by solving the scheme"""
+        # Undocumented properties
+        self.db_handler = database.DataBaseHandler()
+        self.native_scheme_cls = native.HF
+        self.native_inputs_cls = native.Input
+        self.native_scheme_status = None
+
+    @property
+    def run_id(self):
+        """
+        Property that retrieves the run ID from the database handler.
+
+        Returns:
+            str: The run ID associated with the current database handler.
+        """
+        return self.db_handler.run_id
+
+    @mpi.MPI.record_time
+    @mpi.MPI.synchronize_ranks
+    def compute(self, inputs: Input):
+        """
+        Solves the scheme and saves the results.
+
+        Args:
+            inputs: Input parameters.
+        """
+        self.inputs = inputs
+        self._add_run_to_database()
+        self._compute_native()
+        self._save()
+
+    @mpi.MPI.run_only_on_root
+    def compute_rdf(self, rdf_grid: np.ndarray = None):
+        """
+        Computes the radial distribution function (RDF) using the provided RDF grid.
+        If results are available, this method computes the RDF and stores the results
+        in the database.
+
+        Args:
+            rdf_grid: A numpy array representing the RDF grid.
+                If not provided, a default grid will be used.
+        """
+        if self.results is not None:
+            self.results.compute_rdf(rdf_grid)
+            self.db_handler.insert_results(
+                {"rdf": self.results.rdf, "rdf_grid": self.results.rdf_grid}
+            )
+
+    def _add_run_to_database(self):
+        """
+        Adds the current run information to the database.
+
+        This method inserts the run details stored in `self.inputs` into the database
+        using the `db_handler`. It also updates the `database_info` attribute of
+        `self.inputs` with the current `run_id`.
+        """
+        self.db_handler.insert_run(self.inputs)
+        self.inputs.database_info.run_id = self.run_id
+
+    def _compute_native(self):
+        """
+        Computes the native representation of the inputs and processes the results.
+
+        This method performs the following steps:
+        1. Converts the current inputs to their native representation.
+        2. Initializes a native scheme object using the native inputs.
+        3. Computes the native scheme and stores its status.
+        4. Converts the results from the native scheme back to the desired format.
+        """
+        native_inputs = self.native_inputs_cls()
+        self.inputs.to_native(native_inputs)
+        scheme = self.native_scheme_cls(native_inputs)
+        self.native_scheme_status = scheme.compute()
+        self.results.from_native(scheme)
+
+    @mpi.MPI.run_only_on_root
+    def _save(self):
+        """
+        Saves the current state and results to the database.
+
+        This method updates the run status in the database using the current
+        native scheme status and inserts the results into the database.
+        """
+        self.db_handler.update_run_status(self.native_scheme_status)
+        self.db_handler.insert_results(self.results.__dict__)
+
+
+class Input:
+    """
+    Class used to store the inputs for the :obj:`qupled.hf.HF` class.
+    """
+
+    def __init__(self, coupling: float, degeneracy: float):
+        """
+        Initialize the base class with the given parameters.
+
+        Parameters:
+            coupling (float): Coupling parameter.
+            degeneracy (float): Degeneracy parameter.
+        """
+        self.chemical_potential: list[float] = [-10.0, 10.0]
+        """Initial guess for the chemical potential. Default = ``[-10, 10]``"""
+        self.coupling: float = coupling
+        """Coupling parameter."""
+        self.cutoff: float = 10.0
+        """Cutoff for the wave-vector grid. Default =  ``10.0``"""
+        self.degeneracy: float = degeneracy
+        """Degeneracy parameter."""
+        self.frequency_cutoff: float = 10.0
+        """Cutoff for the frequency (applies only in the ground state). Default =  ``10.0``"""
+        self.integral_error: float = 1.0e-5
+        """Accuracy (relative error) in the computation of integrals. Default = ``1.0e-5``"""
+        self.integral_strategy: str = "full"
+        """
+        Scheme used to solve two-dimensional integrals
+        allowed options include:
+
+        - full: the inner integral is evaluated at arbitrary points selected automatically by the quadrature rule
+
+        - segregated: the inner integral is evaluated on a fixed grid that depends on the integrand that is being processed
+
+        Segregated is usually faster than full but it could become
+        less accurate if the fixed points are not chosen correctly. Default =  ``'full'``
+        """
+        self.matsubara: int = 128
+        """Number of Matsubara frequencies. Default = ``128``"""
+        self.resolution: float = 0.1
+        """Resolution of the wave-vector grid. Default =  ``0.1``"""
+        self.threads: int = 1
+        """Number of OMP threads for parallel calculations. Default =  ``1``"""
+        # Undocumented default values
+        self.theory: str = "HF"
+        self.database_info: DatabaseInfo = DatabaseInfo()
+
+    def to_native(self, native_input: any):
+        """
+        Converts the attributes of the current object to their native representations
+        and sets them on the provided `native_input` object.
+
+        This method iterates through the attributes of the current object and checks
+        if the `native_input` object has a corresponding attribute. If it does, the
+        method attempts to convert the attribute's value to its native representation
+        using a `to_native` method, if available. Otherwise, it directly assigns the
+        attribute's value to the `native_input` object.
+
+        Args:
+            native_input (any): The object to which the native representations of the
+                current object's attributes will be assigned.
+        """
+        name = Input.to_native.__name__
+        for attr, value in self.__dict__.items():
+            if hasattr(native_input, attr) and value is not None:
+                value_to_set = (
+                    tonative()
+                    if callable(tonative := getattr(value, name, None))
+                    else value
+                )
+                setattr(native_input, attr, value_to_set)
+
+
+class Result:
+    """
+    Class used to store the results for the :obj:`qupled.hf.HF` class.
+    """
+
+    def __init__(self):
+        self.idr: np.ndarray = None
+        """Ideal density response"""
+        self.lfc: np.ndarray = None
+        """Local field correction"""
+        self.rdf: np.ndarray = None
+        """Radial distribution function"""
+        self.rdf_grid: np.ndarray = None
+        """Radial distribution function grid"""
+        self.sdr: np.ndarray = None
+        """Static density response"""
+        self.ssf: np.ndarray = None
+        """Static structure factor"""
+        self.uint: float = None
+        """Internal energy"""
+        self.wvg: np.ndarray = None
+        """Wave-vector grid"""
+
+    def from_native(self, native_scheme: any):
+        """
+        Updates the attributes of the current object based on the attributes of a given native scheme object.
+
+        Args:
+            native_scheme (any): An object containing attributes to update the current object with.
+
+        Notes:
+            - Only attributes that exist in both the current object and the native_scheme object will be updated.
+            - Attributes with a value of `None` in the native_scheme object will not overwrite the current object's attributes.
+        """
+        for attr in self.__dict__.keys():
+            if hasattr(native_scheme, attr):
+                value = getattr(native_scheme, attr)
+                setattr(self, attr, value) if value is not None else None
+
+    def compute_rdf(self, rdf_grid: np.ndarray | None = None):
+        """
+        Compute the radial distribution function (RDF) for the system.
+
+        Args:
+            rdf_grid (np.ndarray | None, optional): A 1D array specifying the grid points
+                at which the RDF is computed. If None, a default grid ranging from 0.0
+                to 10.0 with a step size of 0.01 is used.
+
+        Returns:
+            None: The computed RDF is stored in the `self.rdf` attribute.
+        """
+        if self.wvg is not None and self.ssf is not None:
+            self.rdf_grid = (
+                rdf_grid if rdf_grid is not None else np.arange(0.0, 10.0, 0.01)
+            )
+            self.rdf = native.compute_rdf(self.rdf_grid, self.wvg, self.ssf)
+
+
+class DatabaseInfo:
+    """
+    Class used to store the database information passed to the native code.
+    """
+
+    def __init__(self):
+        self.name: str = database.DataBaseHandler.DEFAULT_DATABASE_NAME
+        """Database name"""
+        self.run_id: int = None
+        """ID of the run in the database"""
+        self.run_table_name: str = database.DataBaseHandler.RUN_TABLE_NAME
+        """Name of the table used to store the runs in the database"""
+
+    def to_native(self) -> native.DatabaseInfo:
+        """
+        Converts the current object to a native `DatabaseInfo` instance.
+        This method creates a new instance of `native.DatabaseInfo` and copies
+        all non-None attributes from the current object to the new instance.
+        Returns:
+            native.DatabaseInfo: A new instance of `native.DatabaseInfo` with
+            attributes copied from the current object.
+        """
+        native_database_info = native.DatabaseInfo()
+        for attr, value in self.__dict__.items():
+            if value is not None:
+                setattr(native_database_info, attr, value)
+        return native_database_info

qupled/mpi.py

@@ -0,0 +1,69 @@
+import functools
+
+from qupled import native
+
+
+class MPI:
+    """Class to handle the calls to the MPI API"""
+
+    def rank(self):
+        """Get rank of the process"""
+        return native.MPI.rank()
+
+    def is_root(self):
+        """Check if the current process is root (rank 0)"""
+        return native.MPI.is_root()
+
+    def barrier(self):
+        """Setup an MPI barrier"""
+        native.MPI.barrier()
+
+    def timer(self):
+        """Get wall time"""
+        return native.MPI.timer()
+
+    @staticmethod
+    def run_only_on_root(func):
+        """Python decorator for all methods that have to be run only by root"""
+
+        @functools.wraps(func)
+        def wrapper(*args, **kwargs):
+            if native.MPI.is_root():
+                return func(*args, **kwargs)
+
+        return wrapper
+
+    @staticmethod
+    def synchronize_ranks(func):
+        """Python decorator for all methods that need rank synchronization"""
+
+        @functools.wraps(func)
+        def wrapper(*args, **kwargs):
+            func(*args, **kwargs)
+            native.MPI.barrier()
+
+        return wrapper
+
+    @staticmethod
+    def record_time(func):
+        """Python decorator for all methods that have to be timed"""
+
+        @functools.wraps(func)
+        def wrapper(*args, **kwargs):
+            mpi = native.MPI
+            tic = mpi.timer()
+            func(*args, **kwargs)
+            toc = mpi.timer()
+            dt = toc - tic
+            hours = dt // 3600
+            minutes = (dt % 3600) // 60
+            seconds = dt % 60
+            if mpi.is_root():
+                if hours > 0:
+                    print("Elapsed time: %d h, %d m, %d s." % (hours, minutes, seconds))
+                elif minutes > 0:
+                    print("Elapsed time: %d m, %d s." % (minutes, seconds))
+                else:
+                    print("Elapsed time: %.1f s." % seconds)
+
+        return wrapper

qupled/output.py

@@ -0,0 +1,92 @@
+from . import database
+
+
+class DataBase:
+    """Class to read and delete data from a database generated by qupled"""
+
+    # Read runs in the database
+    @staticmethod
+    def inspect_runs(database_name: str | None = None) -> dict:
+        """Reads runs from the database and returns the content in the form of a dictionary.
+
+        Args:
+            database_name: Name of the database to read from. Defaults to None.
+
+        Returns:
+            A dictionary whose keys are the run ids and values are the corresponding runs information.
+        """
+        db_handler = database.DataBaseHandler(database_name)
+        return db_handler.inspect_runs()
+
+    # Read runs in the database
+    @staticmethod
+    def read_run(
+        run_id: int,
+        database_name: str | None = None,
+        input_names: list[str] | None = None,
+        result_names: list[str] | None = None,
+    ) -> dict:
+        """
+        Reads a run from the database.
+
+        Args:
+            run_id: The ID of the run to read.
+            database_name: The name of the database. Defaults to None.
+            input_names: A list of input names to retrieve. Defaults to None.
+            result_names: A list of result names to retrieve. Defaults to None.
+
+        Returns:
+            dict: A dictionary containing the run data.
+        """
+        db_handler = database.DataBaseHandler(database_name)
+        return db_handler.get_run(run_id, input_names, result_names)
+
+    # Read inputs in the database
+    @staticmethod
+    def read_inputs(
+        run_id: int, database_name: str | None = None, names: list[str] | None = None
+    ) -> dict:
+        """Reads inputs from the database and returns the content in the form of a dictionary.
+
+        Args:
+            run_id: Identifier of the run to read input for.
+            database_name: Name of the database to read from (default is None).
+            names: A list of quantities to read (default is None, which reads all available quantities).
+
+        Returns:
+            A dictionary whose keys are the quantities listed in names and values are the corresponding inputs.
+        """
+        db_handler = database.DataBaseHandler(database_name)
+        return db_handler.get_inputs(run_id, names if names is not None else [])
+
+    # Read results in the database
+    @staticmethod
+    def read_results(
+        run_id: int, database_name: str | None = None, names: list[str] | None = None
+    ) -> dict:
+        """Reads results from the database and returns the content in the form of a dictionary.
+
+        Args:
+            run_id: Identifier of the run to read results for.
+            database_name: Name of the database to read from (default is None).
+            names: A list of quantities to read (default is None, which reads all available quantities).
+
+        Returns:
+            A dictionary whose keys are the quantities listed in names and values are the corresponding results.
+        """
+        db_handler = database.DataBaseHandler(database_name)
+        return db_handler.get_results(run_id, names)
+
+    # Delete results from the database
+    @staticmethod
+    def delete_run(run_id: int, database_name: str | None = None):
+        """
+        Deletes a run entry from the database based on the provided run ID.
+
+        Args:
+            run_id: The unique identifier of the run to be deleted.
+            database_name: The name of the database to connect to.
+                If None, the default database will be used.
+        """
+        db_handler = database.DataBaseHandler(database_name)
+        return db_handler.delete_run(run_id)

qupled/qstls.py

@@ -0,0 +1,68 @@
+from __future__ import annotations
+
+from . import database
+from . import native
+from . import stls
+
+
+class Qstls(stls.Stls):
+    """
+    Class used to solve the Qstls scheme.
+    """
+
+    def __init__(self):
+        super().__init__()
+        self.results: stls.Result = stls.Result()
+        # Undocumented properties
+        self.native_scheme_cls = native.Qstls
+        self.native_inputs_cls = native.QstlsInput
+
+    def compute(self, inputs: Input):
+        self.find_fixed_adr_in_database(inputs)
+        super().compute(inputs)
+
+    def find_fixed_adr_in_database(self, inputs: Input):
+        """
+        Searches the database for a run with matching parameters and assigns its ID to the input object.
+
+        This method iterates through all runs in the database and checks if a run matches the given
+        input parameters (degeneracy, theory, cutoff, matsubara, and resolution). If a match is found,
+        the `fixed_run_id` attribute of the input object is updated with the corresponding run ID.
+
+        Args:
+            inputs (Input): The input parameters.
+
+        Returns:
+            None: The method updates the `fixed_run_id` attribute of the `inputs` object if a match is found.
+        """
+        runs = self.db_handler.inspect_runs()
+        inputs.fixed_run_id = None
+        for run in runs:
+            database_keys = database.DataBaseHandler.TableKeys
+            same_degeneracy = run[database_keys.DEGENERACY.value] == inputs.degeneracy
+            same_theory = run[database_keys.THEORY.value] == inputs.theory
+            if not same_theory or not same_degeneracy:
+                continue
+            run_id = run[database_keys.PRIMARY_KEY.value]
+            run_inputs = self.db_handler.get_inputs(run_id)
+            if (
+                run_inputs["cutoff"] == inputs.cutoff
+                and run_inputs["matsubara"] == inputs.matsubara
+                and run_inputs["resolution"] == inputs.resolution
+            ):
+                print(f"Loading fixed ADR from database for run_id = {run_id}")
+                inputs.fixed_run_id = run_id
+                return
+
+
+# Input class
+class Input(stls.Input):
+    """
+    Class used to manage the input for the :obj:`qupled.qstls.Qstls` class.
+    """
+
+    def __init__(self, coupling: float, degeneracy: float):
+        super().__init__(coupling, degeneracy)
+        # Undocumented default values
+        self.fixed_run_id: int | None = None
+        self.theory = "QSTLS"

qupled/qstlsiet.py

@@ -0,0 +1,37 @@
+from __future__ import annotations
+
+from . import native
+from . import qstls
+from . import stlsiet
+
+
+class QstlsIet(qstls.Qstls):
+    """
+    Class used to solve the Qstls-IET schemes.
+    """
+
+    def __init__(self):
+        super().__init__()
+        self.results: stlsiet.Result = stlsiet.Result()
+        self.native_scheme_cls = native.QstlsIet
+        self.native_inputs_cls = native.QstlsIetInput
+
+    @staticmethod
+    def get_initial_guess(
+        run_id: str, database_name: str | None = None
+    ) -> stlsiet.Guess:
+        return stlsiet.StlsIet.get_initial_guess(run_id, database_name)
+
+
+# Input class
+class Input(stlsiet.Input, qstls.Input):
+    """
+    Class used to manage the input for the :obj:`qupled.qstlsiet.QStlsIet` class.
+    Accepted theories: ``QSTLS-HNC``, ``QSTLS-IOI`` and ``QSTLS-LCT``.
+    """
+
+    def __init__(self, coupling: float, degeneracy: float, theory: str):
+        super().__init__(coupling, degeneracy, "STLS-HNC")
+        if theory not in {"QSTLS-HNC", "QSTLS-IOI", "QSTLS-LCT"}:
+            raise ValueError("Invalid dielectric theory")
+        self.theory = theory

qupled/qvsstls.py

@@ -0,0 +1,59 @@
+from __future__ import annotations
+
+from . import native
+from . import qstls
+from . import vsstls
+
+
+class QVSStls(vsstls.VSStls):
+    """
+    Class used to solve the QVStls scheme.
+    """
+
+    def __init__(self):
+        super().__init__()
+        self.results: vsstls.Result = vsstls.Result()
+        # Undocumented properties
+        self.native_scheme_cls = native.QVSStls
+        self.native_inputs_cls = native.QVSStlsInput
+
+    def compute(self, inputs: Input):
+        """
+        Solves the scheme and saves the results.
+
+        Args:
+            inputs: Input parameters.
+        """
+        qstls.Qstls.find_fixed_adr_in_database(self, inputs)
+        super().compute(inputs)
+
+    def _update_input_data(self, inputs: Input):
+        """
+        Updates the input data with additional attributes specific to the current instance.
+
+        This method overrides the parent class's `_update_input_data` method to include
+        logic for setting a default `fixed_run_id` if it is not already provided in the
+        `inputs`.
+
+        Args:
+            inputs: Input parameters.
+
+        Side Effects:
+            - If `inputs.fixed_run_id` is `None`, it is set to the value of `self.run_id`.
+        """
+        super()._update_input_data(inputs)
+        if inputs.fixed_run_id is None:
+            inputs.fixed_run_id = self.run_id
+
+
+# Input class
+class Input(vsstls.Input, qstls.Input):
+    """
+    Class used to manage the input for the :obj:`qupled.qvsstls.QVSStls` class.
+    """
+
+    def __init__(self, coupling: float, degeneracy: float):
+        vsstls.Input.__init__(self, coupling, degeneracy)
+        qstls.Input.__init__(self, coupling, degeneracy)
+        # Undocumented default values
+        self.theory: str = "QVSSTLS"

qupled/rpa.py

@@ -0,0 +1,27 @@
+from __future__ import annotations
+
+from . import hf
+from . import native
+
+
+class Rpa(hf.HF):
+    """
+    Class used to solve the RPA scheme.
+    """
+
+    def __init__(self):
+        super().__init__()
+        self.results: hf.Result = hf.Result()
+        # Undocumented properties
+        self.native_scheme_cls = native.Rpa
+
+
+class Input(hf.Input):
+    """
+    Class used to manage the input for the :obj:`qupled.rpa.Rpa` class.
+    """
+
+    def __init__(self, coupling: float, degeneracy: float):
+        super().__init__(coupling, degeneracy)
+        # Undocumented default values
+        self.theory = "RPA"