if (a%is_triangle()) then if (a%is_upper()) then if (a%is_lower()) then if (a%is_unit()) then
logical
value indicating whether the
matrix is triangular; if is_triangle()
returns .true.
check also if it is lower, upper and with a unit (i.e. assumed)
diagonal.