Sub_Matrix_Decomposition.h

Go to the documentation of this file.

/*
Created on Fri Jun 26 14:13:26 2020
Copyright (C) 2020 Peter Rakyta, Ph.D.

This program is free software: you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.

This program 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 General Public License for more details.

You should have received a copy of the GNU General Public License
along with this program.  If not, see http://www.gnu.org/licenses/.

@author: Peter Rakyta, Ph.D.
*/
#ifndef SUB_MATRIX_DECOMPOSITION_H
#define SUB_MATRIX_DECOMPOSITION_H

#include "qgd/Decomposition_Base.h"
#include "qgd/Sub_Matrix_Decomposition_Cost_Function.h"
#include "qgd/Functor_Cost_Function_Gradient.h"


class Sub_Matrix_Decomposition : public Decomposition_Base {


public:

    bool subdisentaglement_done;

    Matrix subdecomposed_mtx;

    bool optimize_layer_num;

    std::map<int,int> identical_blocks;


public:

Sub_Matrix_Decomposition( Matrix Umtx_in, int qbit_num_in, bool optimize_layer_num_in, guess_type initial_guess_in );

~Sub_Matrix_Decomposition();


void disentangle_submatrices();


void solve_layer_optimization_problem( int num_of_parameters, gsl_vector *solution_guess_gsl);




double optimization_problem( const double* parameters);



static double optimization_problem( const gsl_vector* parameters, void* void_instance );


static void optimization_problem_grad( const gsl_vector* parameters, void* void_instance, gsl_vector* grad  );


static void optimization_problem_combined( const gsl_vector* parameters, void* void_instance, double* f0, gsl_vector* grad );

int set_identical_blocks( int qbit, int identical_blocks_in );


int set_identical_blocks( std::map<int, int> identical_blocks_in );


Sub_Matrix_Decomposition* clone();


};


#endif //SUB_MATRIX_DECOMPOSITION