Add test and a sanity check

Author: pd3

plugins/fill-tags.c

@@ -318,6 +318,9 @@ int ftf_filter_expr(args_t *args, bcf1_t *rec, pop_t *pop, ftf_t *ftf)
     }
     else
     {
+        if ( nval1*rec->n_sample > nval )
+            error("Error: The expression returned %d value(s) total but %d value(s) per sample are required; there are %d samples\n",nval,nval1,rec->n_sample);
+
         int nfill = ftf->len==BCF_VL_FIXED ? ftf->cnt : nval1;
         int ncopy = nval1 < nfill ? nval1 : nfill;
         if ( ftf->type==BCF_HT_REAL )

test/test.pl

@@ -682,6 +682,7 @@
 run_test(\&test_vcf_plugin,$opts,in=>'query.variantkey',out=>'query.add-variantkey.vcf',cmd=>'+add-variantkey',args=>'');
 run_test(\&test_vcf_plugin,$opts,in=>'query.variantkey',out=>'variantkey-hex.out',cmd=>'+variantkey-hex',args=>'test/');
 run_test(\&test_vcf_plugin,$opts,in=>'query.nucleotide',out=>'query.allele-length.tsv',cmd=>'+allele-length',args=>'');
+run_test(\&test_vcf_plugin,$opts,in=>'query.func.1',out=>'fill-tags.func.1.out',cmd=>'+fill-tags --no-version',args=>q[-- -t 'FMT/AB:1=float(FMT/AD[*:0] / ssum(FMT/AD[*]))' | grep -v ^#]);
 run_test(\&test_vcf_plugin,$opts,in=>'fmissing',out=>'fmissing.1.out',cmd=>'+fill-tags --no-version',args=>q[-- -S {PATH}/fmissing.txt -t 'F_MISSING' | grep -v ^#]);
 run_test(\&test_vcf_plugin,$opts,in=>'fmissing',out=>'fmissing.1.out',cmd=>'+fill-tags --no-version',args=>q[-- -S {PATH}/fmissing.txt -t 'F_MISSING:1=F_PASS(GT="mis")' | grep -v ^#]);
 run_test(\&test_vcf_plugin,$opts,in=>'fmissing',out=>'fmissing.2.out',cmd=>'+fill-tags --no-version',args=>q[-- -S {PATH}/fmissing.txt -t 'N_MISSING:1=int(N_PASS(GT="mis"))' | grep -v ^#]);