#======================================================================
#                    . M A T C H _ W A R P . S T R E T C H 
#                    doc: Tue Aug 22 18:31:22 2006
#                    dlm: Wed Aug 23 23:51:25 2006
#                    (c) 2006 A.M. Thurnherr
#                    uE-Info: 52 0 NIL 0 0 72 2 2 4 NIL ofnI
#======================================================================

# stretch (and shift!) monotonically increasing coordinate

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

sub MW_usage(@)
{
	croak("match-warp usage: -w stretch(<field>)\n")
		unless (@_ == 1);
	($MW_cfname) = @_;
	
	$MW_cfnr  =  fnr($MW_cfname);
	$MW_cwfnr = wfnr($MW_cfname);

	croak("$0: field $MW_cfname must be numeric\n")
		unless numberp($ants_[0][$MW_cfnr]);
	for (my($r)=1; $r<@ants_; $r++) {
		croak("$0: field $MW_cfname of rec $r must be numeric ($ants_[$r][$MW_cfnr])\n")
			unless numberp($ants_[$r][$MW_cfnr]);
		croak("$0: field $MW_cfname must be monotonically increasing (" .
				  "$ants_[$r][$MW_cfnr]=>$ants_[$r-1][$MW_cfnr]" .
			  	  ")\n")
			unless ($ants_[$r][$MW_cfnr] >= $ants_[$r-1][$MW_cfnr]);
	}

	croak("$0: warp-file field $MW_cfname must be numeric\n")
		unless numberp($wf_[0][$MW_cwfnr]);
	for (my($r)=1; $r<@wf_; $r++) {
		croak("$0: warp-file field $MW_cfname of rec $r must be numeric ($wf_[$r][$MW_cwfnr])\n")
			unless numberp($wf_[$r][$MW_cwfnr]);
		croak("$0: warp-file field $MW_cfname must be monotonically increasing (" .
				  "$wf_[$r][$MW_cwfnr]=>$wf_[$r-1][$MW_cwfnr]" .
			  	  ")\n")
			unless ($wf_[$r][$MW_cwfnr]  >= $wf_[$r-1][$MW_cwfnr]);
		$MW_unwarped[$r] = $wf_[$r][$MW_cwfnr];
	}
}

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

sub MW_warp(@)
{
	my($offset,$scale) = @_;

	for (my($r)=0; $r<@wf_; $r++) {
		$wf_[$r][$MW_cwfnr] = $MW_unwarped[$r]*$scale + $offset;
	}
}

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

1;