py4sci

KMCLib 1.1 documentation

Plugin Interfaces

«  KMCLib API   ::   Contents   ::   Analysis  »

Plugin Interfaces

KMCLib can easily be extended through the use of plugins, where the user can write extension to the code and get them easily integrated with the core functionality simply by deriving from a specified plugin base class and overloading certain functions on the base class.

The plugins will always have a Python interface class to inherit from, also in the case when the actual code will be called from within the C++ core routines. This allows the user to write plugin extensions in Python without having to recompile the underlying C++ code.

Calls from C++ to Python are managed through SWIG Director classes. These calls inevitably come with an overhead compared to calling compiled code within C++ (see the SWIG documentation for details). The overhead should fortunately in most realistic cases be negligible. On the other hand what happens inside the custom Python functions could very well be time consuming enough to require special attention and implementation in a compiled language. Such techniques will however not be covered here.

KMCRateCalculatorPlugin

The KMCRateCalculatorPluin is used to customize the calculation of rates in the KMC simulation. A user defined RateCalculator plugin must inherit from the KMCRateCalculatorPlugin base class as in this example:

from KMCLib.KMCRateCalculatorPlugin import KMCRateCalculatorPlugin

# Define the derived class.
class CustomPlugin(KMCRateCalculatorPlugin):

    def rate(self, ... ):
        # Implement the custom rate routine.
        ...

    # Override the other base class API functions if needed.
    ...

The derived class is then given to the KMCInteractions object to be used in the KMC simulation.

# Construct the interactions object.
interactions = KMCInteractions( ... )

# Set the custom rate calculator.
interactions.setRateCalculator(rate_calculator=CustomPlugin)

# Construct the KMCLatticeModel.
model = KMCLatticeModel(interactions=interactions, ... )

When running a KMC simulation with a custom rate calculator defined, its rate method will be called each time the rate of an available process is needed. The rate method receives information about the geometry and the types before and after the process, as well as the unmodified rate constant associated with the process. It is then up to the user to handle the geometric and types information to compute a new rate based on the given information. The rate method should return a valid rate for the particular process for the given geometry.

Note that the rate method of the custom rate calculator plugin could be called millions of times during a typical KMC simulation production run. It is therefore highly recommendable that code efficiency is kept in mind when implementing this type of custom extension.

Interface

class KMCLib.KMCRateCalculatorPlugin(configuration)

Class for providing an interface to easily extend and customize the behaviour of the calculation of individual rates in the KMC simulation.

initialize()

Called as the last statement in the base class constructor to allow for custom setup of the object.

rate(coords, types_before, types_after, rate_constant, process_number, global_coordinate)

Called from the base class to get the rate for a particular local geometry. Any class inheriting from the plugin base class must provide an implementation of this function.

Parameters:
  • coords – The coordinates of the configuration as a Nx3 numpy array in fractional units of the primitive cell.
  • types_before – The types before the process, as tuple of strings.
  • types_after – The types after the process, as tuple of strings.
  • rate_constant – The rate constant associated with the process to either update or replace.
  • process_number – The process id number.
  • global_coordinate – The global coordinate of the central index.
Returns:

The custom rate of the process. Note that the returned rate must not be negative or zero.

cutoff()

To determine the radial cutoff of the geometry around the central lattice site to cut out and send down to the rustom rate function. If not implemented by derrived classes the default is to use the cutoff of the largetst process local geometry.

Returns:The desiered cutoff in primitive cell internal coordinates.
Return type:float

KMCAnalysisPlugin

The KMCAnalysisPlugin is used to perform on-the-fly custom analysis during the KMC simulation. A list of instantiated KMCAnalysisPlugin objects can be given to the run function of the KMCLatticeModel. Each user-defined custom analysis object must inherit from the KMCAnalysisPlugin base class. Three functions are available to overload for setting up, handling updates after a step and finalizing the analysis. Here is an example:

from KMCLib.KMCAnalysisPlugin import KMCAnalysisPlugin

# Define the custom analysis class.
class CustomAnalysis(KMCAnalysisPlugin)
    """ Custom analysis class """
    def __init__(self, ... ):
        # Initialize the object.
        ...

    # Overload the interface functions.
    def setup(self, step, time, configuration):
        ...

    def registerStep(self, step, time, configuration):
        ...

    def finalize(self):
        ...

# Instantiate the custom analysis object.
my_analysis1 = CustomAnalysis(...)
my_analysis2 = CustomAnalysis(...)
...

# Setup the list of analysis objects.
analysis = [my_analysis1, my_analysis2, ... ]

# Setup the kmc lattice model to run (see documentation for parameters)
kmc_model = KMCLatticeMode(...)

# Start the simulation with custom on-the-fly analysis every 10:th step.
kmc_mode.run( ...
              control_parameters=KMCControlParameters( ...
                                                       analysis_interval=10)
              analysis=analysis)

Note that how often the registerStep function is called during the simulation is controlled via the analysis_interval parameter given to the KMCControlParameters.

Interface

class KMCLib.KMCAnalysisPlugin

Class for providing an interface to easily extend and customize the behaviour of the on-the-fly analysis functionality of KMCLib.

setup(step, time, configuration)

Function called right before the start of the KMC loop to allow for custom setup of the analysis object.

Parameters:
  • step (int) – The simulation step number.
  • time (float) – The simulation time.
  • configuration (KMCConfiguration) – The up to date configuration of the simulation.
registerStep(step, time, configuration)

Called from the KMC loop after each step.

Parameters:
  • step (int) – The simulation step number.
  • time (float) – The simulation time.
  • configuration (KMCConfiguration) – The up to date configuration of the simulation.
finalize()

Called after the KMC loop to allow for custom finalization and post-processing of collected data.

__init__()

The constructor of the base-class.

«  KMCLib API   ::   Contents   ::   Analysis  »