#======================================================================
#                    C O V S R T . P L 
#                    doc: Sun Sep 26 18:44:11 1999
#                    dlm: Sun Sep 26 18:56:56 1999
#                    (c) 1999 A.M. Thurnherr
#                    uE-Info: 46 2 NIL 0 0 72 2 2 4 ofnI
#======================================================================

# 2nd edition covsrt.c adapted to ANTS

# HISTORY:
#	Sep 26, 1999: - created after confusion about old version [covsrt_old.pl]

sub covsrt($$)
{
	my($covarR,$iaR) = @_;
	my($ma) = $#{$covarR};
	my($mfit) = $#{$iaR};
	my($i,$j,$k);
	my($swap);

	for ($i=$mfit+1;$i<=$ma;$i++) {
		for ($j=1;$j<=$i;$j++) {
			$covarR->[$i][$j] = 0;
			$covarR->[$j][$i] = 0;
		}
	}
	$k=$mfit;
	for ($j=$ma;$j>=1;$j--) {
		if ($iaR->[$j]) {
			for ($i=1;$i<=$ma;$i++) {
				$swap = $covarR->[$i][$k];
				$covarR->[$i][$k] = $covarR->[$i][$j];
				$covarR->[$i][$j] = $swap;
			}
			for ($i=1;$i<=$ma;$i++) {
				$swap = $covarR->[$k][$i];
				$covarR->[$k][$i] = $covarR->[$j][$i];
				$covarR->[$j][$i] = $swap;
			}
			$k--;
		}
	}
}

1;