Source code for pina.model.layers.pod
"""Module for Base Continuous Convolution class."""
from abc import ABCMeta, abstractmethod
import torch
from .stride import Stride
from .utils_convolution import optimizing
[docs]
class PODBlock(torch.nn.Module):
"""
POD layer: it projects the input field on the proper orthogonal
decomposition basis. It needs to be fitted to the data before being used
with the method :meth:`fit`, which invokes the singular value decomposition.
The layer is not trainable.
.. note::
All the POD modes are stored in memory, avoiding to recompute them when the rank changes but increasing the memory usage.
"""
def __init__(self, rank, scale_coefficients=True):
"""
Build the POD layer with the given rank.
:param int rank: The rank of the POD layer.
:param bool scale_coefficients: If True, the coefficients are scaled
after the projection to have zero mean and unit variance.
"""
super().__init__()
self.__scale_coefficients = scale_coefficients
self._basis = None
self._scaler = None
self._rank = rank
@property
def rank(self):
"""
The rank of the POD layer.
:rtype: int
"""
return self._rank
@rank.setter
def rank(self, value):
if value < 1 or not isinstance(value, int):
raise ValueError("The rank must be positive integer")
self._rank = value
@property
def basis(self):
"""
The POD basis. It is a matrix whose columns are the first `self.rank` POD modes.
:rtype: torch.Tensor
"""
if self._basis is None:
return None
return self._basis[: self.rank]
@property
def scaler(self):
"""
The scaler. It is a dictionary with the keys `'mean'` and `'std'` that
store the mean and the standard deviation of the coefficients.
:rtype: dict
"""
if self._scaler is None:
return
return {
"mean": self._scaler["mean"][: self.rank],
"std": self._scaler["std"][: self.rank],
}
@property
def scale_coefficients(self):
"""
If True, the coefficients are scaled after the projection to have zero
mean and unit variance.
:rtype: bool
"""
return self.__scale_coefficients
[docs]
def fit(self, X):
"""
Set the POD basis by performing the singular value decomposition of the
given tensor. If `self.scale_coefficients` is True, the coefficients
are scaled after the projection to have zero mean and unit variance.
:param torch.Tensor X: The tensor to be reduced.
"""
self._fit_pod(X)
if self.__scale_coefficients:
self._fit_scaler(torch.matmul(self._basis, X.T))
def _fit_scaler(self, coeffs):
"""
Private merhod that computes the mean and the standard deviation of the
given coefficients, allowing to scale them to have zero mean and unit
variance. Mean and standard deviation are stored in the private member
`_scaler`.
:param torch.Tensor coeffs: The coefficients to be scaled.
"""
self._scaler = {
"std": torch.std(coeffs, dim=1),
"mean": torch.mean(coeffs, dim=1),
}
def _fit_pod(self, X):
"""
Private method that computes the POD basis of the given tensor and stores it in the private member `_basis`.
:param torch.Tensor X: The tensor to be reduced.
"""
if X.device.type == "mps": # svd_lowrank not arailable for mps
self._basis = torch.svd(X.T)[0].T
else:
self._basis = torch.svd_lowrank(X.T, q=X.shape[0])[0].T
[docs]
def forward(self, X):
"""
The forward pass of the POD layer. By default it executes the
:meth:`reduce` method, reducing the input tensor to its POD
representation. The POD layer needs to be fitted before being used.
:param torch.Tensor X: The input tensor to be reduced.
:return: The reduced tensor.
:rtype: torch.Tensor
"""
return self.reduce(X)
[docs]
def reduce(self, X):
"""
Reduce the input tensor to its POD representation. The POD layer needs
to be fitted before being used.
:param torch.Tensor X: The input tensor to be reduced.
:return: The reduced tensor.
:rtype: torch.Tensor
"""
if self._basis is None:
raise RuntimeError(
"The POD layer needs to be fitted before being used."
)
coeff = torch.matmul(self.basis, X.T)
if coeff.ndim == 1:
coeff = coeff.unsqueeze(1)
coeff = coeff.T
if self.__scale_coefficients:
coeff = (coeff - self.scaler["mean"]) / self.scaler["std"]
return coeff
[docs]
def expand(self, coeff):
"""
Expand the given coefficients to the original space. The POD layer needs
to be fitted before being used.
:param torch.Tensor coeff: The coefficients to be expanded.
:return: The expanded tensor.
:rtype: torch.Tensor
"""
if self._basis is None:
raise RuntimeError(
"The POD layer needs to be trained before being used."
)
if self.__scale_coefficients:
coeff = coeff * self.scaler["std"] + self.scaler["mean"]
predicted = torch.matmul(self.basis.T, coeff.T).T
if predicted.ndim == 1:
predicted = predicted.unsqueeze(0)
return predicted