private void CalculateSiteMagnetic(AtoOCalculator aToO, AnomalyCurrent jQ)
{
foreach (var site in _observationSites)
{
var normalField = CalculateNormalFieldH(site);
// var anomalyField = CalculateAnomalyFieldH(site, jQ);
// OnHFieldsAtSiteCalculated(site, normalField, anomalyField);
}
}
protected void SetModel(OmegaModel model, TensorCache tensors = null)
{
if (IsParallel)
{
Mpi.CheckNumberOfProcesses(model.LateralDimensions);
}
Model = model;
if (_aToOCalculator == null)
{
_aToOCalculator = new AtoOCalculator(this);
}
if (_fgmresSolver == null)
{
_fgmresSolver = new AnomalyCurrentFgmresSolver(this);
}
if (_convolutionOperator == null)
{
_convolutionOperator = new ConvolutionOperator(this);
}
if (tensors == null)
{
_aToOCalculator.CleanGreenTensors();
if (!StoreAtoA)
{
_greenTensorAtoA?.Dispose();
_greenTensorAtoA = null;
}
}
else
{
_aToOCalculator.SetTensors(tensors.eGreenTensors, tensors.hGreenTensors);
SetNewGreenTensor(tensors.gtAtoA);
}
}