Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

BBRKC 2017 Run error #195

Open
adletaw opened this issue Oct 25, 2017 · 5 comments
Open

BBRKC 2017 Run error #195

adletaw opened this issue Oct 25, 2017 · 5 comments
Assignees
Labels
Milestone

Comments

@adletaw
Copy link
Contributor

adletaw commented Oct 25, 2017

+--------------------------+
| Beginning report section |
+--------------------------+
Index bounds do not match in dvector elem_prod(const dvector&, const dvector&)
Makefile:25: recipe for target 'run' failed
make: *** [run] Error 1

@adletaw adletaw added the bug label Oct 25, 2017
@adletaw adletaw self-assigned this Oct 25, 2017
@adletaw
Copy link
Contributor Author

adletaw commented Oct 25, 2017

Error occurs when model reaches line 3947 of src/gmacs.tpl

get_all_sdnr_MAR();

Tracks back to line 3380 from within the get_all_sdnr_MAR() function:

Francis_weights = calc_Francis_weights();

In the calc_Francis_weights() function, the model fails when k = 2 and i = 1 at line 3439:

Obs = sum(elem_prod(d3_obs_size_comps(k,i), mid_points));

Error occurs again here when kk = 2 and ii = 1:

// Compute effective N's
for ( int kk = 1; kk <= nSizeComps; kk++ )
{
	for ( int ii = 1; ii <= nSizeCompRows(kk); ii++ )
	{
		double sdl_tmp     = Sd_length(d3_obs_size_comps(kk,ii));
		effN(kk,ii)        = Eff_N(d3_obs_size_comps(kk,ii), d3_pre_size_comps(kk,ii));
		effN2(kk,ii)       = Eff_N2(d3_obs_size_comps(kk,ii), d3_pre_size_comps(kk,ii));
		pre_mn_size(kk,ii) = mn_length(d3_pre_size_comps(kk,ii));
		obs_mn_size(kk,ii) = mn_length(d3_obs_size_comps(kk,ii));
		lb_mn_size(kk,ii)  = obs_mn_size(kk,ii) - sdl_tmp * 2.0 / sqrt(size_comp_sample_size(kk,ii));
		ub_mn_size(kk,ii)  = obs_mn_size(kk,ii) + sdl_tmp * 2.0 / sqrt(size_comp_sample_size(kk,ii));
		/*
		*/
	}
}

@quantifish
Copy link
Contributor

Ah I think I had this commented out - might be good to start by commenting out and I can take a look at it.

@adletaw
Copy link
Contributor Author

adletaw commented Oct 25, 2017

It looks like a similar problem came up last year in #167

@adletaw adletaw added this to the BBRKC 2017 milestone Oct 25, 2017
@quantifish
Copy link
Contributor

quantifish commented Oct 25, 2017 via email

@adletaw
Copy link
Contributor Author

adletaw commented Oct 25, 2017

Ah okay. Good to know.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants