Method StoreData

kratos/tensor_adaptors/variable_tensor_adaptor.cpp:155–203

}

void VariableTensorAdaptor::StoreData()
{
    std::visit([this](auto pContainer, auto pVariable) {
        using container_type = BareType<decltype(*pContainer)>;

        if constexpr(IsInList<container_type,
                              ModelPart::NodesContainerType,
                              ModelPart::ConditionsContainerType,
                              ModelPart::ElementsContainerType,
                              ModelPart::PropertiesContainerType,
                              ModelPart::MasterSlaveConstraintContainerType,
                              ModelPart::GeometryContainerType>) {

            using variable_type = BareType<decltype(*pVariable)>;
            using data_type = typename variable_type::Type;

            const auto& r_tensor_shape = this->Shape();

            KRATOS_ERROR_IF_NOT(r_tensor_shape[0] == pContainer->size())
                << "Underlying container of the tensor data has changed size [ tensor data = "
                << this->Info() << ", container size = " << pContainer->size() << " ].\n";

            if constexpr(DataTypeTraits<data_type>::IsDynamic) {
                // this zero value may be different from the Variable::Zero()
                // in the case where the variable type is Variable<Vector> or Variable<Matrix>.
                // Because, in these dynamic data type variables, Variable::Zero will create
                // a zero sized vector or matrix, which is useless in the StoreData method.
                // The following method creates correctly sized zero Vector or Matrix
                // accordingly to be assigned to the entities, if they don't have the specified
                // variable in their DataValueContainer.
                const auto& zero = TensorAdaptorUtils::GetZeroValue(*pVariable, this->DataShape());

                ContainerIOUtils::CopyFromContiguousDataArray<data_type>(
                    *pContainer, this->ViewData(), r_tensor_shape.data().begin(),
                    r_tensor_shape.data().begin() + r_tensor_shape.size(),
                    [pVariable, &zero](auto& rEntity) -> auto& {
                        return rEntity.Emplace(*pVariable, zero);
                    });
            } else {
                ContainerIOUtils::CopyFromContiguousDataArray<data_type>(
                    *pContainer, this->ViewData(), r_tensor_shape.data().begin(),
                    r_tensor_shape.data().begin() + r_tensor_shape.size(),
                    [pVariable](auto& rEntity) -> auto& {
                        return rEntity.Emplace(*pVariable);
                    });
            }
        }
    }, this->GetContainer(), mpVariable);
}

std::string VariableTensorAdaptor::Info() const
{