Method: Stick::Matrix::Jacobi.max

Defined in:
lib/stick/matrix/jacobi.rb

.max(a) ⇒ Object

Returns the index pair (p, q) with 1<= p < q <= n and A[p, q] is the maximum in absolute value



18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
 
# File 'lib/stick/matrix/jacobi.rb', line 18

def Jacobi.max(a)
  n = a.row_size
  max = 0
  p = 0
  q = 0
  n.times{|i|
    ((i+1)...n).each{|j|
      val = a[i, j].abs
      if val > max
        max = val
        p = i
        q = j
      end  }}
  return p, q
end