WMD.h

//---------------------------------------------------------------------------
//
// Project: OpenWalnut ( http://www.openwalnut.org )
//
// Copyright 2015 OpenWalnut Community, Nemtics, BSV@Uni-Leipzig
// For more information see http://www.openwalnut.org/copying
//
// This file is part of OpenWalnut.
//
// OpenWalnut is free software: you can redistribute it and/or modify
// it under the terms of the GNU Lesser General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// OpenWalnut is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with OpenWalnut. If not, see <http://www.gnu.org/licenses/>.
//
//---------------------------------------------------------------------------
 
#ifndef WMD_H
#define WMD_H
 
#include "core/common/math/linearAlgebra/WVectorFixed.h"
#include "core/common/math/WTensorSym.h"
#include "core/common/WObjectNDIP.h"
 
#include "WDataSetDTIToScalar_I.h"
 
/**
 * Computes the MD (Mean Diffusivity) of a given Tensor.
 * \f[
 *   MD = \frac{1}{3}\sum_{i=1}^{3} \lambda_i
 * \f]
 */
class WMD : public WObjectNDIP< WDataSetDTIToScalar_I >
{
public:
    /**
     * Creates an object to perform the computation.
     */
    WMD();
 
protected:
    /**
     * Actual MD computation takes place inhere.
     *
     * \param evals With the three given Eigenvalues, we may compute the MD. See the wikipedia article above for the formula.
     * \param tensor Although not needed for MD computation, the API requires us to use the signature.
     *
     * \return The MD of the tensor.
     */
    virtual double tensorToScalar( const WVector3d& evals, const WTensorSym< 2, 3, float >& tensor );
};
 
#endif  // WMD_H