#======================================================================

#                    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;