-// Copyright (C) 2007-2014 CEA/DEN, EDF R&D
+// Copyright (C) 2007-2015 CEA/DEN, EDF R&D
//
// This library is free software; you can redistribute it and/or
// modify it under the terms of the GNU Lesser General Public
#include <sstream>
-using namespace ParaMEDMEM;
+using namespace MEDCoupling;
bool MEDCouplingField::isEqualIfNotWhy(const MEDCouplingField *other, double meshPrec, double valsPrec, std::string& reason) const
{
return _mesh==other->_mesh;
}
+/*!
+ * This method is less strict than MEDCouplingField::areStrictlyCompatible method.
+ * The difference is that the nature is not checked.
+ * This method is used for multiplication and division on fields to operate a first check before attempting operation.
+ */
+bool MEDCouplingField::areStrictlyCompatibleForMulDiv(const MEDCouplingField *other) const
+{
+ if(!other)
+ throw INTERP_KERNEL::Exception("MEDCouplingField::areStrictlyCompatible : input field is NULL !");
+ if(!_type->isEqual(other->_type,1.e-12))
+ return false;
+ return _mesh==other->_mesh;
+}
+
+
void MEDCouplingField::updateTime() const
{
if(_mesh)
return ret;
}
-std::vector<const BigMemoryObject *> MEDCouplingField::getDirectChildren() const
+std::vector<const BigMemoryObject *> MEDCouplingField::getDirectChildrenWithNull() const
{
std::vector<const BigMemoryObject *> ret;
- if(_mesh)
- ret.push_back(_mesh);
- if((const MEDCouplingFieldDiscretization *)_type)
- ret.push_back((const MEDCouplingFieldDiscretization *)_type);
+ ret.push_back(_mesh);
+ ret.push_back((const MEDCouplingFieldDiscretization *)_type);
return ret;
}
/*!
* Returns a type of \ref MEDCouplingSpatialDisc "spatial discretization" of \a this
- * field in terms of enum ParaMEDMEM::TypeOfField.
- * \return ParaMEDMEM::TypeOfField - the type of \a this field.
+ * field in terms of enum MEDCoupling::TypeOfField.
+ * \return MEDCoupling::TypeOfField - the type of \a this field.
* \throw If the geometric type is empty.
*/
TypeOfField MEDCouplingField::getTypeOfField() const
/*!
* Returns the nature of \a this field. This information is very important during
- * interpolation process using ParaMEDMEM::MEDCouplingRemapper or ParaMEDMEM::InterpKernelDEC.
+ * interpolation process using MEDCoupling::MEDCouplingRemapper or MEDCoupling::InterpKernelDEC.
* In other context than the two mentioned above, this attribute is unimportant. This
* attribute is not stored in the MED file.
* For more information of the semantics and the influence of this attribute to the
/*!
* Sets the nature of \a this field. This information is very important during
- * interpolation process using ParaMEDMEM::MEDCouplingRemapper or ParaMEDMEM::InterpKernelDEC.
+ * interpolation process using MEDCoupling::MEDCouplingRemapper or MEDCoupling::InterpKernelDEC.
* In other context than the two mentioned above, this attribute is unimportant. This
* attribute is not stored in the MED file.
* For more information of the semantics and the influence of this attribute to the
* \throw If size of any vector do not match the \a type.
*/
void MEDCouplingField::setGaussLocalizationOnType(INTERP_KERNEL::NormalizedCellType type, const std::vector<double>& refCoo,
- const std::vector<double>& gsCoo, const std::vector<double>& wg) throw(INTERP_KERNEL::Exception)
+ const std::vector<double>& gsCoo, const std::vector<double>& wg)
{
if(!_mesh)
throw INTERP_KERNEL::Exception("Mesh has to be set before calling setGaussLocalizationOnType method !");
* \throw If the range [_begin_,_end_) is empty.
*/
void MEDCouplingField::setGaussLocalizationOnCells(const int *begin, const int *end, const std::vector<double>& refCoo,
- const std::vector<double>& gsCoo, const std::vector<double>& wg) throw(INTERP_KERNEL::Exception)
+ const std::vector<double>& gsCoo, const std::vector<double>& wg)
{
if(!_mesh)
throw INTERP_KERNEL::Exception("Mesh has to be set before calling setGaussLocalizationOnCells method !");
* problem.
* \param [in] locId - the id of the Gauss localization object of interest.
* It must be in range <em> 0 <= locId < getNbOfGaussLocalization() </em>.
- * \return \ref MEDCouplingGaussLocalization & - the Gauss localization object.
+ * \return \ref MEDCoupling::MEDCouplingGaussLocalization "MEDCouplingGaussLocalization" & - the
+ * Gauss localization object.
* \throw If \a this field is not on Gauss points.
* \throw If \a locId is not within the valid range.
* \throw If the spatial discretization of \a this field is NULL.
* \warning This method is const, so the returned object is not apt for modification.
* \param [in] locId - the id of the Gauss localization object of interest.
* It must be in range <em> 0 <= locId < getNbOfGaussLocalization() </em>.
- * \return \ref const MEDCouplingGaussLocalization & - the Gauss localization object.
+ * \return const \ref MEDCouplingGaussLocalization & - the Gauss localization object.
* \throw If \a this field is not on Gauss points.
* \throw If \a locId is not within the valid range.
* \throw If the spatial discretization of \a this field is NULL.
}
MEDCouplingField::MEDCouplingField(const MEDCouplingField& other, bool deepCopy):RefCountObject(other),_name(other._name),_desc(other._desc),_nature(other._nature),
- _mesh(0),_type(0)
+ _mesh(0),_type(0)
{
if(other._mesh)
{
/*!
* Returns a new MEDCouplingMesh constituted by some cells of the underlying mesh of \a
- * this filed, and returns ids of entities (nodes, cells, Gauss points) lying on the
+ * this field, and returns ids of entities (nodes, cells, Gauss points) lying on the
* specified cells. The cells to include to the result mesh are specified by an array of
* cell ids. The new mesh shares the coordinates array with the underlying mesh.
* \param [in] start - an array of cell ids to include to the result mesh.
/*!
* Returns number of tuples expected regarding the spatial discretization of \a this
- * field and number of entities in the underlying mesh.
+ * field and number of entities in the underlying mesh. This method behaves exactly as MEDCouplingFieldDouble::getNumberOfTuples.
* \return int - the number of expected tuples.
* \throw If the spatial discretization of \a this field is NULL.
* \throw If the mesh is not set.
