Branch Coverage

pp-hmmexpect.c

Criterion	Covered	Total	%
branch	66	304	21.7

line	true	false	branch
41	39	45	if (y-x > LOG_BIG) return y;
42	23	22	else if (x-y > LOG_BIG) return x;
44	8	14	else if (x
81	0	6	if (!__privtrans->broadcast.incs) return PDL->make_error(PDL_EUSERERROR, "Error in hmmexpect:" "broadcast.incs NULL");
192	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Float,F)
193	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
202	0	0	{/* Open N */ PDL_EXPAND2(register PDL_Indx N=0, __N_stop=(__N_size)); for(; N<__N_stop; N+=1) { p_o = logadd(p_o, (omega_datap)[0+(__inc_omega_N*(N))] + (alpha_datap)[0+(__inc_alpha_N*(N))+(__inc_alpha_T*(t))]); }} /* Close N */
207	0	0	for (i=0; i<__privtrans->ind_sizes[1]; i++) {
214	0	0	for (t--; t>=0; t--) {
219	0	0	for (i=0; i<__privtrans->ind_sizes[1]; i++) {
223	0	0	for (j=0; j<__privtrans->ind_sizes[1]; j++) {
230	0	0	if (t==0) (epi_datap)[0+(__inc_epi_N*(i))] = logadd(gamma_it, (epi_datap)[0+(__inc_epi_N*(i))]);
236	0	0	}PDL_BROADCASTLOOP_END_hmmexpect_readdata
	0	0	}PDL_BROADCASTLOOP_END_hmmexpect_readdata
239	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	6	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	6	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	6	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	6	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	6	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	6	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	6	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	6	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	6	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	6	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	6	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
	0	6	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_Double,D)
240	0	6	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	6	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	6	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	6	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	6	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	6	6	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	6	6	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
249	12	6	{/* Open N */ PDL_EXPAND2(register PDL_Indx N=0, __N_stop=(__N_size)); for(; N<__N_stop; N+=1) { p_o = logadd(p_o, (omega_datap)[0+(__inc_omega_N*(N))] + (alpha_datap)[0+(__inc_alpha_N*(N))+(__inc_alpha_T*(t))]); }} /* Close N */
254	12	6	for (i=0; i<__privtrans->ind_sizes[1]; i++) {
261	6	6	for (t--; t>=0; t--) {
266	12	6	for (i=0; i<__privtrans->ind_sizes[1]; i++) {
270	24	12	for (j=0; j<__privtrans->ind_sizes[1]; j++) {
277	12	0	if (t==0) (epi_datap)[0+(__inc_epi_N*(i))] = logadd(gamma_it, (epi_datap)[0+(__inc_epi_N*(i))]);
283	0	6	}PDL_BROADCASTLOOP_END_hmmexpect_readdata
	0	6	}PDL_BROADCASTLOOP_END_hmmexpect_readdata
286	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
	0	0	PDL_DECLARE_PARAMS_hmmexpect_1(PDL_LDouble,E)
287	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
	0	0	PDL_BROADCASTLOOP_START_hmmexpect_readdata {
296	0	0	{/* Open N */ PDL_EXPAND2(register PDL_Indx N=0, __N_stop=(__N_size)); for(; N<__N_stop; N+=1) { p_o = logadd(p_o, (omega_datap)[0+(__inc_omega_N*(N))] + (alpha_datap)[0+(__inc_alpha_N*(N))+(__inc_alpha_T*(t))]); }} /* Close N */
301	0	0	for (i=0; i<__privtrans->ind_sizes[1]; i++) {
308	0	0	for (t--; t>=0; t--) {
313	0	0	for (i=0; i<__privtrans->ind_sizes[1]; i++) {
317	0	0	for (j=0; j<__privtrans->ind_sizes[1]; j++) {
324	0	0	if (t==0) (epi_datap)[0+(__inc_epi_N*(i))] = logadd(gamma_it, (epi_datap)[0+(__inc_epi_N*(i))]);
330	0	0	}PDL_BROADCASTLOOP_END_hmmexpect_readdata
	0	0	}PDL_BROADCASTLOOP_END_hmmexpect_readdata
372	0	6	if (!PDL) return (pdl_error){PDL_EFATAL, "PDL core struct is NULL, can't continue",0};
374	0	6	if (!__privtrans) return PDL->make_error_simple(PDL_EFATAL, "Couldn't create trans");
386	0	6	PDL_RETERROR(PDL_err, PDL->type_coerce(__privtrans));
387	0	6	PDL_RETERROR(PDL_err, PDL->make_trans_mutual(__privtrans));