triqs::det_manip::det_manip::try_insert_k
#include <triqs/det_manip/det_manip.hpp>
Synopsis
value_type try_insert_k (std::vector<long> i,std::vector<long> j,std::vector<x_type> x,std::vector<y_type> y)
Double Insert operation at colum j0,j1 and row i0,i1.
- The operation consists in adding :
two columns f(x_i, y_{j0}), f(x_i, y_{j1})
and two rows f(x_{i0}, x_j), f(x_{i1}, x_j)
The new colums/rows will be at col j0, j1, row i0, i1.
0 <= i0,i1,j0,j1 <= N+1, where N is the current size of the matrix. Returns the ratio of det Minv_new / det Minv. This routine does NOT make any modification. It has to be completed with complete_operation().