Method: CArray#lu_refine

Defined in:
lib/carray-gsl/core.rb

#lu_refine(b, x) ⇒ Object 



122
123
124
125
126
127
128
129
 
# File 'lib/carray-gsl/core.rb', line 122

def lu_refine (b, x)
  CArray.attach(self, b, x) {
    a    = self["@A"]
    lu   = self["@LU"]
    perm = self["@perm"]
    return GSL::Linalg::LU.refine(a, lu, perm, b.gv, x.gv).first.ca.v
  }
end