src/ij2r.cpp
int ij2rUpTr(int Nel, int i, int j)
{
  if ((i>j)||(i<0)||(j<0)) return(0);
  else
    return(i*Nel-i*(i-1)/2+(j-i));
}
////////////////
// Returns position of element aij of a the upper triangular 
// part of a Nel x Nel matrix 
Returns position of element aij of a the upper triangular
part of a Nel x Nel matrix
Returns position of element aij of a the upper triangular part of a Nel x Nel matrix
int ij2rUpTr(int Nel, int i, int j)
{
  if ((i>j)||(i<0)||(j<0)) return(0);
  else
    return(i*Nel-i*(i-1)/2+(j-i));
}
int ij2rUpTr ( int Nel , int i , int j ) { if ( ( i > j ) || ( i < 0 ) || ( j < 0 ) ) return ( 0 ) ; else return ( i * Nel - i * ( i - 1 ) / 2 + ( j - i ) ) ; }
Returns position of element aij of a the upper triangular
part of a Nel x Nel matrix
Returns position of element aij of a the upper triangular part of a Nel x Nel matrix