tspacefigofmerit.h

Go to the documentation of this file.

/* This file is part of the Tomographer project, which is distributed under the
 * terms of the MIT license.
 *
 * The MIT License (MIT)
 *
 * Copyright (c) 2016 ETH Zurich, Institute for Theoretical Physics, Philippe Faist
 *
 * Permission is hereby granted, free of charge, to any person obtaining a copy
 * of this software and associated documentation files (the "Software"), to deal
 * in the Software without restriction, including without limitation the rights
 * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
 * copies of the Software, and to permit persons to whom the Software is
 * furnished to do so, subject to the following conditions:
 *
 * The above copyright notice and this permission notice shall be included in
 * all copies or substantial portions of the Software.
 *
 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
 * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
 * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
 * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
 * SOFTWARE.
 */


#ifndef TOMOGRAPHER_DENSEDM_TSPACEFIGOFMERIT_H
#define TOMOGRAPHER_DENSEDM_TSPACEFIGOFMERIT_H


#include <tomographer2/tools/needownoperatornew.h>
#include <tomographer2/densedm/dmtypes.h>
#include <tomographer2/densedm/distmeasures.h>
#include <tomographer2/densedm/param_herm_x.h>


namespace Tomographer {
namespace DenseDM {
namespace TSpace {


template<typename DMTypes_, typename ValueType_ = double>
class FidelityToRefCalculator
  : public virtual Tools::NeedOwnOperatorNew<typename DMTypes_::MatrixType>::ProviderType
{
public:
  typedef DMTypes_ DMTypes;
  typedef typename DMTypes::MatrixType MatrixType;
  typedef typename DMTypes::MatrixTypeConstRef MatrixTypeConstRef;

  typedef ValueType_ ValueType;

private:
  MatrixType _ref_T;

public:
  FidelityToRefCalculator(MatrixTypeConstRef T_ref)
    : _ref_T(T_ref)
  {
  }

  inline ValueType getValue(MatrixTypeConstRef T) const
  {
    return fidelityT<ValueType>(T, _ref_T);
  }
};




template<typename DMTypes_, typename ValueType_ = double>
class PurifDistToRefCalculator
  : public virtual Tools::NeedOwnOperatorNew<typename DMTypes_::MatrixType>::ProviderType
{
public:
  typedef DMTypes_ DMTypes;
  typedef typename DMTypes::MatrixType MatrixType;
  typedef typename DMTypes::MatrixTypeConstRef MatrixTypeConstRef;

  typedef ValueType_ ValueType;

private:
  MatrixType _ref_T;

public:
  PurifDistToRefCalculator(MatrixTypeConstRef T_ref)
    : _ref_T(T_ref)
  {
  }

  inline ValueType getValue(MatrixTypeConstRef T) const
  {
    ValueType F = fidelityT<ValueType>(T, _ref_T);
    if (F >= ValueType(1)) {
      return 0;
    }
    return std::sqrt(ValueType(1) - F*F);
  }
};

template<typename DMTypes_, typename ValueType_ = double>
class TrDistToRefCalculator
  : public virtual Tools::NeedOwnOperatorNew<typename DMTypes_::MatrixType>::ProviderType
{
public:
  typedef DMTypes_ DMTypes;
  typedef typename DMTypes::MatrixType MatrixType;
  typedef typename DMTypes::MatrixTypeConstRef MatrixTypeConstRef;

  typedef ValueType_ ValueType;

private:
  MatrixType _ref_rho;

public:
  TrDistToRefCalculator(MatrixTypeConstRef rho_ref)
    : _ref_rho(rho_ref)
  {
  }

  inline ValueType getValue(MatrixTypeConstRef T) const
  {
    return traceDistance<ValueType>(T*T.adjoint(), _ref_rho);
  }
};



template<typename DMTypes_>
class ObservableValueCalculator
  : public virtual Tools::NeedOwnOperatorNew<typename DMTypes_::VectorParamType>::ProviderType
{
public:
  typedef DMTypes_ DMTypes;
  typedef typename DMTypes::MatrixType MatrixType;
  typedef typename DMTypes::MatrixTypeConstRef MatrixTypeConstRef;
  typedef typename DMTypes::VectorParamType VectorParamType;
  typedef typename DMTypes::VectorParamTypeConstRef VectorParamTypeConstRef;

  typedef typename DMTypes::RealScalar ValueType;

private:
  ParamX<DMTypes> _param_x;

  VectorParamType _A_x;

public:
  ObservableValueCalculator(DMTypes dmt, MatrixTypeConstRef A)
    : _param_x(dmt), _A_x(_param_x.HermToX(A))
  {
  }
  ObservableValueCalculator(DMTypes dmt, VectorParamTypeConstRef A_x)
    : _param_x(dmt), _A_x(A_x)
  {
  }

  inline ValueType getValue(MatrixTypeConstRef T) const
  {
    return _A_x.transpose() * _param_x.HermToX(T*T.adjoint());
  }
};


} // namespace TSpace
} // namespace DenseDM
} // namespace Tomographer


#endif