method

unitary?

v1_9_3_125 - Show latest stable - Class: Matrix

unitary?()

public

Returns true is this is a unitary matrix Raises an error if matrix is not square.

# File lib/matrix.rb, line 729
  def unitary?
    Matrix.Raise ErrDimensionMismatch unless square?
    rows.each_with_index do |row, i|
      column_size.times do |j|
        s = 0
        row_size.times do |k|
          s += row[k].conj * rows[k][j]
        end
        return false unless s == (i == j ? 1 : 0)
      end
    end
    true
  end

unitary?

Related methods