analysis.MSTE package¶
Submodules¶
analysis.MSTE.MSTE module¶
MSTE calculation
-
analysis.MSTE.MSTE.cluster(x, v, t, box, energy)[source]¶ Get either MST or MSTE clusters from the configuration.
Parameters: x : numpy float64 array
Positions of the particles in the system
v : numpy float64 array
Velocities of the particles
t : numpy int32 array
Types of the particles
box : numpy float64 array
Box
energy : boolean
Whether to do energy considerations in the cluster computation.
Returns: index : numpy array
An array with the cluster index of each particle.
-
analysis.MSTE.MSTE.connections(x, v, t, index, box, energy, expansion)[source]¶ Find the connections with a given cluster distribution.
Parameters: x : numpy float64 array
Positions of the particles in the system
v : numpy float64 array
Velocities of the particles
t : numpy int32 array
Types of the particles
index : numpy array
Cluster index of the particles
box : numpy float64 array
Box
energy : boolean
Whether to do energy considerations in the cluster computation.
expansion : float, optional
The expansion velocity of the walls of the box, in box units.
Returns: conn : 2D numpy array
The connections between the clusters, in the format: [[cluster1, cluster2, wall], [cluster1, cluster2, wall]...]
Notes
- So far, when energy is “True”, this only works for systems with medium potential.
- The box needs to be cubic for this to work
-
analysis.MSTE.MSTE.mste(x, v, t, box, energy, expansion=0.0)[source]¶ Calculate MSTE.
Parameters: x : numpy float64 array
Positions of the particles in the system
v : numpy float64 array
Velocities of the particles
t : numpy int32 array
Types of the particles
box : numpy float64 array
Box
energy : boolean
Whether to do energy considerations in the cluster computation.
expansion : float, optional
The expansion velocity of the walls of the box, in box units.
Returns: value, (mst, inf) : numpy array, numpy array, list
value is the [mass, occupancy, fraction] histogram mst is the array of indices to which each particle belongs. inf is the list of infinite clusters.
Notes
1. So far, when energy is “True”, this only works for systems with medium potential.
- The box needs to be cubic for this to work