#include <triqs/det_manip/det_manip.hpp>


det_manip::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().