#======================================================================
#                    . M A T C H _ M I N I M I Z E . M A D 
#                    doc: Tue Aug 22 18:31:22 2006
#                    dlm: Wed Aug 23 23:15:13 2006
#                    (c) 2006 A.M. Thurnherr
#                    uE-Info: 39 26 NIL 0 0 72 2 2 4 NIL ofnI
#======================================================================

# sum of absolute deviations

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

sub MM_usage(@)
{
	@MM_fname = @_;

	croak("match-minimize usage: -m SAD(<field>[,...])\n")
		unless (@MM_fname > 0);
	for (my($fi)=0; $fi<@MM_fname; $fi++) {
		$MM_fnr[$fi]  =  &fnr($MM_fname[$fi]);
		$MM_wfnr[$fi] = &wfnr($MM_fname[$fi]);
		&IS_init($MM_wfnr[$fi],$MM_cwfnr);			# prepare interpolation
	}
	for (my($r)=0; $r<@ants_; $r++) {
		for (my($fi)=0; $fi<@MM_fname; $fi++) {
			croak("$0: $MM_fname[$fi] must be numeric\n")
				unless numberp($ants_[$r][$MM_fnr[$fi]]) &&
					   numberp($ants_[$r][$MM_wfnr[$fi]]);
		}
	}
}

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

sub MM_eval(@)
{
	&MW_warp(@_);
	
	my($sad) = my($nad) = 0;
	for (my($r)=0; $r<@ants_; $r++) {
		my($xi) = &bSearch($MW_cwfnr,$ants_[$r][$MW_cfnr]);
		next unless defined($xi);
		$nad++;
		for (my($fi)=0; $fi<@MM_fnr; $fi++) {
			my($wv) = &IS_interpolate(\@wf_,$MW_cwfnr,$ants_[$r][$MW_cfnr],$xi,$MM_wfnr[$fi]);
			$sad += abs($ants_[$r][$MM_fnr[$fi]] - $wv);
		}
	}
	return $sad/$nad;
}

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

1;