*  
*  NOTE:  This procedure requires RATS v.5.11 or higher.  
*  
*  This procedure generates simulated out-of-sample forecasts, forecast errors, 
*  and tests of equal mean square error (MSE) and encompassing for a pair of 
*  nested models (the first model is a restricted version of the second) of 
*  a scalar predictand.  Critical values are reported along with the test statistics.
*  For detail greater than that provided, a (PDF) documentation file can be obtained
*  from the following web pages (look for a link to the procedure and the documentation):
*    http://www.missouri.edu/~mwmd4f/
*    http://www.kansascityfed.org/econres/staff/tec.htm
*
*  Following the nomenclature of the references below, the tests are referred to as 
*  MSE-t, MSE-F, ENC_t, and ENC-F.  Note that, as explained in the references,
*  with nested models, the tests are one-sided (rejecting only if the test statistic
*  is larger than the alpha-th percentile, for a significance level of 1-alpha).
*
*  References:
*   Clark, Todd. E., and Michael W. McCracken, 2001, "Tests of Equal Forecast Accuracy 
*   and Encompassing for Nested Models," Journal of Econometrics 105 (Nov.), pp. 85-110.
*
*   McCracken, Michael W., 2004, "Asymptotics for Out of Sample Tests of Granger 
*   Causality," manuscript, University of Missouri.
*
*  Written by Taisuke Nakata and Todd Clark, Federal Reserve Bank of Kansas City
*  Created September 2004.
*
**** Syntax:
*
*  forecastproc(options) LHS fcst start end fc_usr error_usr
*  # variables_list(1) 
*  # variables_list(2)
*  
*   -Parameters:
*  LHS: Dependent Variable
*  fcst (required) : the first period for which forecast is generated.    
*  start (optional) : the beginning of the estimation sample.
*  end (optional) : the last period for which forecast is generated.
*  fc_usr (optional) : the series (a vector of 2 series) to store the forecast values
*  error_usr (optional) : the series (a vector of 2 series) to store the forecast errors
*                          (if you only want the errors and not the forecasts, use a *
*                           in the forecast series position)
*

*   -Options:
*     print/[noprint] :print forecasts (default is noprint)
*     scheme = 1 fixed  (in proc options, user can refer to schemes by either name or #)
*              2 rolling 
*              3 recursive (Default)
*
*  Note on use of the rolling scheme:  the number of observations used in the rolling
*   sample will be (fcst-1)-start+1.  This means that, for a given forecast starting
*   point, to use a particular window size the user must appropriately set the starting
*   point for estimation (start).
*   
*    
*   -Supplementary Card # (required): 
*  variables_list(1) : Variables to be included in the restricted model.
*  variables_list(2) : Additional variables to be included in the unrestricted model.
*    *Users have to include "constant" among the variables if they want it to be included in the regression. 
*
*   -Comments:  
*       If users DO NOT specify the sample start and end, then this procedure will find and use 
*  the period for which all series used in the procedure have non-missing data. 
*       Even when users DO specify the sample start and end, it still will find and use, within the period    
*  specified by users, the sub-period for which all series used in the procedure have non-missing data.
*
*  -Structure of the procedure (after preliminaries)
*  Step 1 generates forecast for restricted model
*  Step 2 generates forecast for unrestricted model
*  Step 3 computes summary stats and MSE-F, MSE-t, ENC-F, and ENC-t statistics. 
*  Step 4 pulls the appropriate critical values from tables read in at the start of the procedure
*  Step 5 displays the results.

proc forecastproc LHS fcst start end fc_usr error_usr

type series LHS
type integer start end fcst
type vector[series] *fc_usr *error_usr

option choice scheme 3 fixed rolling recursive
option switch print 0

local vec[integer] res_ser unres_ser modelcheck
local integer stser endser stshared endshared P R maxnsig maxnpi maxnk cmpos1 cmpos2 ndf
local vec[real] me mae rmse mse
local vec[series] fc error error2
local vec cvpis
local real pival pidist
local series assignser
local vec[rec] cv_msef cv_mset cv_encf cv_enct

***********************************
*********************************** reading in matrices of critical values
***********************************
fixed vec[rec] reccv_mset(3)(20,20)
0.885	0.780	0.657	0.598	0.512	0.443	0.402	0.370	0.330	0.306	0.281	0.172	0.111	0.055	0.000	-0.022	-0.047	-0.075	-0.100	-0.239
0.932	0.786	0.614	0.541	0.455	0.361	0.295	0.253	0.235	0.194	0.160	0.026	-0.064	-0.144	-0.218	-0.261	-0.294	-0.336	-0.374	-0.529
0.939	0.751	0.551	0.454	0.356	0.279	0.222	0.175	0.108	0.074	0.035	-0.131	-0.249	-0.320	-0.397	-0.448	-0.491	-0.537	-0.583	-0.772
0.898	0.742	0.562	0.419	0.263	0.169	0.094	0.052	-0.014	-0.054	-0.106	-0.278	-0.413	-0.483	-0.557	-0.624	-0.676	-0.734	-0.779	-1.010
0.866	0.694	0.461	0.315	0.179	0.062	-0.021	-0.083	-0.145	-0.174	-0.228	-0.421	-0.562	-0.656	-0.736	-0.802	-0.853	-0.904	-0.950	-1.209
0.823	0.642	0.394	0.256	0.108	-0.011	-0.101	-0.164	-0.218	-0.266	-0.319	-0.518	-0.696	-0.793	-0.864	-0.945	-1.002	-1.053	-1.094	-1.369
0.811	0.615	0.359	0.213	0.062	-0.088	-0.152	-0.230	-0.305	-0.363	-0.449	-0.629	-0.796	-0.891	-0.993	-1.069	-1.123	-1.181	-1.225	-1.533
0.773	0.574	0.329	0.139	0.003	-0.131	-0.203	-0.293	-0.383	-0.452	-0.516	-0.725	-0.889	-1.014	-1.110	-1.181	-1.237	-1.289	-1.344	-1.667
0.733	0.561	0.273	0.096	-0.068	-0.187	-0.286	-0.377	-0.437	-0.518	-0.579	-0.812	-0.987	-1.102	-1.207	-1.277	-1.348	-1.391	-1.448	-1.813
0.749	0.529	0.226	0.032	-0.130	-0.248	-0.355	-0.454	-0.524	-0.591	-0.651	-0.928	-1.100	-1.233	-1.333	-1.413	-1.496	-1.545	-1.611	-1.970
0.734	0.503	0.195	-0.016	-0.187	-0.323	-0.418	-0.507	-0.602	-0.657	-0.738	-0.994	-1.192	-1.334	-1.439	-1.522	-1.589	-1.656	-1.717	-2.099
0.706	0.451	0.155	-0.061	-0.242	-0.375	-0.469	-0.575	-0.660	-0.726	-0.798	-1.103	-1.294	-1.431	-1.551	-1.622	-1.706	-1.779	-1.852	-2.238
0.694	0.451	0.121	-0.118	-0.291	-0.421	-0.536	-0.631	-0.718	-0.801	-0.870	-1.158	-1.370	-1.512	-1.616	-1.725	-1.796	-1.872	-1.951	-2.333
0.674	0.410	0.090	-0.152	-0.317	-0.480	-0.599	-0.718	-0.810	-0.873	-0.939	-1.249	-1.475	-1.614	-1.722	-1.807	-1.897	-1.974	-2.046	-2.443
0.653	0.379	0.048	-0.191	-0.384	-0.532	-0.669	-0.774	-0.857	-0.936	-1.001	-1.318	-1.557	-1.703	-1.828	-1.899	-2.001	-2.077	-2.153	-2.557
0.636	0.376	0.027	-0.228	-0.427	-0.564	-0.692	-0.825	-0.933	-0.997	-1.063	-1.379	-1.614	-1.762	-1.891	-1.989	-2.076	-2.153	-2.229	-2.668
0.597	0.352	-0.003	-0.268	-0.455	-0.631	-0.746	-0.887	-0.978	-1.074	-1.131	-1.436	-1.687	-1.851	-1.966	-2.065	-2.164	-2.230	-2.323	-2.760
0.574	0.293	-0.059	-0.315	-0.517	-0.692	-0.802	-0.924	-1.038	-1.129	-1.188	-1.509	-1.748	-1.924	-2.055	-2.143	-2.264	-2.328	-2.410	-2.863
0.554	0.267	-0.079	-0.373	-0.563	-0.739	-0.842	-0.974	-1.091	-1.187	-1.253	-1.584	-1.835	-1.997	-2.141	-2.234	-2.342	-2.418	-2.498	-2.966
0.539	0.262	-0.136	-0.390	-0.610	-0.776	-0.893	-1.026	-1.148	-1.239	-1.305	-1.654	-1.900	-2.078	-2.204	-2.326	-2.438	-2.512	-2.595	-3.064
																			
1.245	1.111	0.994	0.971	0.863	0.771	0.740	0.705	0.671	0.638	0.610	0.509	0.425	0.378	0.309	0.290	0.255	0.215	0.192	0.045
1.274	1.140	0.986	0.868	0.782	0.704	0.623	0.596	0.537	0.507	0.478	0.340	0.235	0.152	0.094	0.054	0.009	-0.027	-0.077	-0.265
1.300	1.120	0.968	0.808	0.685	0.610	0.552	0.496	0.438	0.419	0.386	0.225	0.066	-0.023	-0.065	-0.140	-0.191	-0.248	-0.271	-0.519
1.264	1.101	0.914	0.772	0.609	0.502	0.419	0.345	0.285	0.239	0.221	0.026	-0.110	-0.187	-0.268	-0.340	-0.393	-0.438	-0.496	-0.746
1.222	1.061	0.849	0.689	0.491	0.386	0.308	0.224	0.148	0.107	0.081	-0.104	-0.277	-0.372	-0.470	-0.534	-0.579	-0.614	-0.672	-0.952
1.192	0.998	0.768	0.615	0.429	0.328	0.259	0.141	0.078	0.055	-0.019	-0.238	-0.394	-0.488	-0.561	-0.637	-0.706	-0.756	-0.816	-1.122
1.199	0.976	0.742	0.546	0.372	0.279	0.191	0.072	-0.002	-0.034	-0.105	-0.313	-0.498	-0.584	-0.677	-0.758	-0.811	-0.876	-0.958	-1.271
1.193	0.928	0.677	0.462	0.302	0.198	0.105	0.020	-0.058	-0.101	-0.176	-0.413	-0.593	-0.687	-0.793	-0.877	-0.944	-1.027	-1.078	-1.398
1.112	0.912	0.617	0.397	0.276	0.121	0.030	-0.055	-0.122	-0.193	-0.257	-0.518	-0.691	-0.821	-0.906	-1.010	-1.087	-1.145	-1.223	-1.541
1.082	0.890	0.566	0.358	0.205	0.043	-0.072	-0.162	-0.222	-0.296	-0.339	-0.633	-0.821	-0.950	-1.065	-1.163	-1.228	-1.286	-1.369	-1.717
1.061	0.862	0.544	0.278	0.137	0.000	-0.109	-0.229	-0.300	-0.350	-0.408	-0.710	-0.899	-1.043	-1.164	-1.257	-1.336	-1.401	-1.472	-1.847
1.037	0.822	0.498	0.246	0.062	-0.058	-0.169	-0.288	-0.366	-0.423	-0.502	-0.784	-1.000	-1.165	-1.273	-1.358	-1.449	-1.511	-1.587	-1.967
1.017	0.771	0.470	0.214	0.013	-0.093	-0.223	-0.333	-0.427	-0.477	-0.565	-0.863	-1.078	-1.240	-1.354	-1.435	-1.517	-1.594	-1.657	-2.061
0.996	0.731	0.432	0.168	-0.028	-0.157	-0.273	-0.390	-0.485	-0.561	-0.641	-0.941	-1.149	-1.314	-1.456	-1.530	-1.653	-1.704	-1.782	-2.202
1.003	0.737	0.388	0.139	-0.082	-0.204	-0.335	-0.448	-0.560	-0.647	-0.711	-1.019	-1.240	-1.398	-1.533	-1.623	-1.740	-1.789	-1.870	-2.309
0.984	0.714	0.342	0.096	-0.103	-0.253	-0.381	-0.496	-0.600	-0.706	-0.762	-1.099	-1.305	-1.497	-1.597	-1.706	-1.816	-1.889	-1.954	-2.439
0.951	0.715	0.309	0.056	-0.159	-0.318	-0.442	-0.544	-0.657	-0.733	-0.829	-1.160	-1.384	-1.560	-1.679	-1.789	-1.896	-1.978	-2.032	-2.532
0.927	0.653	0.266	-0.010	-0.208	-0.372	-0.483	-0.619	-0.734	-0.815	-0.917	-1.214	-1.457	-1.633	-1.756	-1.857	-1.975	-2.048	-2.137	-2.647
0.917	0.631	0.227	-0.068	-0.258	-0.412	-0.557	-0.670	-0.812	-0.883	-0.962	-1.269	-1.546	-1.716	-1.846	-1.945	-2.071	-2.141	-2.223	-2.710
0.894	0.602	0.237	-0.082	-0.281	-0.470	-0.589	-0.710	-0.834	-0.938	-1.015	-1.335	-1.606	-1.780	-1.925	-2.027	-2.157	-2.221	-2.322	-2.823
																			
1.921	1.784	1.625	1.515	1.462	1.436	1.413	1.343	1.316	1.274	1.238	1.133	1.059	0.982	0.952	0.912	0.848	0.834	0.789	0.650
1.986	1.856	1.563	1.436	1.387	1.312	1.276	1.196	1.158	1.127	1.074	0.903	0.815	0.702	0.631	0.580	0.535	0.490	0.460	0.273
1.840	1.737	1.542	1.448	1.359	1.252	1.148	1.071	0.976	0.978	0.953	0.797	0.648	0.465	0.418	0.335	0.276	0.245	0.233	-0.013
1.872	1.731	1.581	1.365	1.195	1.119	1.108	1.041	0.902	0.861	0.854	0.614	0.462	0.330	0.269	0.190	0.141	0.107	0.084	-0.235
1.849	1.679	1.468	1.242	1.095	0.995	0.979	0.913	0.795	0.732	0.677	0.436	0.301	0.178	0.102	0.052	0.010	-0.058	-0.115	-0.436
1.836	1.639	1.390	1.200	1.042	0.943	0.859	0.755	0.686	0.610	0.593	0.350	0.181	0.064	-0.017	-0.092	-0.151	-0.228	-0.228	-0.645
1.836	1.649	1.341	1.154	0.994	0.872	0.810	0.637	0.549	0.476	0.438	0.225	0.039	-0.097	-0.173	-0.227	-0.304	-0.342	-0.422	-0.791
1.789	1.659	1.298	1.090	0.879	0.788	0.728	0.503	0.444	0.401	0.359	0.153	-0.029	-0.207	-0.293	-0.343	-0.434	-0.434	-0.535	-0.928
1.813	1.607	1.268	1.112	0.804	0.724	0.634	0.523	0.427	0.391	0.305	0.022	-0.085	-0.305	-0.385	-0.470	-0.552	-0.588	-0.657	-1.083
1.743	1.534	1.193	1.035	0.758	0.621	0.506	0.419	0.347	0.285	0.185	-0.046	-0.247	-0.409	-0.469	-0.568	-0.664	-0.708	-0.766	-1.188
1.736	1.464	1.135	0.917	0.694	0.535	0.439	0.287	0.262	0.193	0.124	-0.103	-0.316	-0.456	-0.555	-0.673	-0.732	-0.808	-0.888	-1.294
1.694	1.453	1.056	0.868	0.617	0.452	0.404	0.262	0.142	0.135	0.048	-0.199	-0.432	-0.588	-0.662	-0.783	-0.889	-0.941	-1.023	-1.464
1.676	1.458	0.999	0.818	0.585	0.372	0.332	0.197	0.093	0.062	0.007	-0.315	-0.538	-0.657	-0.790	-0.887	-0.992	-1.032	-1.111	-1.599
1.611	1.383	1.019	0.759	0.538	0.371	0.256	0.148	0.038	-0.032	-0.094	-0.372	-0.614	-0.811	-0.903	-1.018	-1.091	-1.155	-1.197	-1.718
1.631	1.352	0.967	0.690	0.463	0.283	0.158	0.024	-0.012	-0.089	-0.145	-0.485	-0.711	-0.862	-0.997	-1.095	-1.214	-1.275	-1.334	-1.816
1.657	1.353	0.922	0.672	0.428	0.251	0.105	-0.007	-0.063	-0.154	-0.195	-0.535	-0.801	-0.942	-1.054	-1.203	-1.274	-1.359	-1.472	-1.917
1.587	1.316	0.895	0.626	0.411	0.238	0.118	-0.023	-0.094	-0.230	-0.283	-0.601	-0.885	-1.042	-1.153	-1.282	-1.381	-1.441	-1.550	-2.021
1.630	1.288	0.916	0.648	0.379	0.181	0.030	-0.076	-0.170	-0.241	-0.335	-0.717	-0.957	-1.110	-1.230	-1.361	-1.494	-1.524	-1.660	-2.167
1.605	1.261	0.867	0.603	0.375	0.133	-0.007	-0.170	-0.238	-0.296	-0.387	-0.809	-1.024	-1.192	-1.331	-1.459	-1.587	-1.635	-1.758	-2.292
1.578	1.282	0.826	0.540	0.318	0.090	-0.055	-0.200	-0.273	-0.346	-0.444	-0.805	-1.048	-1.259	-1.416	-1.561	-1.648	-1.753	-1.836	-2.319

fixed vec[rec] reccv_msef(3)(20,20)
0.530	0.659	0.814	0.796	0.798	0.751	0.759	0.698	0.685	0.687	0.616	0.422	0.321	0.163	0.000	-0.072	-0.155	-0.243	-0.362	-0.916
0.794	0.912	1.029	1.077	1.008	0.880	0.785	0.697	0.666	0.587	0.506	0.088	-0.238	-0.597	-0.926	-1.150	-1.333	-1.560	-1.724	-2.948
0.970	1.064	1.117	1.121	0.960	0.857	0.691	0.599	0.386	0.276	0.127	-0.583	-1.156	-1.576	-2.019	-2.356	-2.656	-2.968	-3.397	-5.116
1.060	1.225	1.313	1.184	0.829	0.545	0.354	0.197	-0.058	-0.234	-0.456	-1.329	-2.234	-2.741	-3.300	-3.766	-4.239	-4.678	-5.104	-7.549
1.130	1.277	1.228	0.958	0.614	0.241	-0.099	-0.361	-0.656	-0.820	-1.072	-2.254	-3.425	-4.102	-4.822	-5.431	-5.963	-6.499	-6.851	-10.055
1.220	1.313	1.164	0.890	0.419	-0.044	-0.405	-0.776	-1.072	-1.395	-1.664	-3.152	-4.451	-5.309	-6.274	-7.035	-7.608	-8.220	-8.722	-12.393
1.272	1.333	1.118	0.799	0.242	-0.363	-0.728	-1.194	-1.657	-2.033	-2.507	-4.090	-5.612	-6.620	-7.649	-8.595	-9.202	-9.893	-10.552	-14.844
1.338	1.369	1.058	0.552	0.014	-0.632	-1.076	-1.633	-2.174	-2.731	-3.160	-4.995	-6.734	-8.005	-9.133	-10.167	-10.916	-11.643	-12.444	-17.451
1.354	1.432	0.920	0.393	-0.327	-1.007	-1.595	-2.229	-2.666	-3.250	-3.794	-6.002	-7.716	-9.235	-10.643	-11.635	-12.677	-13.506	-14.260	-19.907
1.458	1.401	0.884	0.155	-0.600	-1.341	-2.008	-2.782	-3.348	-3.839	-4.437	-7.100	-9.149	-10.804	-12.223	-13.402	-14.492	-15.364	-16.557	-22.716
1.450	1.404	0.756	-0.073	-0.968	-1.782	-2.517	-3.219	-3.975	-4.597	-5.310	-8.089	-10.317	-12.240	-13.780	-15.137	-16.293	-17.452	-18.462	-25.490
1.492	1.358	0.601	-0.309	-1.371	-2.179	-2.980	-3.845	-4.636	-5.233	-5.980	-9.211	-11.758	-13.827	-15.550	-16.996	-18.294	-19.611	-20.656	-28.295
1.524	1.360	0.499	-0.583	-1.606	-2.628	-3.506	-4.418	-5.195	-6.055	-6.758	-10.276	-12.963	-15.226	-16.934	-18.576	-19.989	-21.417	-22.675	-30.910
1.556	1.286	0.380	-0.767	-1.882	-3.049	-4.004	-5.167	-6.144	-6.805	-7.745	-11.357	-14.257	-16.762	-18.551	-20.328	-21.831	-23.441	-24.710	-33.567
1.556	1.277	0.212	-1.078	-2.360	-3.436	-4.573	-5.773	-6.730	-7.600	-8.623	-12.392	-12.721	-15.501	-17.549	-19.094	-20.901	-22.092	-23.899	-33.303
1.552	1.296	0.115	-1.224	-2.604	-3.899	-4.965	-6.300	-7.478	-8.352	-9.239	-13.379	-16.770	-19.620	-21.835	-23.673	-25.385	-27.114	-28.710	-38.943
1.518	1.236	-0.020	-1.536	-2.944	-4.385	-5.549	-7.020	-8.175	-9.048	-10.102	-14.418	-18.177	-21.053	-23.632	-25.460	-27.386	-29.019	-30.477	-41.532
1.508	1.107	-0.295	-1.849	-3.436	-4.932	-6.222	-7.642	-8.833	-9.897	-10.922	-15.575	-19.430	-22.610	-25.146	-27.239	-29.376	-31.081	-32.617	-44.179
1.504	0.983	-0.421	-2.225	-3.788	-5.499	-6.752	-8.190	-9.618	-10.646	-11.815	-16.779	-20.907	-24.090	-26.786	-29.047	-31.164	-33.241	-34.948	-46.747
1.480	1.047	-0.722	-2.432	-4.203	-5.831	-7.332	-8.780	-10.249	-11.500	-12.465	-17.852	-22.209	-25.603	-28.634	-30.938	-33.312	-35.323	-37.109	-49.731
																			
0.850	1.038	1.298	1.554	1.567	1.548	1.583	1.623	1.599	1.553	1.518	1.389	1.357	1.292	1.096	0.994	0.873	0.829	0.680	0.236
1.184	1.453	1.733	1.891	1.820	1.802	1.819	1.752	1.734	1.692	1.706	1.259	0.948	0.748	0.415	0.265	0.053	-0.140	-0.393	-1.504
1.434	1.710	2.062	2.073	1.978	1.909	1.930	1.795	1.715	1.710	1.612	0.977	0.336	-0.116	-0.405	-0.789	-1.030	-1.421	-1.662	-3.511
1.566	1.964	2.246	2.194	1.900	1.809	1.578	1.376	1.256	1.122	1.029	0.150	-0.582	-1.103	-1.673	-2.183	-2.504	-2.923	-3.378	-5.803
1.688	2.082	2.235	2.242	1.773	1.449	1.316	1.045	0.718	0.502	0.459	-0.614	-1.687	-2.369	-3.152	-3.712	-4.121	-4.569	-5.068	-8.088
1.828	2.124	2.217	2.121	1.660	1.360	1.181	0.761	0.413	0.299	-0.109	-1.455	-2.570	-3.396	-4.217	-4.976	-5.657	-6.233	-6.782	-10.512
2.000	2.239	2.424	2.057	1.604	1.282	0.928	0.378	-0.008	-0.199	-0.591	-2.158	-3.571	-4.553	-5.390	-6.370	-7.082	-7.630	-8.442	-12.790
2.136	2.312	2.373	1.895	1.390	0.943	0.587	0.131	-0.372	-0.680	-1.140	-3.001	-4.520	-5.707	-6.805	-7.749	-8.583	-9.455	-10.257	-15.105
2.168	2.440	2.219	1.714	1.286	0.631	0.198	-0.356	-0.851	-1.241	-1.696	-3.886	-5.592	-6.913	-8.236	-9.268	-10.280	-11.245	-11.945	-17.917
2.202	2.489	2.157	1.536	1.055	0.205	-0.431	-1.071	-1.459	-1.988	-2.378	-4.924	-6.911	-8.449	-9.944	-11.147	-12.287	-13.182	-14.101	-20.469
2.218	2.487	2.204	1.320	0.725	0.002	-0.689	-1.507	-2.042	-2.529	-3.013	-5.824	-8.014	-9.702	-11.387	-12.804	-14.014	-15.019	-16.097	-23.018
2.244	2.474	2.063	1.225	0.352	-0.336	-1.060	-1.957	-2.615	-3.178	-3.855	-6.794	-9.289	-11.295	-13.091	-14.502	-15.697	-17.009	-17.984	-25.723
2.312	2.437	1.993	1.164	0.061	-0.623	-1.507	-2.460	-3.233	-3.689	-4.431	-7.777	-10.259	-12.589	-14.310	-15.843	-17.214	-18.552	-19.679	-28.242
2.398	2.441	1.897	0.885	-0.161	-1.001	-1.939	-2.914	-3.743	-4.602	-5.435	-8.914	-11.559	-13.972	-15.957	-17.413	-19.201	-20.552	-21.706	-30.991
2.454	2.453	1.814	0.712	-0.540	-1.385	-2.322	-3.538	-4.420	-5.329	-6.208	-9.927	-12.721	-15.501	-17.549	-19.094	-20.901	-22.092	-23.899	-33.303
2.434	2.486	1.670	0.570	-0.599	-1.799	-2.884	-3.857	-5.064	-5.926	-6.771	-10.843	-13.942	-16.786	-18.948	-20.724	-22.572	-24.265	-25.859	-35.805
2.436	2.541	1.539	0.304	-1.073	-2.197	-3.371	-4.577	-5.666	-6.463	-7.489	-11.874	-15.387	-18.285	-20.254	-22.281	-24.399	-26.016	-27.785	-38.537
2.428	2.432	1.373	-0.061	-1.362	-2.812	-3.736	-5.142	-6.345	-7.461	-8.446	-12.987	-16.609	-19.725	-22.088	-24.233	-26.307	-28.027	-29.868	-41.109
2.430	2.377	1.229	-0.449	-1.858	-3.163	-4.411	-5.761	-6.991	-8.135	-9.132	-13.903	-18.024	-21.072	-23.709	-25.705	-28.298	-30.011	-31.849	-43.657
2.430	2.346	1.240	-0.575	-1.936	-3.695	-4.864	-6.328	-7.553	-8.911	-10.113	-14.801	-19.277	-22.384	-24.859	-27.415	-30.155	-31.959	-33.754	-46.250
																			
1.608	2.129	2.768	3.179	3.459	3.584	3.771	3.589	3.838	3.882	3.951	4.056	4.251	3.813	3.900	3.665	3.724	3.741	3.659	3.187
1.996	2.691	3.426	3.907	4.129	4.200	4.362	4.304	4.309	4.278	4.250	3.991	3.799	3.582	3.300	3.032	2.930	2.744	2.759	1.737
2.418	3.092	4.080	4.136	4.322	4.341	4.337	4.192	4.089	4.365	4.184	3.869	3.277	2.685	2.561	2.159	1.826	1.637	1.690	-0.091
2.714	3.440	4.541	4.609	4.378	4.202	4.586	4.477	4.337	4.247	4.096	3.428	2.681	2.282	1.728	1.187	0.946	0.725	0.651	-2.051
2.902	3.673	4.466	4.434	4.249	4.351	4.349	4.187	3.945	3.783	3.783	2.618	1.932	1.307	0.769	0.427	0.071	-0.386	-1.009	-4.068
3.212	3.846	4.545	4.676	4.637	4.703	4.286	4.144	3.981	3.525	3.321	2.392	1.394	0.508	-0.152	-0.790	-1.349	-1.879	-2.171	-6.055
3.450	4.098	4.508	4.419	4.271	4.312	4.150	3.677	3.155	3.090	2.880	1.521	0.322	-0.845	-1.432	-2.133	-2.898	-3.262	-4.112	-8.336
3.408	4.130	4.645	4.625	4.202	4.147	3.912	3.185	2.933	2.952	2.484	1.111	-0.186	-1.852	-2.664	-3.252	-4.378	-4.374	-5.321	-10.384
3.540	4.388	4.703	4.873	4.122	4.066	3.753	3.027	2.925	2.802	2.186	0.150	-0.804	-2.994	-3.789	-4.587	-5.420	-6.044	-6.894	-12.523
3.646	4.433	4.813	4.718	3.944	3.645	3.194	2.578	2.282	2.152	1.436	-0.361	-2.249	-4.056	-4.865	-5.836	-7.034	-8.124	-8.731	-14.987
3.794	4.466	4.845	4.438	3.746	3.199	2.755	2.220	1.657	1.529	0.919	-0.898	-3.235	-4.535	-5.923	-7.388	-8.105	-9.442	-10.245	-17.161
4.052	4.619	4.686	4.157	3.548	2.791	2.533	1.868	1.179	1.065	0.417	-1.781	-4.528	-6.167	-7.378	-9.072	-10.076	-11.105	-12.255	-19.858
4.070	4.906	4.602	4.276	3.676	2.487	2.401	1.413	0.740	0.426	0.061	-2.834	-5.360	-7.240	-9.058	-10.300	-11.906	-12.721	-13.965	-22.478
3.972	4.891	4.701	4.312	3.479	2.420	1.852	1.007	0.347	-0.266	-0.776	-3.592	-6.422	-8.713	-10.538	-12.091	-13.795	-15.033	-15.777	-25.160
4.152	4.775	4.630	3.951	3.020	1.953	1.254	0.196	-0.098	-0.836	-1.388	-4.707	-7.570	-9.996	-12.097	-13.649	-15.336	-16.813	-18.139	-27.748
4.378	4.953	4.620	4.018	2.785	1.927	0.797	-0.056	-0.564	-1.373	-1.742	-5.658	-8.689	-11.054	-13.427	-15.016	-16.885	-18.588	-19.889	-30.393
4.362	4.924	4.591	3.951	2.902	1.670	0.897	-0.176	-0.829	-2.128	-2.738	-6.586	-10.765	-12.653	-14.245	-16.441	-18.821	-20.381	-22.052	-32.639
4.496	5.005	4.769	4.030	2.667	1.364	0.208	-0.650	-1.391	-2.328	-3.485	-7.754	-11.808	-13.999	-15.997	-18.578	-20.605	-22.316	-23.965	-35.546
4.396	4.941	4.648	3.990	2.553	1.010	-0.054	-1.475	-2.145	-2.913	-3.824	-8.929	-12.577	-15.450	-17.441	-20.141	-22.399	-24.204	-25.865	-38.295
4.542	5.117	4.750	3.579	2.279	0.730	-0.466	-1.980	-2.645	-3.484	-4.751	-9.391	-13.152	-16.179	-18.938	-21.657	-23.712	-25.916	-28.024	-40.968

fixed vec[rec] reccv_enct(3)(20,20)
1.056	1.002	1.005	0.995	0.968	0.955	0.942	0.939	0.959	0.956	0.939	0.937	0.928	0.922	0.923	0.922	0.935	0.931	0.912	0.914
1.166	1.101	1.086	1.096	1.077	1.066	1.060	1.044	1.049	1.053	1.035	1.034	1.031	1.028	1.015	1.012	1.020	0.997	1.015	1.021
1.227	1.138	1.105	1.136	1.118	1.113	1.124	1.111	1.098	1.094	1.114	1.083	1.052	1.074	1.060	1.053	1.060	1.061	1.064	1.074
1.219	1.175	1.192	1.177	1.136	1.132	1.111	1.099	1.104	1.112	1.111	1.091	1.084	1.090	1.089	1.078	1.076	1.078	1.057	1.090
1.205	1.192	1.170	1.172	1.117	1.092	1.116	1.092	1.084	1.102	1.100	1.080	1.073	1.074	1.088	1.057	1.069	1.079	1.078	1.081
1.222	1.187	1.162	1.169	1.125	1.109	1.125	1.140	1.132	1.126	1.137	1.104	1.081	1.122	1.100	1.107	1.100	1.110	1.118	1.119
1.223	1.178	1.180	1.195	1.164	1.142	1.168	1.169	1.165	1.137	1.134	1.136	1.140	1.127	1.134	1.124	1.134	1.132	1.120	1.124
1.211	1.195	1.205	1.193	1.181	1.159	1.187	1.169	1.169	1.152	1.151	1.169	1.133	1.144	1.118	1.139	1.141	1.155	1.129	1.143
1.213	1.223	1.203	1.208	1.187	1.166	1.181	1.190	1.179	1.172	1.168	1.184	1.153	1.157	1.159	1.172	1.201	1.177	1.149	1.182
1.244	1.226	1.208	1.194	1.176	1.175	1.177	1.178	1.162	1.175	1.170	1.152	1.130	1.142	1.129	1.151	1.190	1.162	1.139	1.169
1.250	1.236	1.224	1.213	1.199	1.169	1.165	1.192	1.182	1.185	1.185	1.162	1.157	1.175	1.143	1.160	1.167	1.171	1.167	1.175
1.245	1.232	1.230	1.199	1.202	1.193	1.183	1.180	1.186	1.180	1.183	1.148	1.161	1.167	1.159	1.170	1.180	1.182	1.160	1.176
1.267	1.238	1.236	1.209	1.215	1.193	1.198	1.225	1.215	1.203	1.210	1.169	1.156	1.169	1.176	1.191	1.193	1.186	1.186	1.192
1.270	1.230	1.232	1.199	1.214	1.187	1.193	1.193	1.204	1.203	1.193	1.174	1.165	1.170	1.180	1.193	1.186	1.189	1.181	1.180
1.261	1.241	1.234	1.194	1.194	1.166	1.191	1.184	1.187	1.184	1.198	1.168	1.156	1.167	1.172	1.176	1.172	1.186	1.177	1.200
1.263	1.260	1.244	1.212	1.232	1.213	1.213	1.198	1.183	1.206	1.202	1.183	1.168	1.163	1.186	1.195	1.176	1.182	1.200	1.172
1.247	1.262	1.245	1.221	1.217	1.225	1.247	1.219	1.212	1.199	1.208	1.201	1.163	1.157	1.195	1.202	1.194	1.186	1.208	1.184
1.257	1.240	1.237	1.221	1.221	1.199	1.241	1.222	1.194	1.190	1.206	1.204	1.163	1.185	1.186	1.215	1.208	1.211	1.196	1.218
1.257	1.239	1.231	1.212	1.232	1.196	1.222	1.230	1.224	1.198	1.219	1.195	1.180	1.198	1.203	1.207	1.201	1.196	1.199	1.212
1.245	1.252	1.226	1.237	1.229	1.212	1.237	1.229	1.215	1.203	1.218	1.204	1.176	1.181	1.176	1.215	1.196	1.189	1.189	1.212
																			
1.422	1.360	1.338	1.399	1.350	1.331	1.350	1.360	1.332	1.329	1.322	1.329	1.332	1.336	1.313	1.313	1.335	1.317	1.331	1.329
1.505	1.467	1.445	1.441	1.413	1.413	1.425	1.439	1.427	1.440	1.443	1.409	1.385	1.380	1.368	1.368	1.361	1.369	1.354	1.383
1.574	1.525	1.529	1.496	1.462	1.476	1.491	1.475	1.478	1.501	1.473	1.469	1.464	1.436	1.468	1.419	1.422	1.412	1.417	1.447
1.594	1.596	1.552	1.532	1.469	1.463	1.482	1.472	1.483	1.473	1.481	1.474	1.441	1.445	1.442	1.426	1.411	1.403	1.411	1.449
1.567	1.583	1.544	1.557	1.470	1.460	1.474	1.475	1.464	1.459	1.472	1.456	1.429	1.447	1.447	1.443	1.443	1.474	1.444	1.425
1.576	1.548	1.536	1.546	1.498	1.491	1.477	1.509	1.507	1.496	1.478	1.472	1.454	1.485	1.493	1.504	1.483	1.507	1.474	1.470
1.630	1.577	1.574	1.551	1.543	1.492	1.525	1.527	1.500	1.514	1.498	1.523	1.470	1.491	1.508	1.517	1.509	1.536	1.527	1.502
1.653	1.566	1.592	1.548	1.513	1.519	1.554	1.534	1.522	1.534	1.530	1.536	1.511	1.503	1.517	1.533	1.539	1.548	1.556	1.520
1.613	1.595	1.557	1.532	1.561	1.530	1.543	1.546	1.539	1.549	1.544	1.540	1.546	1.537	1.568	1.533	1.539	1.560	1.545	1.532
1.602	1.609	1.556	1.510	1.546	1.496	1.514	1.519	1.539	1.531	1.518	1.505	1.477	1.506	1.535	1.499	1.513	1.534	1.531	1.516
1.587	1.606	1.570	1.521	1.534	1.546	1.565	1.531	1.555	1.559	1.537	1.540	1.492	1.498	1.508	1.524	1.501	1.534	1.548	1.505
1.588	1.589	1.576	1.553	1.537	1.523	1.556	1.551	1.548	1.558	1.548	1.535	1.502	1.469	1.533	1.537	1.509	1.534	1.560	1.537
1.602	1.583	1.591	1.564	1.548	1.552	1.556	1.561	1.565	1.564	1.575	1.548	1.519	1.501	1.541	1.542	1.528	1.580	1.568	1.556
1.618	1.583	1.577	1.559	1.526	1.544	1.583	1.571	1.561	1.549	1.554	1.540	1.503	1.526	1.530	1.519	1.502	1.547	1.562	1.548
1.639	1.589	1.573	1.566	1.562	1.568	1.584	1.567	1.536	1.554	1.552	1.559	1.516	1.503	1.537	1.542	1.507	1.532	1.532	1.542
1.610	1.615	1.580	1.592	1.563	1.547	1.582	1.567	1.550	1.560	1.565	1.530	1.540	1.537	1.528	1.565	1.533	1.580	1.542	1.558
1.618	1.634	1.598	1.565	1.574	1.582	1.584	1.553	1.578	1.580	1.559	1.569	1.534	1.551	1.555	1.569	1.551	1.579	1.575	1.559
1.603	1.605	1.585	1.565	1.546	1.559	1.573	1.564	1.572	1.584	1.574	1.572	1.532	1.538	1.571	1.594	1.542	1.570	1.576	1.571
1.593	1.606	1.582	1.552	1.555	1.565	1.600	1.563	1.559	1.574	1.570	1.561	1.545	1.540	1.558	1.569	1.543	1.567	1.580	1.567
1.578	1.609	1.603	1.576	1.579	1.571	1.576	1.557	1.561	1.564	1.580	1.589	1.577	1.560	1.583	1.563	1.550	1.586	1.590	1.594
																			
2.102	2.030	1.997	2.006	2.054	2.052	2.036	2.033	2.030	2.077	2.082	2.100	2.120	2.077	2.077	2.074	2.097	2.087	2.057	2.044
2.187	2.214	2.073	2.089	2.178	2.195	2.169	2.151	2.164	2.135	2.115	2.137	2.119	2.086	2.042	2.040	2.030	2.034	2.066	2.011
2.155	2.144	2.203	2.180	2.146	2.143	2.151	2.127	2.099	2.148	2.134	2.137	2.130	2.086	2.072	2.076	2.096	2.147	2.166	2.076
2.230	2.200	2.273	2.232	2.158	2.137	2.181	2.195	2.158	2.178	2.183	2.118	2.098	2.139	2.139	2.182	2.189	2.166	2.208	2.227
2.233	2.215	2.245	2.162	2.102	2.172	2.163	2.204	2.168	2.136	2.179	2.125	2.132	2.109	2.105	2.134	2.146	2.144	2.151	2.196
2.265	2.240	2.216	2.188	2.190	2.212	2.218	2.198	2.182	2.182	2.207	2.187	2.189	2.186	2.164	2.163	2.162	2.194	2.192	2.207
2.267	2.281	2.200	2.138	2.152	2.184	2.216	2.212	2.143	2.168	2.206	2.175	2.145	2.139	2.154	2.177	2.225	2.237	2.236	2.204
2.253	2.284	2.231	2.210	2.106	2.156	2.199	2.195	2.181	2.182	2.225	2.186	2.125	2.197	2.171	2.191	2.217	2.249	2.230	2.202
2.243	2.296	2.194	2.248	2.142	2.224	2.202	2.223	2.227	2.243	2.212	2.188	2.191	2.186	2.227	2.241	2.179	2.274	2.251	2.213
2.274	2.231	2.220	2.223	2.143	2.146	2.197	2.232	2.233	2.251	2.243	2.200	2.308	2.232	2.294	2.264	2.264	2.338	2.268	2.234
2.280	2.267	2.228	2.189	2.143	2.130	2.202	2.213	2.223	2.248	2.327	2.263	2.285	2.298	2.277	2.263	2.264	2.308	2.297	2.234
2.294	2.256	2.168	2.142	2.133	2.126	2.168	2.189	2.189	2.175	2.240	2.238	2.291	2.273	2.260	2.253	2.257	2.270	2.285	2.266
2.294	2.289	2.177	2.205	2.170	2.151	2.178	2.176	2.188	2.241	2.241	2.238	2.296	2.263	2.280	2.223	2.290	2.300	2.264	2.284
2.250	2.259	2.198	2.247	2.172	2.157	2.173	2.151	2.182	2.204	2.162	2.208	2.237	2.223	2.314	2.226	2.294	2.314	2.291	2.299
2.256	2.276	2.208	2.172	2.158	2.159	2.204	2.166	2.203	2.252	2.205	2.230	2.240	2.240	2.277	2.177	2.230	2.264	2.273	2.234
2.293	2.273	2.202	2.206	2.177	2.184	2.217	2.195	2.243	2.232	2.234	2.251	2.259	2.265	2.250	2.251	2.289	2.315	2.273	2.267
2.288	2.254	2.204	2.220	2.226	2.212	2.187	2.208	2.215	2.270	2.238	2.242	2.246	2.268	2.210	2.220	2.243	2.270	2.272	2.256
2.284	2.281	2.266	2.261	2.216	2.219	2.209	2.225	2.263	2.255	2.250	2.229	2.253	2.230	2.231	2.257	2.246	2.281	2.216	2.249
2.275	2.257	2.263	2.242	2.213	2.225	2.211	2.237	2.224	2.249	2.230	2.187	2.219	2.183	2.239	2.220	2.234	2.273	2.227	2.257
2.281	2.286	2.274	2.236	2.193	2.230	2.217	2.234	2.297	2.275	2.343	2.310	2.234	2.251	2.292	2.248	2.230	2.289	2.238	2.316

fixed vec[rec] reccv_encf(3)(20,20)
0.335	0.473	0.685	0.791	0.925	0.984	1.076	1.137	1.208	1.253	1.280	1.442	1.540	1.609	1.710	1.807	1.890	1.904	1.904	2.132
0.524	0.716	1.019	1.228	1.371	1.471	1.573	1.694	1.791	1.893	1.914	2.074	2.268	2.428	2.523	2.650	2.723	2.797	2.788	3.224
0.686	0.890	1.285	1.532	1.718	1.905	2.097	2.183	2.238	2.285	2.366	2.664	2.900	3.132	3.186	3.309	3.451	3.518	3.530	4.054
0.776	1.062	1.528	1.789	1.951	2.169	2.307	2.399	2.542	2.635	2.727	3.032	3.296	3.513	3.654	3.739	3.882	3.948	3.993	4.632
0.864	1.198	1.639	1.970	2.166	2.346	2.561	2.660	2.700	2.921	3.033	3.289	3.622	3.863	4.010	4.072	4.232	4.358	4.438	5.030
0.934	1.321	1.800	2.161	2.376	2.632	2.794	2.891	3.095	3.250	3.389	3.678	3.895	4.262	4.491	4.563	4.794	4.856	4.994	5.619
1.024	1.439	2.026	2.398	2.669	2.912	3.074	3.220	3.403	3.468	3.662	4.075	4.372	4.696	4.962	5.071	5.185	5.239	5.327	6.077
1.115	1.535	2.144	2.544	2.865	3.073	3.339	3.474	3.696	3.805	3.998	4.467	4.700	4.976	5.165	5.386	5.512	5.641	5.755	6.515
1.171	1.650	2.299	2.693	3.072	3.271	3.576	3.769	3.925	4.041	4.328	4.797	5.098	5.415	5.624	5.884	6.090	6.243	6.259	7.110
1.250	1.749	2.435	2.834	3.174	3.413	3.715	3.914	4.035	4.247	4.388	4.900	5.197	5.521	5.727	6.007	6.278	6.420	6.406	7.374
1.332	1.847	2.563	2.996	3.279	3.573	3.863	4.111	4.290	4.450	4.561	5.171	5.567	5.828	6.032	6.392	6.582	6.689	6.767	7.808
1.405	1.958	2.688	3.071	3.487	3.743	4.033	4.290	4.472	4.638	4.725	5.333	5.742	6.089	6.381	6.717	6.890	7.068	7.124	8.024
1.454	2.023	2.805	3.248	3.671	3.937	4.220	4.510	4.699	4.927	5.051	5.661	6.012	6.394	6.650	7.007	7.213	7.424	7.544	8.517
1.498	2.050	2.870	3.353	3.824	4.031	4.358	4.695	4.888	5.036	5.165	5.878	6.237	6.592	6.959	7.267	7.364	7.617	7.821	8.628
1.553	2.120	2.950	3.434	3.910	4.133	4.490	4.776	4.944	5.139	5.361	6.005	6.342	6.877	7.107	7.457	7.614	7.774	8.031	8.980
1.612	2.241	3.052	3.623	4.084	4.379	4.761	4.942	5.151	5.427	5.550	6.177	6.598	6.979	7.445	7.638	7.795	8.045	8.280	9.155
1.644	2.300	3.199	3.730	4.246	4.536	4.906	5.177	5.369	5.561	5.701	6.443	6.798	7.148	7.701	8.005	8.073	8.286	8.533	9.700
1.688	2.386	3.237	3.818	4.336	4.633	5.063	5.364	5.491	5.719	5.840	6.688	7.010	7.441	7.693	8.198	8.473	8.651	8.814	10.013
1.749	2.439	3.329	3.887	4.425	4.686	5.189	5.492	5.729	5.840	6.072	6.726	7.255	7.666	8.097	8.473	8.690	8.829	8.974	10.198
1.780	2.501	3.421	4.080	4.566	4.889	5.339	5.605	5.853	6.017	6.259	6.933	7.403	7.879	8.132	8.610	8.715	8.857	9.110	10.455
																			
0.520	0.744	1.079	1.312	1.467	1.584	1.724	1.845	1.967	2.084	2.085	2.374	2.561	2.685	2.802	2.920	3.068	3.129	3.189	3.540
0.766	1.028	1.481	1.779	2.016	2.234	2.510	2.580	2.663	2.785	2.889	3.293	3.426	3.627	3.795	3.827	4.025	4.191	4.211	4.832
0.940	1.273	1.865	2.256	2.442	2.709	2.909	3.065	3.250	3.463	3.564	3.989	4.213	4.384	4.722	4.713	4.869	4.946	5.086	5.762
1.060	1.526	2.181	2.576	2.743	3.007	3.313	3.463	3.636	3.769	3.894	4.542	4.772	4.957	5.243	5.338	5.408	5.650	5.735	6.670
1.186	1.671	2.349	2.788	2.972	3.283	3.592	3.808	4.009	4.121	4.307	4.810	5.025	5.469	5.745	5.939	6.094	6.208	6.372	7.116
1.281	1.790	2.514	3.008	3.260	3.586	3.966	4.236	4.407	4.524	4.701	5.277	5.603	6.157	6.369	6.580	6.878	7.063	7.120	7.937
1.417	1.936	2.768	3.280	3.606	3.892	4.270	4.617	4.788	4.910	5.103	5.751	6.115	6.455	6.878	7.068	7.455	7.711	7.770	8.653
1.518	2.121	2.988	3.470	3.821	4.226	4.594	4.842	4.999	5.231	5.289	6.140	6.623	6.942	7.133	7.524	7.920	8.219	8.376	9.223
1.584	2.300	3.100	3.667	4.144	4.555	4.886	5.159	5.435	5.556	5.630	6.545	7.096	7.477	7.912	8.037	8.352	8.633	8.635	9.764
1.695	2.394	3.218	3.729	4.294	4.576	4.949	5.283	5.522	5.728	5.923	6.634	7.136	7.587	8.058	8.215	8.570	8.974	9.040	10.129
1.733	2.521	3.357	3.930	4.506	4.839	5.342	5.532	5.830	6.143	6.196	7.037	7.402	7.763	8.223	8.599	8.933	9.213	9.397	10.436
1.803	2.566	3.476	4.104	4.649	5.088	5.486	5.747	6.043	6.376	6.502	7.355	7.860	8.090	8.887	9.103	9.323	9.652	9.942	10.876
1.895	2.646	3.681	4.308	4.819	5.334	5.797	6.150	6.422	6.578	6.796	7.563	8.015	8.554	9.156	9.480	9.870	10.251	10.484	11.689
2.000	2.748	3.770	4.451	5.037	5.455	6.055	6.336	6.562	6.852	7.027	7.845	8.490	8.941	9.372	9.655	9.933	10.419	10.585	12.000
2.057	2.826	3.872	4.648	5.114	5.727	6.120	6.406	6.708	6.888	7.200	8.220	8.608	9.129	9.555	10.019	10.216	10.834	10.699	12.236
2.104	2.970	4.113	4.855	5.345	5.864	6.341	6.639	6.900	7.167	7.480	8.398	9.025	9.469	9.962	10.427	10.709	11.186	11.143	12.691
2.182	3.110	4.247	4.947	5.590	6.108	6.637	6.885	7.146	7.418	7.723	8.730	9.202	9.781	10.345	10.812	11.002	11.508	11.558	13.151
2.231	3.144	4.276	5.099	5.678	6.156	6.756	7.038	7.393	7.647	7.986	9.065	9.633	9.985	10.604	11.182	11.185	11.646	12.008	13.311
2.266	3.239	4.369	5.192	5.812	6.313	6.894	7.175	7.596	7.825	8.095	9.252	9.894	10.289	10.938	11.256	11.411	11.940	12.200	13.696
2.318	3.329	4.570	5.342	6.084	6.543	7.046	7.304	7.766	8.120	8.467	9.420	10.112	10.621	11.267	11.488	11.852	12.326	12.509	14.382
																			
0.954	1.397	2.098	2.662	2.904	3.209	3.436	3.611	3.722	3.920	4.134	4.769	5.087	5.367	5.521	5.696	5.797	5.924	6.072	6.558
1.258	1.854	2.604	3.172	3.543	4.102	4.446	4.647	4.705	4.971	5.107	5.671	6.282	6.479	6.539	6.593	6.812	7.111	7.212	8.160
1.529	2.115	3.098	3.658	4.128	4.574	4.936	5.221	5.179	5.484	5.805	6.699	7.074	7.401	7.960	8.090	8.536	8.799	8.964	9.764
1.711	2.479	3.462	4.060	4.613	5.066	5.435	5.690	6.002	6.221	6.507	7.474	8.064	8.475	8.892	9.296	9.709	10.207	9.948	11.335
1.870	2.654	3.651	4.398	4.798	5.517	5.808	6.352	6.667	6.870	7.074	7.977	8.396	8.965	9.489	9.670	10.508	10.532	10.655	12.218
2.053	2.866	3.976	4.699	5.416	5.957	6.526	6.868	7.029	7.321	7.309	8.735	9.226	9.835	10.212	10.699	11.160	11.544	11.862	13.001
2.269	3.047	4.203	4.783	5.813	6.170	6.850	7.249	7.611	7.922	8.283	9.207	9.541	10.380	11.055	11.475	11.745	12.373	12.680	14.218
2.289	3.109	4.356	5.327	5.915	6.529	7.037	7.700	7.985	8.385	8.577	9.658	10.428	11.006	11.649	12.177	12.590	13.141	13.330	14.788
2.473	3.311	4.747	5.625	6.116	6.959	7.601	8.279	8.544	8.949	9.318	10.143	10.871	11.601	12.516	12.956	13.031	13.586	13.942	15.678
2.548	3.442	4.758	5.815	6.433	7.224	8.076	8.450	8.832	9.367	9.928	10.850	11.944	12.571	13.313	13.662	13.794	14.627	14.502	16.102
2.666	3.624	5.080	6.000	6.841	7.677	8.451	8.863	9.253	9.897	10.286	11.487	12.268	13.278	13.554	14.055	14.803	15.321	15.260	16.832
2.762	3.732	5.165	6.153	6.965	7.772	8.658	8.961	9.305	10.049	10.504	11.698	12.970	13.369	13.752	14.563	14.780	15.778	16.092	17.699
2.846	3.994	5.371	6.601	7.287	7.911	8.832	9.100	9.856	10.080	10.701	12.265	13.162	13.837	14.456	14.941	15.469	16.174	16.600	18.334
2.916	4.181	5.745	6.919	7.507	8.225	8.745	9.424	10.271	10.385	10.935	12.458	13.065	14.077	14.789	15.059	15.905	16.227	16.809	19.431
2.995	4.224	5.783	6.911	7.910	8.590	9.269	9.717	10.405	10.829	11.261	12.828	13.697	14.461	14.988	15.182	15.855	16.440	17.237	19.343
3.194	4.377	6.065	7.172	7.993	8.926	9.495	10.094	10.462	11.290	11.637	13.182	14.004	14.676	15.434	16.128	16.854	17.705	17.898	19.783
3.283	4.506	6.273	7.245	8.516	9.177	9.841	10.364	10.953	11.552	11.944	13.844	14.739	15.454	16.072	16.316	17.110	17.701	18.128	20.237
3.341	4.645	6.531	7.620	8.447	9.450	9.913	10.620	11.198	11.903	12.480	13.924	15.166	15.530	15.998	16.685	17.333	17.754	17.650	20.556
3.385	4.642	6.617	7.818	8.578	9.648	10.205	10.591	11.296	12.021	12.667	13.954	15.287	15.712	16.367	16.658	17.431	17.725	18.102	21.022
3.468	4.859	6.800	7.940	8.970	10.084	10.627	11.511	12.056	12.629	13.342	14.745	15.554	16.332	17.310	17.455	18.091	18.618	18.919	21.565

fixed vec[rec] rollcv_mset(3)(20,20)
0.903	0.776	0.637	0.530	0.401	0.317	0.246	0.180	0.136	0.116	0.078	-0.068	-0.191	-0.292	-0.386	-0.491	-0.548	-0.627	-0.721	-1.336
0.915	0.755	0.569	0.425	0.280	0.155	0.111	0.026	-0.050	-0.094	-0.140	-0.357	-0.534	-0.691	-0.844	-0.966	-1.088	-1.222	-1.324	-2.236
0.938	0.718	0.521	0.346	0.201	0.064	-0.042	-0.137	-0.224	-0.302	-0.346	-0.593	-0.857	-1.043	-1.214	-1.387	-1.523	-1.662	-1.812	-2.909
0.866	0.731	0.494	0.248	0.098	-0.047	-0.164	-0.262	-0.362	-0.434	-0.505	-0.802	-1.090	-1.305	-1.517	-1.701	-1.883	-2.056	-2.197	-3.506
0.825	0.694	0.402	0.154	-0.025	-0.168	-0.305	-0.399	-0.508	-0.589	-0.674	-1.008	-1.305	-1.561	-1.807	-2.023	-2.212	-2.419	-2.592	-4.032
0.811	0.602	0.319	0.088	-0.095	-0.262	-0.423	-0.523	-0.638	-0.732	-0.821	-1.184	-1.527	-1.808	-2.072	-2.299	-2.517	-2.727	-2.931	-4.518
0.791	0.573	0.279	0.038	-0.157	-0.326	-0.497	-0.611	-0.750	-0.841	-0.933	-1.338	-1.692	-2.034	-2.292	-2.547	-2.775	-3.022	-3.236	-4.961
0.758	0.541	0.244	-0.042	-0.244	-0.408	-0.576	-0.727	-0.838	-0.957	-1.040	-1.482	-1.863	-2.219	-2.518	-2.774	-3.016	-3.282	-3.490	-5.370
0.742	0.520	0.193	-0.105	-0.322	-0.491	-0.657	-0.803	-0.951	-1.049	-1.153	-1.633	-2.045	-2.409	-2.723	-3.003	-3.267	-3.525	-3.776	-5.751
0.727	0.500	0.138	-0.144	-0.374	-0.568	-0.757	-0.902	-1.045	-1.167	-1.288	-1.791	-2.197	-2.573	-2.915	-3.203	-3.485	-3.769	-4.038	-6.103
0.701	0.474	0.094	-0.217	-0.434	-0.626	-0.849	-0.990	-1.128	-1.267	-1.385	-1.935	-2.362	-2.754	-3.099	-3.388	-3.726	-4.016	-4.285	-6.454
0.679	0.453	0.079	-0.276	-0.516	-0.713	-0.914	-1.083	-1.229	-1.363	-1.480	-2.042	-2.517	-2.923	-3.273	-3.590	-3.923	-4.231	-4.506	-6.800
0.674	0.421	0.032	-0.317	-0.564	-0.760	-0.966	-1.144	-1.305	-1.454	-1.575	-2.159	-2.654	-3.070	-3.455	-3.799	-4.129	-4.443	-4.739	-7.109
0.658	0.394	-0.006	-0.337	-0.620	-0.818	-1.039	-1.243	-1.403	-1.537	-1.676	-2.293	-2.791	-3.238	-3.621	-3.978	-4.320	-4.647	-4.957	-7.407
0.628	0.356	-0.059	-0.401	-0.678	-0.888	-1.127	-1.315	-1.488	-1.633	-1.770	-2.405	-2.937	-3.377	-3.806	-4.154	-4.513	-4.846	-5.158	-7.715
0.609	0.331	-0.096	-0.432	-0.720	-0.944	-1.172	-1.390	-1.550	-1.724	-1.874	-2.492	-3.058	-3.519	-3.954	-4.307	-4.690	-5.043	-5.349	-7.985
0.574	0.296	-0.130	-0.487	-0.763	-1.008	-1.222	-1.462	-1.620	-1.798	-1.938	-2.604	-3.179	-3.633	-4.086	-4.462	-4.855	-5.228	-5.549	-8.264
0.565	0.272	-0.172	-0.541	-0.829	-1.089	-1.296	-1.524	-1.695	-1.878	-2.021	-2.713	-3.292	-3.768	-4.233	-4.597	-5.018	-5.392	-5.730	-8.534
0.552	0.248	-0.198	-0.581	-0.895	-1.159	-1.360	-1.585	-1.762	-1.952	-2.097	-2.810	-3.413	-3.893	-4.365	-4.774	-5.176	-5.580	-5.915	-8.788
0.519	0.206	-0.232	-0.642	-0.938	-1.202	-1.419	-1.643	-1.827	-2.025	-2.174	-2.911	-3.509	-4.012	-4.510	-4.918	-5.345	-5.744	-6.105	-9.040
																			
1.251	1.117	0.970	0.859	0.722	0.651	0.575	0.510	0.455	0.382	0.334	0.171	0.087	-0.026	-0.129	-0.220	-0.291	-0.366	-0.442	-1.084
1.280	1.105	0.884	0.753	0.631	0.484	0.401	0.304	0.235	0.166	0.103	-0.073	-0.279	-0.437	-0.596	-0.707	-0.821	-0.957	-1.066	-1.971
1.274	1.088	0.842	0.667	0.490	0.381	0.251	0.146	0.066	-0.016	-0.084	-0.331	-0.581	-0.761	-0.939	-1.111	-1.242	-1.398	-1.535	-2.639
1.267	1.087	0.852	0.585	0.376	0.274	0.136	0.024	-0.080	-0.173	-0.222	-0.525	-0.794	-1.027	-1.254	-1.423	-1.590	-1.777	-1.936	-3.236
1.229	1.034	0.716	0.479	0.280	0.155	-0.019	-0.090	-0.219	-0.329	-0.385	-0.737	-1.033	-1.303	-1.550	-1.771	-1.935	-2.117	-2.300	-3.748
1.176	0.966	0.621	0.407	0.225	0.058	-0.119	-0.218	-0.336	-0.428	-0.535	-0.903	-1.236	-1.524	-1.813	-2.022	-2.213	-2.454	-2.662	-4.234
1.154	0.936	0.628	0.346	0.171	-0.011	-0.182	-0.320	-0.433	-0.531	-0.663	-1.057	-1.423	-1.773	-2.060	-2.275	-2.501	-2.748	-2.941	-4.676
1.157	0.924	0.562	0.258	0.081	-0.099	-0.281	-0.432	-0.552	-0.672	-0.785	-1.235	-1.602	-1.952	-2.227	-2.489	-2.755	-3.004	-3.239	-5.114
1.117	0.892	0.504	0.213	0.021	-0.156	-0.374	-0.529	-0.623	-0.785	-0.885	-1.369	-1.762	-2.131	-2.455	-2.704	-2.981	-3.262	-3.487	-5.471
1.068	0.872	0.443	0.133	-0.038	-0.258	-0.466	-0.605	-0.765	-0.909	-1.011	-1.518	-1.936	-2.305	-2.655	-2.914	-3.217	-3.490	-3.752	-5.841
1.042	0.823	0.400	0.058	-0.137	-0.340	-0.533	-0.689	-0.846	-0.986	-1.121	-1.642	-2.114	-2.473	-2.820	-3.115	-3.426	-3.726	-4.015	-6.169
1.026	0.799	0.362	0.022	-0.231	-0.426	-0.612	-0.768	-0.935	-1.089	-1.237	-1.786	-2.241	-2.630	-3.022	-3.340	-3.639	-3.948	-4.248	-6.514
0.999	0.781	0.346	-0.003	-0.240	-0.483	-0.686	-0.842	-1.036	-1.180	-1.320	-1.882	-2.365	-2.814	-3.171	-3.499	-3.822	-4.140	-4.464	-6.854
0.986	0.720	0.320	-0.039	-0.307	-0.535	-0.751	-0.911	-1.104	-1.267	-1.413	-1.981	-2.521	-2.983	-3.354	-3.685	-4.039	-4.356	-4.684	-7.130
0.957	0.708	0.254	-0.069	-0.346	-0.613	-0.806	-1.015	-1.183	-1.350	-1.494	-2.073	-2.658	-3.105	-3.514	-3.874	-4.219	-4.551	-4.894	-7.425
0.964	0.692	0.216	-0.097	-0.398	-0.663	-0.868	-1.098	-1.281	-1.434	-1.594	-2.191	-2.772	-3.255	-3.658	-4.027	-4.397	-4.746	-5.092	-7.699
0.933	0.667	0.189	-0.154	-0.448	-0.711	-0.946	-1.156	-1.341	-1.525	-1.634	-2.316	-2.878	-3.377	-3.797	-4.183	-4.585	-4.922	-5.275	-8.002
0.910	0.608	0.136	-0.219	-0.534	-0.783	-1.011	-1.232	-1.420	-1.608	-1.744	-2.409	-3.011	-3.491	-3.930	-4.343	-4.746	-5.111	-5.462	-8.239
0.896	0.580	0.083	-0.265	-0.596	-0.860	-1.068	-1.286	-1.484	-1.662	-1.824	-2.486	-3.115	-3.631	-4.082	-4.487	-4.920	-5.295	-5.661	-8.502
0.864	0.562	0.075	-0.313	-0.640	-0.904	-1.120	-1.351	-1.532	-1.747	-1.878	-2.589	-3.227	-3.750	-4.247	-4.668	-5.093	-5.483	-5.845	-8.759
																			
1.875	1.799	1.604	1.447	1.340	1.221	1.179	1.098	1.021	0.969	0.882	0.735	0.505	0.448	0.346	0.260	0.173	0.104	0.034	-0.622
1.959	1.757	1.504	1.325	1.180	1.165	0.996	0.953	0.883	0.744	0.640	0.392	0.203	0.004	-0.158	-0.279	-0.401	-0.490	-0.577	-1.509
1.860	1.669	1.473	1.271	1.076	0.984	0.896	0.773	0.614	0.504	0.431	0.157	-0.052	-0.311	-0.516	-0.624	-0.802	-0.964	-1.055	-2.169
1.905	1.700	1.503	1.183	1.003	0.903	0.755	0.656	0.455	0.342	0.234	-0.005	-0.262	-0.546	-0.805	-0.935	-1.117	-1.276	-1.465	-2.735
1.881	1.627	1.347	1.112	0.927	0.790	0.657	0.504	0.307	0.193	0.123	-0.223	-0.550	-0.818	-1.041	-1.203	-1.383	-1.622	-1.786	-3.245
1.826	1.680	1.312	1.007	0.850	0.641	0.558	0.336	0.195	0.069	0.017	-0.418	-0.708	-1.036	-1.310	-1.509	-1.698	-1.913	-2.083	-3.708
1.842	1.620	1.233	0.989	0.751	0.526	0.485	0.227	0.055	-0.039	-0.127	-0.582	-0.873	-1.210	-1.546	-1.807	-1.992	-2.239	-2.398	-4.175
1.819	1.582	1.178	0.918	0.702	0.466	0.349	0.132	-0.018	-0.176	-0.302	-0.722	-1.062	-1.436	-1.745	-1.977	-2.250	-2.428	-2.678	-4.571
1.768	1.510	1.110	0.845	0.600	0.408	0.235	0.036	-0.099	-0.277	-0.407	-0.834	-1.265	-1.623	-1.902	-2.194	-2.381	-2.704	-2.946	-4.970
1.713	1.428	1.075	0.808	0.536	0.298	0.122	-0.064	-0.248	-0.381	-0.482	-0.994	-1.416	-1.801	-2.128	-2.412	-2.669	-2.952	-3.239	-5.359
1.709	1.389	1.001	0.763	0.451	0.198	0.049	-0.142	-0.352	-0.484	-0.587	-1.084	-1.534	-1.942	-2.309	-2.625	-2.868	-3.187	-3.435	-5.652
1.647	1.377	0.944	0.660	0.388	0.121	-0.072	-0.284	-0.468	-0.563	-0.680	-1.220	-1.731	-2.119	-2.487	-2.774	-3.093	-3.389	-3.669	-5.974
1.632	1.393	0.934	0.640	0.352	0.047	-0.128	-0.387	-0.487	-0.662	-0.815	-1.367	-1.836	-2.229	-2.695	-2.977	-3.350	-3.629	-3.891	-6.272
1.597	1.357	0.912	0.569	0.272	-0.065	-0.200	-0.459	-0.599	-0.752	-0.840	-1.481	-1.947	-2.417	-2.832	-3.164	-3.526	-3.859	-4.161	-6.567
1.633	1.280	0.845	0.517	0.143	-0.141	-0.322	-0.524	-0.672	-0.831	-0.984	-1.594	-2.077	-2.535	-2.961	-3.319	-3.675	-4.050	-4.377	-6.896
1.592	1.368	0.838	0.491	0.162	-0.136	-0.361	-0.550	-0.781	-0.918	-1.045	-1.687	-2.285	-2.700	-3.117	-3.485	-3.867	-4.267	-4.618	-7.165
1.581	1.281	0.794	0.435	0.134	-0.191	-0.407	-0.646	-0.837	-1.018	-1.128	-1.737	-2.384	-2.824	-3.285	-3.659	-4.039	-4.435	-4.809	-7.448
1.595	1.248	0.791	0.377	0.058	-0.227	-0.442	-0.716	-0.894	-1.080	-1.207	-1.844	-2.459	-2.925	-3.438	-3.830	-4.221	-4.603	-4.979	-7.729
1.536	1.261	0.759	0.349	-0.020	-0.285	-0.474	-0.799	-0.944	-1.171	-1.294	-1.905	-2.593	-3.033	-3.561	-3.990	-4.403	-4.816	-5.191	-7.990
1.519	1.267	0.689	0.259	-0.027	-0.323	-0.579	-0.823	-1.006	-1.203	-1.367	-2.026	-2.722	-3.164	-3.671	-4.198	-4.553	-4.989	-5.374	-8.241

fixed vec[rec] rollcv_msef(3)(20,20)
0.536	0.667	0.838	0.865	0.773	0.693	0.602	0.482	0.390	0.355	0.276	-0.286	-0.875	-1.506	-2.246	-2.949	-3.532	-4.272	-5.187	-13.049
0.812	0.920	1.028	1.004	0.806	0.468	0.399	0.095	-0.198	-0.394	-0.623	-1.993	-3.447	-4.913	-6.572	-8.043	-9.725	-11.325	-13.113	-30.130
1.006	1.074	1.135	0.944	0.617	0.224	-0.174	-0.600	-1.080	-1.529	-1.847	-3.914	-6.576	-8.888	-11.323	-13.863	-16.194	-18.839	-21.416	-47.629
1.096	1.284	1.263	0.777	0.356	-0.200	-0.788	-1.341	-1.973	-2.528	-3.182	-6.127	-9.485	-12.721	-16.354	-19.491	-22.936	-26.582	-29.796	-66.046
1.152	1.315	1.143	0.541	-0.114	-0.774	-1.583	-2.300	-3.102	-3.896	-4.649	-8.585	-12.784	-17.001	-21.287	-25.763	-29.963	-34.596	-38.747	-84.385
1.224	1.311	0.981	0.365	-0.428	-1.348	-2.392	-3.270	-4.292	-5.180	-6.114	-11.040	-16.077	-21.463	-26.638	-31.806	-36.862	-42.561	-47.968	-103.064
1.292	1.355	0.939	0.151	-0.783	-1.833	-3.061	-4.153	-5.487	-6.474	-7.583	-13.380	-19.388	-25.585	-31.927	-37.903	-44.106	-50.467	-56.648	-121.808
1.386	1.357	0.886	-0.186	-1.278	-2.492	-3.824	-5.225	-6.556	-7.690	-8.939	-15.798	-22.643	-29.690	-37.081	-43.992	-51.114	-58.675	-65.549	-140.949
1.378	1.433	0.772	-0.521	-1.835	-3.139	-4.586	-6.133	-7.734	-9.173	-10.558	-18.228	-25.996	-34.174	-42.471	-50.545	-58.373	-66.600	-74.992	-159.911
1.456	1.411	0.555	-0.701	-2.189	-3.790	-5.566	-7.200	-9.046	-10.574	-12.294	-20.945	-29.770	-38.535	-47.791	-56.336	-65.678	-74.860	-84.236	-178.525
1.478	1.421	0.387	-1.140	-2.708	-4.517	-6.566	-8.306	-10.226	-11.915	-13.979	-23.518	-33.168	-42.934	-53.072	-62.565	-73.171	-83.193	-93.524	-197.625
1.494	1.445	0.342	-1.529	-3.383	-5.221	-7.339	-9.269	-11.431	-13.445	-15.523	-25.920	-36.975	-47.635	-58.520	-69.103	-80.574	-91.819	-102.910 -216.499
1.514	1.384	0.146	-1.840	-3.827	-5.881	-8.052	-10.310	-12.608	-14.865	-17.096	-28.279	-40.075	-52.057	-64.057	-75.533	-87.883	-99.987	-112.424 -235.673
1.532	1.318	-0.028	-1.980	-4.337	-6.490	-9.031	-11.543	-13.998	-16.354	-18.750	-31.166	-43.912	-56.740	-69.583	-82.448	-95.317	-108.397 -121.537 -255.123
1.534	1.252	-0.295	-2.432	-4.813	-7.238	-9.981	-12.673	-15.258	-17.947	-20.461	-33.836	-47.598	-61.190	-74.968	-88.882	-102.787 -116.771 -130.799 -274.306
1.500	1.188	-0.524	-2.757	-5.289	-7.934	-10.866	-13.753	-16.616	-19.581	-22.149	-36.462	-51.297	-65.645	-80.608	-95.309	-110.186 -125.317 -140.108 -293.420
1.504	1.129	-0.708	-3.214	-5.896	-8.696	-11.686	-14.894	-17.915	-20.918	-23.751	-39.033	-54.748	-70.279	-86.177	-101.467 -117.455 -133.796 -149.606 -312.050
1.482	1.032	-0.938	-3.657	-6.589	-9.700	-12.668	-15.966	-19.140	-22.482	-25.541	-41.674	-58.276	-74.731	-91.453	-108.142 -124.755 -141.945 -158.876 -331.571
1.512	0.964	-1.146	-4.079	-7.372	-10.665	-13.670	-17.151	-20.447	-24.093	-27.151	-44.572	-61.951	-79.256	-96.868	-114.553 -132.496 -150.486 -168.448 -350.989
1.454	0.850	-1.339	-4.569	-7.847	-11.303	-14.548	-18.252	-21.907	-25.584	-28.859	-47.226	-65.687	-83.619	-102.646 -121.162 -139.867 -158.986 -178.143 -369.797
																			
0.854	1.112	1.394	1.644	1.627	1.583	1.574	1.469	1.488	1.378	1.215	0.742	0.451	-0.135	-0.751	-1.354	-1.979	-2.706	-3.269	-11.039
1.232	1.481	1.802	1.889	1.841	1.695	1.495	1.264	1.015	0.783	0.504	-0.433	-1.934	-3.249	-4.768	-6.154	-7.591	-9.109	-10.593	-27.333
1.472	1.752	2.089	2.042	1.700	1.532	1.100	0.694	0.340	-0.071	-0.471	-2.224	-4.623	-6.831	-9.187	-11.336	-13.760	-15.985	-18.475	-44.235
1.600	2.078	2.332	1.979	1.536	1.228	0.701	0.116	-0.491	-1.112	-1.487	-4.295	-7.122	-10.431	-13.807	-16.903	-19.937	-23.601	-26.618	-61.973
1.768	2.191	2.164	1.783	1.175	0.764	-0.121	-0.542	-1.454	-2.172	-2.765	-6.423	-10.403	-14.423	-18.647	-22.859	-26.860	-30.766	-35.403	-79.950
1.926	2.181	2.117	1.652	1.062	0.318	-0.745	-1.424	-2.302	-3.162	-4.256	-8.690	-13.437	-18.415	-23.869	-28.748	-33.687	-38.756	-44.256	-98.141
1.994	2.273	2.321	1.478	0.868	-0.076	-1.162	-2.243	-3.224	-4.261	-5.620	-10.816	-16.786	-22.722	-29.063	-34.702	-40.410	-46.749	-53.028	-116.555
2.080	2.436	2.169	1.213	0.468	-0.626	-1.885	-3.233	-4.430	-5.681	-6.991	-13.305	-19.983	-26.721	-33.692	-40.360	-47.271	-54.594	-61.548	-135.371
2.164	2.518	1.966	1.075	0.138	-1.113	-2.620	-4.036	-5.258	-6.931	-8.345	-15.623	-23.121	-31.169	-38.774	-46.314	-54.546	-62.531	-70.346	-153.722
2.224	2.520	1.893	0.730	-0.235	-1.733	-3.496	-4.940	-6.512	-8.255	-9.863	-18.293	-26.322	-35.339	-44.060	-52.562	-61.526	-70.511	-79.623	-172.501
2.220	2.518	1.782	0.312	-0.879	-2.421	-4.243	-5.979	-7.757	-9.390	-11.427	-20.679	-29.906	-39.596	-49.054	-58.604	-68.435	-78.363	-88.995	-190.741
2.288	2.524	1.633	0.129	-1.560	-3.180	-4.905	-6.884	-9.040	-11.048	-13.161	-23.176	-33.366	-44.056	-54.604	-65.017	-75.415	-86.306	-97.768	-209.587
2.338	2.580	1.637	-0.020	-1.702	-3.708	-5.835	-7.854	-10.080	-12.247	-14.431	-25.356	-36.822	-48.373	-59.927	-71.205	-82.686	-94.350	-106.812 -229.218
2.386	2.502	1.577	-0.248	-2.237	-4.319	-6.724	-8.794	-11.311	-13.616	-16.197	-27.947	-40.352	-52.948	-65.237	-77.257	-89.994	-103.149 -116.026 -248.078
2.414	2.520	1.297	-0.464	-2.575	-5.123	-7.502	-10.137	-12.559	-15.149	-17.730	-30.131	-44.012	-57.284	-70.579	-84.165	-97.584	-111.607 -125.715 -266.248
2.472	2.510	1.208	-0.639	-3.082	-5.699	-8.208	-11.309	-14.139	-16.645	-19.357	-32.531	-47.186	-61.597	-75.874	-89.979	-104.437 -119.873 -134.853 -284.690
2.470	2.589	1.038	-0.965	-3.512	-6.275	-9.143	-12.270	-15.219	-18.067	-20.661	-35.606	-50.779	-66.089	-80.988	-96.245	-112.175 -127.581 -143.619 -304.487
2.452	2.417	0.772	-1.494	-4.344	-7.101	-9.933	-13.393	-16.460	-19.623	-22.229	-38.051	-54.403	-70.123	-86.275	-102.292 -119.747 -136.277 -153.161 -323.628
2.472	2.336	0.487	-1.886	-4.910	-7.965	-10.801	-14.309	-17.628	-20.870	-23.920	-40.412	-57.844	-74.619	-91.742	-108.662 -126.747 -144.554 -162.849 -342.607
2.486	2.319	0.435	-2.259	-5.430	-8.603	-11.596	-15.314	-18.702	-22.402	-25.265	-43.162	-61.562	-78.663	-97.572	-115.704 -134.340 -153.434 -172.230 -361.550
																			
1.638	2.230	2.812	3.300	3.634	3.811	3.688	3.721	3.924	3.612	3.765	3.757	2.991	2.876	2.495	1.897	1.397	0.845	0.292	-6.776
2.036	2.819	3.544	3.988	4.066	4.398	4.403	4.109	4.293	4.046	3.566	2.678	1.612	0.035	-1.520	-2.671	-3.858	-5.144	-6.340	-21.685
2.476	3.128	4.135	4.120	4.264	4.519	4.386	4.123	3.373	3.089	2.685	1.242	-0.415	-2.895	-5.472	-6.941	-9.189	-11.489	-13.588	-37.884
2.724	3.649	4.474	4.586	4.432	4.459	4.296	3.621	2.905	2.337	1.699	-0.044	-2.587	-5.888	-9.515	-11.685	-14.543	-17.426	-20.787	-53.761
3.008	3.721	4.504	4.710	4.508	4.199	4.042	3.216	2.167	1.370	1.055	-2.272	-5.524	-9.801	-13.339	-16.874	-20.439	-24.461	-28.531	-70.501
3.214	4.093	4.532	4.786	4.456	3.899	3.473	2.324	1.500	0.615	0.159	-4.130	-8.171	-13.029	-17.869	-22.801	-26.979	-32.109	-36.623	-88.626
3.552	4.340	4.568	4.566	4.181	3.450	3.167	1.696	0.442	-0.300	-1.151	-6.148	-10.416	-16.838	-22.332	-28.014	-32.931	-38.460	-44.518	-106.942
3.480	4.537	4.707	4.681	4.041	3.065	2.488	1.013	-0.144	-1.562	-2.729	-8.479	-14.019	-21.144	-27.515	-33.224	-39.911	-45.815	-52.898	-123.822
3.552	4.438	4.711	4.443	3.754	2.622	1.723	0.327	-0.877	-2.543	-3.923	-10.505	-17.288	-25.036	-31.855	-38.396	-45.705	-53.420	-61.324	-141.699
3.728	4.413	4.815	4.589	3.460	2.145	1.121	-0.624	-2.313	-3.795	-5.166	-12.539	-20.858	-28.333	-37.185	-44.266	-52.839	-61.576	-70.365	-159.548
3.726	4.468	4.706	4.490	2.917	1.454	0.402	-1.201	-3.586	-4.815	-6.579	-14.559	-23.510	-32.534	-41.914	-50.393	-58.777	-69.487	-79.141	-179.414
3.936	4.657	4.388	3.829	2.635	0.880	-0.635	-2.625	-4.873	-5.983	-7.674	-16.826	-27.125	-36.439	-47.604	-56.157	-66.401	-76.566	-87.531	-197.327
3.994	4.930	4.569	3.860	2.296	0.367	-1.143	-3.888	-5.040	-7.231	-9.648	-19.187	-29.454	-39.692	-52.488	-62.133	-73.832	-84.790	-96.119	-214.979
4.064	4.882	4.893	3.532	2.091	-0.529	-1.944	-4.571	-6.493	-8.599	-10.499	-21.606	-32.571	-44.527	-57.372	-68.416	-80.484	-92.287	-105.412 -232.448
4.222	4.825	4.753	3.474	1.145	-1.231	-3.009	-5.436	-7.567	-9.942	-12.285	-23.910	-35.910	-48.704	-61.637	-73.995	-87.422	-102.298 -114.464 -252.097
4.522	5.132	4.567	3.368	1.287	-1.189	-3.634	-6.106	-8.655	-11.272	-13.375	-26.306	-39.420	-53.419	-66.881	-80.342	-94.306	-109.982 -124.542 -269.905
4.458	5.050	4.337	3.041	1.103	-1.868	-4.129	-7.166	-9.952	-12.363	-15.048	-27.802	-43.210	-56.350	-72.102	-86.878	-100.778 -117.425 -133.193 -287.873
4.632	5.148	4.505	2.921	0.472	-2.089	-4.441	-7.726	-10.965	-14.206	-16.328	-30.590	-46.097	-60.420	-77.783	-92.906	-108.023 -124.373 -143.035 -306.482
4.550	5.144	4.276	2.545	-0.172	-2.678	-5.065	-8.716	-11.804	-15.430	-17.949	-32.796	-49.673	-64.342	-82.202	-99.516	-115.815 -133.011 -151.264 -325.465
4.672	5.283	4.241	1.964	-0.278	-3.254	-5.984	-9.463	-12.514	-16.103	-19.498	-34.764	-53.196	-69.299	-87.519	-105.649 -122.830 -141.534 -161.212 -343.651

fixed vec[rec] rollcv_enct(3)(20,20)
1.070	1.028	1.019	0.999	0.951	0.949	0.948	0.964	0.953	0.965	0.992	1.032	1.058	1.079	1.098	1.115	1.147	1.145	1.136	1.178
1.150	1.091	1.067	1.083	1.064	1.022	1.044	1.059	1.029	1.053	1.073	1.124	1.118	1.149	1.161	1.167	1.189	1.186	1.189	1.208
1.223	1.121	1.104	1.110	1.082	1.089	1.086	1.067	1.062	1.075	1.089	1.138	1.137	1.160	1.186	1.196	1.224	1.221	1.201	1.251
1.203	1.184	1.185	1.104	1.100	1.100	1.092	1.081	1.103	1.123	1.117	1.179	1.170	1.194	1.215	1.220	1.252	1.225	1.232	1.248
1.177	1.201	1.160	1.125	1.063	1.106	1.105	1.127	1.110	1.119	1.127	1.175	1.197	1.218	1.213	1.229	1.239	1.246	1.247	1.260
1.214	1.170	1.153	1.149	1.140	1.119	1.117	1.137	1.125	1.115	1.161	1.197	1.169	1.222	1.226	1.227	1.239	1.236	1.245	1.243
1.218	1.201	1.198	1.187	1.171	1.141	1.132	1.152	1.153	1.167	1.161	1.233	1.217	1.205	1.227	1.244	1.228	1.230	1.233	1.262
1.232	1.198	1.177	1.150	1.168	1.175	1.163	1.187	1.169	1.177	1.184	1.225	1.229	1.213	1.244	1.253	1.252	1.238	1.256	1.256
1.221	1.228	1.204	1.186	1.185	1.216	1.183	1.198	1.190	1.192	1.198	1.214	1.229	1.207	1.252	1.257	1.263	1.271	1.265	1.271
1.240	1.241	1.206	1.168	1.190	1.181	1.169	1.183	1.170	1.179	1.190	1.196	1.227	1.221	1.241	1.281	1.279	1.269	1.260	1.261
1.239	1.257	1.221	1.180	1.207	1.201	1.158	1.191	1.169	1.170	1.197	1.199	1.216	1.230	1.251	1.278	1.256	1.245	1.248	1.267
1.241	1.275	1.224	1.182	1.176	1.186	1.190	1.187	1.193	1.185	1.201	1.200	1.207	1.233	1.242	1.280	1.256	1.251	1.253	1.253
1.251	1.270	1.246	1.184	1.203	1.224	1.207	1.196	1.188	1.188	1.192	1.208	1.218	1.241	1.246	1.277	1.254	1.263	1.260	1.254
1.249	1.264	1.231	1.202	1.197	1.218	1.187	1.193	1.176	1.174	1.191	1.212	1.212	1.224	1.252	1.258	1.268	1.251	1.243	1.276
1.259	1.259	1.232	1.215	1.215	1.223	1.201	1.191	1.181	1.176	1.194	1.224	1.210	1.225	1.249	1.266	1.261	1.274	1.265	1.263
1.239	1.270	1.230	1.215	1.228	1.233	1.223	1.183	1.174	1.182	1.188	1.234	1.219	1.218	1.258	1.270	1.251	1.266	1.269	1.259
1.246	1.263	1.248	1.229	1.225	1.217	1.221	1.207	1.188	1.194	1.213	1.215	1.216	1.250	1.250	1.275	1.253	1.254	1.256	1.253
1.235	1.255	1.245	1.203	1.200	1.191	1.224	1.201	1.197	1.193	1.203	1.235	1.235	1.268	1.268	1.288	1.266	1.264	1.264	1.248
1.254	1.255	1.255	1.224	1.201	1.184	1.228	1.203	1.218	1.213	1.212	1.218	1.234	1.269	1.251	1.287	1.268	1.229	1.263	1.249
1.237	1.250	1.243	1.230	1.233	1.208	1.217	1.198	1.231	1.209	1.230	1.222	1.232	1.295	1.257	1.288	1.259	1.245	1.247	1.239
																			
1.413	1.383	1.348	1.372	1.341	1.338	1.334	1.325	1.342	1.350	1.344	1.398	1.437	1.440	1.452	1.493	1.483	1.544	1.525	1.530
1.508	1.456	1.419	1.439	1.410	1.407	1.421	1.394	1.401	1.410	1.394	1.478	1.497	1.496	1.484	1.524	1.548	1.529	1.531	1.579
1.575	1.504	1.496	1.467	1.440	1.447	1.426	1.430	1.455	1.434	1.432	1.529	1.494	1.524	1.529	1.546	1.584	1.576	1.596	1.603
1.587	1.584	1.570	1.490	1.450	1.461	1.469	1.499	1.497	1.484	1.499	1.575	1.593	1.582	1.577	1.611	1.625	1.610	1.595	1.651
1.600	1.577	1.504	1.476	1.449	1.480	1.460	1.501	1.483	1.480	1.495	1.566	1.574	1.565	1.540	1.573	1.631	1.608	1.593	1.623
1.576	1.544	1.516	1.483	1.484	1.516	1.485	1.513	1.510	1.509	1.524	1.571	1.570	1.571	1.577	1.609	1.618	1.623	1.599	1.613
1.584	1.546	1.580	1.506	1.520	1.539	1.529	1.557	1.528	1.522	1.509	1.555	1.600	1.543	1.558	1.598	1.608	1.590	1.621	1.608
1.611	1.597	1.578	1.508	1.527	1.549	1.518	1.523	1.529	1.511	1.531	1.585	1.595	1.574	1.596	1.633	1.643	1.591	1.632	1.586
1.606	1.619	1.554	1.543	1.580	1.582	1.561	1.540	1.569	1.519	1.560	1.614	1.587	1.588	1.597	1.637	1.660	1.626	1.650	1.629
1.604	1.613	1.541	1.507	1.555	1.551	1.517	1.553	1.509	1.502	1.523	1.578	1.568	1.577	1.608	1.639	1.637	1.655	1.629	1.633
1.574	1.610	1.522	1.532	1.548	1.538	1.549	1.552	1.525	1.540	1.546	1.573	1.576	1.620	1.595	1.639	1.625	1.637	1.613	1.630
1.587	1.612	1.560	1.517	1.534	1.538	1.538	1.553	1.549	1.528	1.535	1.549	1.584	1.600	1.604	1.633	1.645	1.623	1.614	1.636
1.592	1.631	1.589	1.569	1.576	1.568	1.569	1.575	1.532	1.535	1.563	1.606	1.567	1.579	1.597	1.658	1.660	1.654	1.633	1.633
1.603	1.612	1.593	1.544	1.548	1.556	1.553	1.560	1.516	1.516	1.553	1.611	1.612	1.570	1.596	1.647	1.609	1.650	1.614	1.630
1.610	1.611	1.579	1.577	1.561	1.544	1.579	1.571	1.574	1.556	1.569	1.636	1.577	1.600	1.581	1.635	1.649	1.633	1.636	1.633
1.621	1.626	1.588	1.570	1.583	1.558	1.580	1.553	1.561	1.534	1.559	1.600	1.603	1.599	1.626	1.621	1.637	1.628	1.627	1.650
1.599	1.646	1.601	1.573	1.607	1.591	1.592	1.571	1.585	1.550	1.597	1.624	1.579	1.596	1.622	1.648	1.618	1.610	1.616	1.611
1.598	1.634	1.569	1.585	1.604	1.576	1.581	1.548	1.576	1.534	1.594	1.624	1.577	1.613	1.640	1.652	1.616	1.635	1.605	1.604
1.595	1.628	1.551	1.568	1.573	1.575	1.604	1.562	1.572	1.540	1.568	1.620	1.592	1.614	1.612	1.640	1.626	1.630	1.617	1.621
1.593	1.622	1.572	1.560	1.566	1.573	1.627	1.570	1.590	1.570	1.610	1.628	1.595	1.627	1.586	1.626	1.599	1.595	1.588	1.626
																			
2.069	2.055	1.967	2.016	2.066	2.049	1.998	1.978	2.039	2.028	1.996	2.074	2.040	2.105	2.116	2.144	2.147	2.147	2.162	2.182
2.207	2.117	2.059	2.051	2.100	2.140	2.151	2.149	2.161	2.126	2.109	2.104	2.137	2.159	2.120	2.121	2.140	2.157	2.184	2.188
2.143	2.138	2.182	2.126	2.107	2.208	2.197	2.194	2.093	2.112	2.127	2.207	2.158	2.141	2.125	2.206	2.196	2.205	2.215	2.243
2.226	2.184	2.242	2.188	2.150	2.194	2.226	2.222	2.130	2.175	2.152	2.291	2.314	2.248	2.227	2.238	2.235	2.253	2.245	2.260
2.195	2.181	2.197	2.198	2.168	2.157	2.226	2.218	2.148	2.162	2.167	2.216	2.252	2.212	2.314	2.318	2.329	2.244	2.326	2.301
2.296	2.316	2.169	2.211	2.186	2.199	2.216	2.191	2.118	2.194	2.163	2.200	2.259	2.278	2.261	2.256	2.326	2.281	2.359	2.364
2.326	2.300	2.192	2.193	2.215	2.210	2.252	2.184	2.146	2.207	2.184	2.179	2.203	2.256	2.270	2.259	2.324	2.295	2.309	2.279
2.249	2.268	2.208	2.221	2.218	2.213	2.215	2.227	2.191	2.183	2.145	2.264	2.254	2.305	2.279	2.273	2.325	2.334	2.286	2.285
2.262	2.217	2.203	2.199	2.221	2.191	2.201	2.250	2.225	2.222	2.161	2.287	2.343	2.329	2.289	2.288	2.395	2.355	2.376	2.348
2.266	2.212	2.236	2.261	2.221	2.173	2.233	2.232	2.179	2.216	2.199	2.216	2.267	2.317	2.281	2.268	2.279	2.343	2.347	2.324
2.242	2.180	2.201	2.251	2.186	2.178	2.236	2.233	2.164	2.209	2.188	2.277	2.323	2.357	2.318	2.314	2.309	2.338	2.364	2.291
2.208	2.247	2.152	2.201	2.201	2.150	2.214	2.201	2.118	2.177	2.130	2.231	2.257	2.322	2.344	2.349	2.343	2.381	2.310	2.348
2.274	2.292	2.206	2.221	2.183	2.142	2.208	2.154	2.199	2.230	2.233	2.230	2.263	2.341	2.319	2.337	2.344	2.333	2.358	2.423
2.224	2.275	2.268	2.223	2.245	2.162	2.193	2.167	2.218	2.223	2.229	2.263	2.291	2.314	2.322	2.289	2.352	2.321	2.308	2.313
2.287	2.228	2.223	2.218	2.163	2.129	2.175	2.161	2.183	2.229	2.201	2.281	2.283	2.360	2.343	2.291	2.348	2.275	2.286	2.369
2.322	2.325	2.224	2.273	2.209	2.207	2.193	2.183	2.133	2.214	2.261	2.247	2.259	2.368	2.355	2.345	2.337	2.259	2.260	2.322
2.299	2.272	2.237	2.252	2.247	2.175	2.236	2.223	2.192	2.227	2.272	2.323	2.288	2.325	2.370	2.321	2.353	2.299	2.266	2.312
2.305	2.302	2.261	2.227	2.227	2.170	2.232	2.173	2.190	2.230	2.273	2.336	2.326	2.341	2.372	2.320	2.334	2.319	2.256	2.338
2.287	2.292	2.236	2.232	2.158	2.227	2.234	2.201	2.221	2.243	2.290	2.387	2.274	2.357	2.321	2.279	2.308	2.210	2.223	2.333
2.305	2.313	2.268	2.261	2.161	2.231	2.240	2.242	2.267	2.277	2.279	2.398	2.294	2.347	2.326	2.238	2.277	2.202	2.212	2.298

fixed vec[rec] rollcv_encf(3)(20,20)
0.339	0.502	0.764	0.919	1.082	1.210	1.355	1.493	1.598	1.703	1.808	2.322	2.704	2.971	3.284	3.585	3.863	4.160	4.280	5.960
0.539	0.760	1.109	1.402	1.632	1.824	1.991	2.200	2.302	2.471	2.652	3.357	3.815	4.290	4.689	5.066	5.503	5.777	6.075	8.306
0.695	0.941	1.409	1.746	2.003	2.275	2.475	2.616	2.750	3.051	3.241	3.991	4.579	5.097	5.699	6.233	6.658	7.096	7.318	10.506
0.791	1.137	1.637	1.974	2.265	2.642	2.872	3.066	3.275	3.516	3.655	4.738	5.358	6.122	6.694	7.160	7.822	8.036	8.534	11.951
0.871	1.249	1.791	2.203	2.464	2.834	3.124	3.435	3.620	3.845	4.095	5.157	6.018	6.865	7.346	7.931	8.490	9.065	9.632	13.389
0.984	1.380	1.950	2.432	2.869	3.169	3.429	3.782	4.011	4.282	4.576	5.738	6.347	7.386	8.125	8.704	9.362	9.748	10.390	14.372
1.066	1.533	2.188	2.691	3.138	3.461	3.800	4.086	4.455	4.741	4.951	6.320	7.151	7.847	8.654	9.417	10.013	10.532	11.037	15.746
1.148	1.616	2.316	2.900	3.377	3.789	4.124	4.478	4.836	4.958	5.396	6.722	7.624	8.417	9.356	10.085	10.789	11.225	11.953	16.684
1.226	1.779	2.471	3.096	3.605	4.124	4.422	4.773	5.149	5.449	5.693	7.004	8.115	8.785	9.990	10.735	11.542	12.267	12.935	17.766
1.284	1.838	2.627	3.149	3.763	4.186	4.580	4.937	5.273	5.538	5.968	7.227	8.463	9.478	10.349	11.595	12.251	12.894	13.366	18.645
1.343	1.955	2.777	3.371	3.948	4.444	4.757	5.244	5.479	5.745	6.244	7.530	8.794	9.861	10.887	12.000	12.707	13.214	13.942	19.642
1.415	2.100	2.945	3.535	4.091	4.592	5.042	5.555	5.841	6.115	6.529	7.797	9.114	10.320	11.296	12.416	13.110	13.791	14.554	20.152
1.489	2.156	3.063	3.684	4.296	4.821	5.292	5.761	6.096	6.300	6.731	8.315	9.555	10.758	11.760	12.918	13.531	14.400	15.140	20.999
1.532	2.224	3.090	3.785	4.445	4.986	5.414	5.820	6.231	6.401	6.921	8.521	9.864	10.941	12.215	13.198	14.215	14.934	15.559	22.123
1.568	2.278	3.248	4.008	4.584	5.182	5.682	6.052	6.354	6.658	7.185	8.855	10.091	11.321	12.551	13.713	14.567	15.578	16.373	22.773
1.628	2.375	3.361	4.177	4.793	5.409	5.920	6.092	6.547	6.929	7.296	9.212	10.437	11.707	12.961	14.198	15.016	15.937	16.840	23.399
1.674	2.431	3.510	4.169	4.978	5.452	6.103	6.449	6.846	7.202	7.684	9.432	10.786	12.256	13.295	14.838	15.303	16.328	17.136	23.841
1.697	2.486	3.578	4.280	5.032	5.511	6.226	6.680	6.971	7.435	7.859	9.745	11.256	12.830	14.033	15.252	15.928	16.874	17.684	24.391
1.776	2.573	3.683	4.433	5.181	5.641	6.388	6.795	7.316	7.774	8.158	9.914	11.468	13.233	14.150	15.622	16.364	16.804	18.099	25.190
1.824	2.655	3.790	4.566	5.422	5.896	6.573	6.957	7.575	7.858	8.520	10.252	11.761	13.754	14.564	15.973	16.686	17.299	18.271	25.568
																			
0.534	0.788	1.161	1.506	1.722	1.946	2.198	2.374	2.498	2.683	2.836	3.444	3.926	4.404	4.792	5.202	5.472	5.851	6.112	8.055
0.793	1.089	1.640	2.031	2.380	2.704	2.949	3.212	3.366	3.593	3.723	4.678	5.327	5.881	6.259	6.899	7.440	7.772	8.133	11.094
0.962	1.391	1.955	2.500	2.847	3.271	3.549	3.759	4.108	4.219	4.435	5.630	6.426	7.061	7.759	8.276	8.955	9.478	10.013	13.746
1.087	1.607	2.310	2.779	3.220	3.681	4.079	4.510	4.811	5.051	5.289	6.654	7.505	8.346	9.030	9.823	10.631	11.010	11.419	16.201
1.211	1.773	2.545	3.073	3.522	4.082	4.468	4.899	5.199	5.470	5.840	7.116	8.222	9.011	9.774	10.593	11.682	11.975	12.564	17.535
1.334	1.870	2.766	3.413	3.939	4.473	4.881	5.326	5.675	5.979	6.288	7.917	8.886	9.828	10.642	11.689	12.526	13.073	13.575	18.880
1.465	2.082	3.029	3.696	4.299	4.904	5.364	5.805	6.199	6.389	6.763	8.346	9.694	10.324	11.207	12.478	13.318	13.904	14.768	20.317
1.557	2.241	3.217	3.865	4.610	5.306	5.644	6.049	6.484	6.759	7.231	8.863	10.287	11.088	12.236	13.576	14.497	14.764	15.905	21.456
1.657	2.403	3.397	4.223	4.972	5.585	5.996	6.432	6.860	7.111	7.670	9.474	10.836	11.724	13.078	14.396	15.568	16.214	17.062	23.168
1.709	2.526	3.519	4.382	5.138	5.719	6.129	6.694	6.975	7.349	7.867	9.777	11.268	12.392	13.679	15.124	16.006	16.942	17.752	24.402
1.772	2.611	3.650	4.554	5.359	5.919	6.600	7.069	7.542	7.902	8.340	10.258	11.644	13.226	14.153	15.717	16.553	17.599	18.149	25.433
1.854	2.695	3.838	4.694	5.478	6.233	6.826	7.297	7.867	8.142	8.663	10.519	12.145	13.697	14.828	16.278	17.510	18.391	19.190	27.001
1.957	2.828	4.022	5.007	5.832	6.558	7.150	7.595	7.977	8.572	9.143	11.252	12.581	13.866	15.386	17.151	18.349	19.170	19.968	28.000
2.021	2.853	4.148	5.055	5.904	6.621	7.433	7.939	8.236	8.792	9.281	11.714	13.408	14.525	15.794	17.494	18.347	19.813	20.274	28.719
2.109	2.986	4.275	5.303	6.185	6.826	7.610	8.102	8.661	9.119	9.751	12.096	13.607	14.973	16.161	18.030	19.484	20.293	21.311	29.578
2.193	3.184	4.468	5.549	6.451	7.054	7.861	8.314	9.027	9.220	10.007	12.201	14.078	15.650	17.147	18.547	19.951	20.833	21.863	30.678
2.237	3.324	4.618	5.748	6.641	7.406	8.028	8.643	9.389	9.539	10.377	12.852	14.329	16.054	17.626	19.179	20.216	21.288	22.566	30.872
2.282	3.341	4.695	5.869	6.799	7.582	8.357	8.729	9.455	9.833	10.727	13.165	14.711	16.664	18.305	19.892	20.659	21.927	22.891	31.823
2.329	3.450	4.733	5.925	6.905	7.701	8.615	9.078	9.713	10.069	10.791	13.479	15.261	17.150	18.550	20.413	21.270	22.664	23.693	32.707
2.386	3.520	4.906	6.044	7.031	7.889	8.961	9.331	10.095	10.547	11.160	13.926	15.645	17.642	18.687	20.458	21.476	22.661	23.871	33.728
																			
0.997	1.466	2.278	2.854	3.264	3.676	3.852	4.410	4.614	4.835	5.064	5.976	6.673	7.169	7.625	8.278	8.743	9.165	9.382	12.436
1.282	1.999	2.929	3.596	4.282	4.799	5.389	5.568	6.062	6.338	6.502	7.713	8.445	9.371	9.996	10.556	11.146	11.657	12.392	16.575
1.561	2.237	3.401	4.218	4.893	5.609	6.055	6.442	6.646	7.198	7.513	9.112	9.943	10.894	11.529	12.840	13.423	14.235	14.563	19.995
1.783	2.588	3.681	4.603	5.312	6.192	6.845	7.326	7.514	8.173	8.529	10.205	11.835	12.709	13.356	14.748	15.608	16.272	16.988	23.073
1.921	2.830	4.039	5.028	5.918	6.775	7.423	7.848	8.310	8.578	9.163	10.884	12.495	13.924	15.160	16.584	17.399	17.741	19.061	26.182
2.120	2.945	4.424	5.303	6.480	7.193	8.140	8.407	8.768	9.404	9.885	11.636	13.851	15.141	16.172	17.767	19.090	19.198	20.819	29.041
2.324	3.213	4.580	5.799	6.670	7.609	8.313	8.759	9.160	10.097	10.305	12.616	14.027	16.018	16.996	18.585	20.212	20.732	22.140	29.920
2.376	3.461	4.945	6.097	7.165	7.964	8.792	9.237	9.982	10.789	10.985	13.769	15.540	17.160	18.362	19.807	21.587	22.318	22.977	31.732
2.538	3.531	5.151	6.382	7.546	8.277	9.132	10.226	10.724	11.154	11.710	14.426	16.771	18.596	19.701	20.940	23.305	23.823	25.638	34.297
2.569	3.600	5.365	6.699	7.705	8.624	9.630	10.375	10.928	11.573	12.043	14.784	16.860	19.439	20.924	21.356	22.919	24.524	26.531	36.029
2.663	3.754	5.572	6.884	8.089	9.170	10.221	11.022	11.351	12.349	12.537	15.764	17.809	19.900	21.868	22.906	24.435	26.439	27.427	36.947
2.786	4.006	5.708	7.052	8.346	9.340	10.429	11.210	11.423	12.648	12.663	15.794	18.115	20.854	22.553	24.255	25.555	27.727	28.112	38.917
2.879	4.327	5.856	7.429	8.464	9.635	10.753	11.313	12.216	13.128	13.597	16.670	18.926	21.456	23.009	25.162	26.692	27.646	29.564	41.710
2.926	4.449	6.257	7.697	9.079	10.068	10.898	11.699	12.758	13.447	14.211	17.050	19.530	21.930	24.135	25.370	27.658	28.870	29.806	41.324
2.970	4.519	6.232	7.944	9.080	10.314	11.207	11.843	12.898	13.737	14.526	18.101	20.347	23.032	25.113	26.501	28.263	29.245	30.937	43.848
3.271	4.626	6.588	8.074	9.261	10.562	11.802	12.623	12.996	14.237	15.335	18.105	20.568	24.236	26.390	27.513	29.049	30.013	31.032	44.676
3.288	4.760	6.813	8.458	9.920	10.985	12.224	12.875	13.512	14.460	15.740	19.317	21.576	24.559	26.556	27.911	29.993	30.653	31.881	45.911
3.421	4.833	7.074	8.411	9.809	11.118	12.365	13.331	14.142	15.094	16.396	19.819	22.334	24.566	27.303	28.820	30.451	32.186	32.686	47.530
3.373	4.924	7.273	8.697	9.981	11.521	12.680	13.475	14.354	15.690	16.709	20.294	22.800	25.658	27.639	28.728	30.755	31.440	33.539	48.531
3.556	5.041	7.458	9.082	10.233	11.903	13.208	13.906	14.924	15.971	16.632	21.010	23.361	26.603	28.238	29.016	31.073	32.646	33.646	48.706

fixed vec[rec] fixcv_mset(3)(20,20)
1.149	1.079	1.042	1.040	0.976	0.917	0.896	0.893	0.908	0.834	0.862	0.755	0.689	0.637	0.545	0.553	0.525	0.485	0.450	0.306
1.096	0.999	0.901	0.873	0.798	0.711	0.680	0.639	0.578	0.556	0.520	0.375	0.270	0.150	0.046	0.015	-0.050	-0.108	-0.225	-0.617
1.063	0.922	0.793	0.705	0.621	0.540	0.511	0.455	0.386	0.373	0.306	0.089	-0.057	-0.216	-0.300	-0.397	-0.533	-0.625	-0.715	-1.314
1.004	0.895	0.764	0.575	0.476	0.367	0.340	0.273	0.204	0.171	0.081	-0.141	-0.351	-0.525	-0.651	-0.740	-0.922	-1.056	-1.155	-1.909
0.945	0.838	0.636	0.487	0.374	0.258	0.193	0.115	0.020	-0.022	-0.085	-0.417	-0.604	-0.808	-0.969	-1.070	-1.236	-1.393	-1.526	-2.417
0.912	0.764	0.552	0.400	0.299	0.169	0.103	0.003	-0.106	-0.146	-0.235	-0.599	-0.819	-1.024	-1.201	-1.350	-1.542	-1.743	-1.842	-2.871
0.895	0.731	0.513	0.332	0.215	0.060	-0.003	-0.147	-0.252	-0.308	-0.386	-0.756	-1.021	-1.266	-1.440	-1.594	-1.782	-1.992	-2.113	-3.343
0.868	0.663	0.467	0.247	0.153	-0.029	-0.115	-0.227	-0.343	-0.440	-0.502	-0.924	-1.217	-1.470	-1.714	-1.866	-2.042	-2.269	-2.445	-3.754
0.844	0.655	0.396	0.182	0.034	-0.111	-0.224	-0.303	-0.437	-0.543	-0.625	-1.062	-1.359	-1.676	-1.903	-2.077	-2.309	-2.526	-2.689	-4.119
0.797	0.616	0.348	0.125	-0.055	-0.210	-0.305	-0.398	-0.559	-0.645	-0.729	-1.248	-1.568	-1.871	-2.120	-2.322	-2.575	-2.787	-2.979	-4.548
0.794	0.579	0.287	0.055	-0.151	-0.299	-0.401	-0.511	-0.665	-0.745	-0.847	-1.385	-1.724	-2.037	-2.316	-2.516	-2.810	-3.058	-3.277	-4.884
0.761	0.562	0.254	-0.001	-0.196	-0.362	-0.477	-0.597	-0.736	-0.849	-0.967	-1.512	-1.883	-2.228	-2.533	-2.775	-2.988	-3.284	-3.522	-5.159
0.733	0.531	0.225	-0.065	-0.229	-0.423	-0.557	-0.659	-0.834	-0.977	-1.077	-1.610	-2.004	-2.359	-2.686	-2.939	-3.184	-3.485	-3.747	-5.551
0.719	0.494	0.181	-0.123	-0.299	-0.486	-0.630	-0.767	-0.954	-1.050	-1.200	-1.731	-2.131	-2.490	-2.840	-3.128	-3.388	-3.697	-3.947	-5.841
0.698	0.457	0.134	-0.172	-0.350	-0.536	-0.717	-0.860	-1.009	-1.150	-1.293	-1.838	-2.278	-2.661	-2.993	-3.300	-3.597	-3.879	-4.132	-6.089
0.670	0.449	0.106	-0.202	-0.391	-0.620	-0.760	-0.920	-1.081	-1.220	-1.371	-1.961	-2.388	-2.750	-3.127	-3.430	-3.752	-4.034	-4.344	-6.434
0.648	0.411	0.073	-0.238	-0.454	-0.657	-0.838	-1.017	-1.175	-1.299	-1.447	-2.045	-2.534	-2.922	-3.310	-3.641	-3.947	-4.234	-4.542	-6.692
0.637	0.387	0.020	-0.316	-0.509	-0.720	-0.921	-1.103	-1.263	-1.403	-1.530	-2.146	-2.644	-3.033	-3.453	-3.790	-4.135	-4.403	-4.720	-6.979
0.617	0.340	-0.026	-0.367	-0.583	-0.800	-0.985	-1.178	-1.348	-1.505	-1.623	-2.236	-2.754	-3.138	-3.609	-3.955	-4.303	-4.633	-4.903	-7.212
0.578	0.309	-0.055	-0.409	-0.632	-0.826	-1.029	-1.248	-1.420	-1.566	-1.723	-2.335	-2.858	-3.301	-3.763	-4.088	-4.459	-4.800	-5.110	-7.483
																			
1.506	1.416	1.364	1.428	1.346	1.252	1.301	1.293	1.249	1.235	1.218	1.138	1.120	1.041	0.977	0.993	0.965	0.956	0.903	0.745
1.468	1.342	1.301	1.265	1.164	1.072	1.034	1.046	0.977	0.982	0.955	0.800	0.661	0.548	0.458	0.397	0.386	0.362	0.340	-0.154
1.432	1.277	1.195	1.095	1.014	0.909	0.893	0.851	0.761	0.735	0.733	0.504	0.353	0.209	0.105	0.035	-0.093	-0.222	-0.275	-0.823
1.369	1.281	1.110	0.997	0.883	0.755	0.689	0.650	0.607	0.566	0.509	0.257	0.052	-0.073	-0.201	-0.293	-0.470	-0.551	-0.658	-1.300
1.333	1.193	1.009	0.863	0.725	0.646	0.570	0.486	0.410	0.365	0.291	0.029	-0.190	-0.334	-0.498	-0.582	-0.776	-0.918	-1.020	-1.785
1.269	1.122	0.936	0.771	0.652	0.538	0.487	0.367	0.314	0.222	0.152	-0.145	-0.383	-0.602	-0.716	-0.862	-1.054	-1.211	-1.346	-2.350
1.263	1.086	0.878	0.692	0.557	0.446	0.346	0.254	0.191	0.074	0.014	-0.335	-0.576	-0.822	-0.953	-1.076	-1.300	-1.473	-1.658	-2.758
1.249	1.064	0.807	0.623	0.481	0.363	0.268	0.151	0.054	-0.042	-0.120	-0.524	-0.771	-1.053	-1.185	-1.362	-1.548	-1.783	-1.919	-3.131
1.197	1.031	0.754	0.537	0.416	0.305	0.162	0.050	-0.067	-0.171	-0.242	-0.655	-0.964	-1.252	-1.417	-1.604	-1.804	-2.027	-2.159	-3.524
1.143	1.007	0.688	0.455	0.337	0.167	0.040	-0.057	-0.174	-0.246	-0.358	-0.786	-1.151	-1.453	-1.645	-1.820	-2.040	-2.307	-2.466	-3.939
1.135	0.948	0.677	0.395	0.259	0.102	-0.045	-0.136	-0.270	-0.354	-0.471	-0.970	-1.307	-1.646	-1.841	-2.072	-2.337	-2.543	-2.777	-4.283
1.118	0.912	0.618	0.321	0.206	0.046	-0.090	-0.177	-0.377	-0.439	-0.563	-1.144	-1.423	-1.791	-2.065	-2.224	-2.520	-2.777	-3.012	-4.616
1.090	0.876	0.603	0.298	0.128	-0.045	-0.158	-0.279	-0.469	-0.530	-0.643	-1.220	-1.566	-1.939	-2.173	-2.393	-2.707	-2.981	-3.244	-4.880
1.079	0.845	0.537	0.226	0.051	-0.111	-0.258	-0.403	-0.569	-0.674	-0.771	-1.294	-1.743	-2.071	-2.329	-2.623	-2.924	-3.176	-3.447	-5.203
1.074	0.811	0.479	0.187	-0.015	-0.192	-0.322	-0.463	-0.680	-0.764	-0.866	-1.405	-1.879	-2.225	-2.481	-2.840	-3.092	-3.404	-3.688	-5.471
1.054	0.812	0.456	0.181	-0.033	-0.233	-0.365	-0.528	-0.759	-0.856	-0.948	-1.511	-1.950	-2.364	-2.629	-2.927	-3.243	-3.535	-3.846	-5.780
1.041	0.792	0.428	0.116	-0.080	-0.287	-0.428	-0.605	-0.795	-0.953	-1.047	-1.645	-2.112	-2.488	-2.801	-3.159	-3.451	-3.748	-4.009	-6.055
0.965	0.748	0.402	0.063	-0.152	-0.362	-0.513	-0.673	-0.872	-1.014	-1.114	-1.728	-2.215	-2.606	-2.900	-3.292	-3.603	-3.912	-4.193	-6.322
0.963	0.722	0.361	-0.001	-0.215	-0.420	-0.552	-0.769	-0.955	-1.108	-1.210	-1.815	-2.340	-2.724	-3.128	-3.457	-3.773	-4.099	-4.403	-6.558
0.943	0.699	0.308	-0.059	-0.272	-0.459	-0.635	-0.846	-1.051	-1.178	-1.291	-1.936	-2.446	-2.829	-3.226	-3.628	-3.968	-4.279	-4.590	-6.815
																			
2.201	2.051	1.974	2.061	2.037	2.024	1.992	2.018	1.996	2.016	1.993	1.887	1.826	1.738	1.741	1.742	1.696	1.668	1.679	1.555
2.145	2.089	1.923	1.947	1.964	1.749	1.751	1.665	1.725	1.646	1.613	1.493	1.534	1.305	1.160	1.346	1.246	1.133	1.095	0.782
2.045	1.977	1.957	1.805	1.739	1.602	1.520	1.597	1.463	1.513	1.407	1.201	1.139	1.003	0.903	0.892	0.769	0.627	0.574	0.098
2.013	1.883	1.829	1.687	1.528	1.467	1.475	1.422	1.318	1.255	1.277	1.028	0.859	0.704	0.625	0.489	0.345	0.253	0.307	-0.330
1.930	1.878	1.716	1.596	1.405	1.254	1.301	1.230	1.171	1.115	1.034	0.750	0.566	0.429	0.426	0.297	0.038	-0.058	-0.009	-0.832
1.933	1.874	1.628	1.481	1.382	1.146	1.188	1.091	1.016	1.007	0.878	0.544	0.356	0.284	0.150	-0.033	-0.253	-0.334	-0.370	-1.338
1.925	1.859	1.556	1.377	1.257	1.105	1.103	0.987	0.896	0.828	0.765	0.409	0.169	0.072	-0.104	-0.278	-0.448	-0.590	-0.724	-1.663
1.856	1.827	1.467	1.245	1.146	1.029	0.980	0.860	0.786	0.762	0.666	0.230	0.051	-0.185	-0.366	-0.549	-0.669	-0.845	-0.972	-2.020
1.878	1.697	1.440	1.198	1.124	0.902	0.791	0.683	0.644	0.595	0.507	0.082	-0.160	-0.378	-0.568	-0.739	-0.972	-1.034	-1.158	-2.400
1.824	1.604	1.354	1.126	0.998	0.797	0.659	0.557	0.550	0.505	0.415	-0.065	-0.389	-0.581	-0.837	-1.047	-1.283	-1.321	-1.577	-2.875
1.800	1.580	1.359	1.095	0.925	0.800	0.623	0.488	0.470	0.447	0.305	-0.188	-0.517	-0.735	-1.031	-1.208	-1.410	-1.591	-1.752	-3.224
1.758	1.564	1.218	1.073	0.807	0.723	0.560	0.398	0.351	0.333	0.263	-0.282	-0.667	-0.915	-1.174	-1.410	-1.612	-1.868	-2.029	-3.556
1.764	1.542	1.141	1.017	0.776	0.639	0.507	0.360	0.275	0.253	0.120	-0.407	-0.808	-1.077	-1.323	-1.558	-1.857	-2.129	-2.272	-3.769
1.679	1.479	1.168	0.992	0.709	0.566	0.350	0.300	0.176	0.152	0.033	-0.591	-0.953	-1.266	-1.511	-1.752	-2.004	-2.302	-2.479	-4.131
1.715	1.462	1.108	0.890	0.638	0.459	0.316	0.222	0.072	0.011	-0.077	-0.741	-1.139	-1.399	-1.635	-1.894	-2.259	-2.472	-2.707	-4.348
1.694	1.451	1.097	0.872	0.622	0.441	0.291	0.189	0.005	-0.085	-0.118	-0.815	-1.274	-1.518	-1.716	-2.082	-2.413	-2.634	-2.788	-4.713
1.672	1.432	1.041	0.837	0.588	0.422	0.248	0.119	-0.083	-0.149	-0.214	-0.886	-1.368	-1.682	-1.935	-2.232	-2.586	-2.817	-2.971	-5.023
1.674	1.414	1.019	0.790	0.534	0.382	0.140	0.058	-0.129	-0.247	-0.322	-0.980	-1.531	-1.844	-1.990	-2.387	-2.690	-2.936	-3.164	-5.208
1.668	1.381	1.018	0.763	0.479	0.299	0.052	0.001	-0.209	-0.298	-0.434	-1.134	-1.599	-1.996	-2.151	-2.589	-2.850	-3.108	-3.373	-5.442
1.663	1.391	1.002	0.742	0.408	0.185	0.016	-0.045	-0.299	-0.389	-0.539	-1.293	-1.702	-2.087	-2.327	-2.723	-3.045	-3.317	-3.621	-5.659

fixed vec[rec] fixcv_msef(3)(20,20)
0.514	0.649	0.835	0.885	0.933	0.964	1.009	0.986	1.050	1.080	1.037	0.957	0.923	0.861	0.722	0.765	0.693	0.587	0.550	0.308
0.784	0.914	1.140	1.237	1.299	1.268	1.330	1.250	1.126	1.189	1.151	0.800	0.691	0.376	0.120	0.031	-0.116	-0.296	-0.687	-2.605
1.000	1.106	1.328	1.367	1.375	1.291	1.289	1.225	1.073	1.043	0.954	0.265	-0.200	-0.796	-1.248	-1.827	-2.551	-3.169	-3.792	-9.279
1.100	1.317	1.472	1.362	1.195	1.109	1.008	0.860	0.667	0.544	0.289	-0.604	-1.696	-2.730	-3.701	-4.521	-5.681	-7.041	-8.082	-18.289
1.178	1.387	1.414	1.341	1.038	0.865	0.696	0.445	0.083	-0.075	-0.347	-2.069	-3.561	-5.026	-6.806	-7.857	-9.601	-11.455	-13.472	-28.749
1.242	1.417	1.428	1.167	0.962	0.634	0.410	0.014	-0.449	-0.666	-1.113	-3.418	-5.657	-7.586	-9.613	-11.722	-13.882	-17.019	-18.341	-39.491
1.342	1.536	1.457	1.057	0.817	0.254	-0.015	-0.668	-1.314	-1.610	-2.230	-5.111	-7.818	-10.307	-12.774	-15.063	-18.230	-21.317	-24.030	-51.222
1.434	1.501	1.374	0.900	0.622	-0.146	-0.593	-1.157	-1.925	-2.678	-3.109	-6.961	-9.965	-13.680	-16.896	-19.936	-23.395	-26.936	-30.062	-64.846
1.474	1.564	1.329	0.727	0.181	-0.578	-1.168	-1.751	-2.653	-3.457	-4.122	-8.539	-12.122	-16.987	-20.567	-23.677	-27.674	-32.490	-36.047	-76.473
1.504	1.583	1.176	0.520	-0.246	-1.125	-1.722	-2.449	-3.606	-4.312	-5.252	-10.586	-15.099	-20.027	-24.210	-28.406	-33.185	-38.608	-43.478	-91.402
1.566	1.541	1.080	0.251	-0.753	-1.611	-2.439	-3.266	-4.483	-5.422	-6.360	-12.672	-17.996	-23.134	-28.287	-32.924	-39.199	-45.102	-50.862	-106.128
1.566	1.562	0.968	-0.006	-1.030	-2.150	-2.877	-3.960	-5.289	-6.447	-7.715	-14.328	-20.487	-27.196	-33.381	-38.311	-45.213	-51.426	-58.336	-119.656
1.580	1.545	0.957	-0.299	-1.365	-2.626	-3.689	-4.726	-6.272	-7.509	-9.007	-16.348	-23.257	-30.388	-36.584	-43.373	-50.603	-57.290	-65.122	-132.123
1.594	1.498	0.792	-0.648	-1.770	-3.148	-4.440	-5.718	-7.488	-8.646	-10.395	-17.999	-25.517	-34.032	-41.050	-48.327	-56.580	-64.914	-72.925	-148.285
1.622	1.489	0.588	-0.940	-2.081	-3.645	-5.146	-6.612	-8.600	-9.970	-11.804	-20.106	-28.919	-37.653	-45.706	-52.938	-62.103	-71.476	-80.935	-164.664
1.620	1.443	0.464	-1.097	-2.548	-4.343	-5.803	-7.524	-9.490	-11.263	-13.214	-22.300	-31.717	-40.832	-49.181	-58.075	-67.410	-77.936	-87.885	-179.386
1.628	1.416	0.339	-1.390	-2.997	-4.831	-6.562	-8.453	-10.475	-12.263	-14.445	-24.392	-34.979	-44.682	-55.186	-63.559	-74.458	-85.295	-95.432	-195.274
1.596	1.339	0.094	-1.882	-3.479	-5.399	-7.493	-9.469	-11.611	-13.466	-15.798	-26.499	-37.837	-47.966	-59.345	-68.965	-80.168	-92.782	-102.542 -210.560
1.600	1.248	-0.137	-2.285	-4.003	-6.152	-8.135	-10.306	-12.680	-15.123	-17.111	-28.189	-40.515	-51.424	-64.158	-74.913	-86.851	-98.571	-110.813 -224.065
1.558	1.143	-0.302	-2.549	-4.490	-6.736	-8.928	-11.443	-13.952	-16.618	-18.724	-30.654	-43.074	-55.741	-68.981	-79.489	-92.305	-105.992 -117.984 -238.681
																			
0.784	1.015	1.345	1.534	1.677	1.667	1.738	1.807	1.812	1.857	1.862	1.770	1.841	1.753	1.649	1.632	1.554	1.544	1.478	1.262
1.132	1.421	1.765	1.999	2.077	2.116	2.169	2.232	2.275	2.260	2.195	2.073	1.759	1.496	1.242	1.132	1.085	0.959	0.894	-0.535
1.408	1.653	2.050	2.322	2.308	2.319	2.325	2.336	2.187	2.283	2.275	1.814	1.054	0.849	0.392	0.164	-0.417	-0.934	-1.163	-4.594
1.536	1.947	2.374	2.478	2.426	2.238	2.310	2.175	2.063	1.891	1.784	1.047	0.262	-0.367	-0.916	-1.557	-2.594	-3.330	-3.927	-10.835
1.664	2.018	2.368	2.504	2.337	2.167	2.109	1.862	1.593	1.540	1.249	0.133	-1.113	-1.983	-3.288	-3.951	-5.222	-6.432	-7.664	-18.743
1.826	2.164	2.406	2.422	2.267	2.010	1.995	1.654	1.302	1.107	0.744	-0.829	-2.314	-4.007	-5.415	-6.667	-8.843	-10.700	-12.226	-27.849
1.962	2.282	2.441	2.410	2.198	1.886	1.535	1.263	0.943	0.356	0.071	-2.188	-3.928	-6.471	-8.003	-9.508	-12.391	-15.151	-16.979	-38.138
2.078	2.394	2.580	2.244	2.007	1.666	1.266	0.744	0.293	-0.244	-0.658	-3.442	-6.198	-8.998	-10.438	-13.121	-16.672	-19.548	-22.064	-48.672
2.090	2.525	2.533	2.064	1.881	1.434	0.892	0.292	-0.344	-0.960	-1.536	-4.878	-8.212	-11.668	-14.490	-17.140	-20.966	-24.581	-27.326	-60.424
2.158	2.611	2.481	1.967	1.601	0.936	0.282	-0.360	-1.068	-1.467	-2.404	-6.509	-10.288	-15.156	-17.776	-20.844	-25.136	-30.016	-33.123	-74.059
2.202	2.594	2.526	1.801	1.219	0.561	-0.306	-0.928	-1.872	-2.308	-3.393	-8.034	-12.471	-17.537	-21.213	-25.084	-29.762	-35.993	-39.885	-86.071
2.298	2.622	2.418	1.516	0.944	0.272	-0.559	-1.218	-2.666	-3.140	-4.102	-10.225	-14.956	-20.680	-25.598	-29.134	-34.751	-42.005	-46.922	-97.657
2.336	2.567	2.433	1.449	0.714	-0.266	-1.113	-1.980	-3.411	-4.017	-5.302	-11.821	-17.560	-23.417	-28.646	-33.110	-39.468	-47.133	-53.385	-108.740
2.452	2.624	2.294	1.215	0.293	-0.739	-1.783	-2.906	-4.590	-5.499	-6.323	-13.173	-20.216	-26.382	-32.391	-37.973	-45.684	-52.946	-60.102	-122.480
2.574	2.610	2.151	0.973	-0.079	-1.230	-2.266	-3.545	-5.644	-6.460	-7.474	-15.480	-23.010	-29.530	-35.791	-43.345	-50.717	-58.338	-67.068	-135.357
2.572	2.721	1.991	1.013	-0.225	-1.575	-2.627	-4.245	-6.298	-7.584	-8.606	-17.078	-25.278	-32.761	-38.869	-47.723	-55.546	-63.665	-72.742	-150.310
2.566	2.711	1.997	0.617	-0.563	-1.978	-3.315	-4.958	-6.971	-8.601	-9.984	-18.903	-28.154	-36.373	-44.278	-53.283	-61.521	-70.683	-79.521	-165.811
2.470	2.570	1.920	0.389	-1.016	-2.656	-4.025	-5.907	-7.889	-9.542	-11.311	-20.608	-30.694	-39.697	-48.679	-57.323	-66.864	-76.922	-87.265	-178.445
2.536	2.630	1.830	-0.004	-1.488	-3.195	-4.450	-6.706	-8.920	-10.610	-12.449	-22.753	-33.771	-42.016	-52.872	-62.082	-72.466	-83.876	-94.647	-190.485
2.556	2.661	1.555	-0.382	-1.911	-3.606	-5.306	-7.623	-9.935	-11.587	-13.583	-24.310	-35.778	-45.202	-56.598	-66.612	-78.954	-90.799	-101.847 -203.429
																			
1.480	1.981	2.681	3.055	3.230	3.377	3.562	3.619	3.816	3.812	3.838	3.960	3.927	4.337	4.244	4.107	4.247	4.018	4.012	3.658
1.840	2.554	3.241	3.514	3.944	4.019	4.173	4.364	4.251	4.556	4.414	4.622	4.466	4.207	3.995	4.251	4.224	3.910	3.991	2.561
2.324	2.985	3.854	4.103	4.272	4.233	4.549	4.764	4.687	4.915	4.900	4.631	4.313	3.896	3.850	4.026	3.183	2.617	2.255	0.538
2.576	3.283	3.999	4.339	4.349	4.629	4.602	5.002	4.793	4.984	5.028	4.480	3.668	3.179	2.979	2.444	1.558	1.303	1.596	-2.003
2.748	3.437	4.212	4.454	4.330	4.396	4.739	5.044	4.761	4.731	4.560	3.803	2.921	2.624	2.344	1.717	0.194	-0.345	-0.068	-6.705
3.084	3.755	4.467	4.754	4.559	4.715	4.836	4.515	4.561	4.303	4.365	3.329	2.227	1.853	1.023	-0.236	-1.779	-2.592	-2.720	-13.120
3.294	3.980	4.599	4.683	4.704	5.000	4.828	4.667	4.489	4.367	4.155	2.779	1.091	0.355	-0.822	-2.005	-3.751	-5.411	-6.243	-18.836
3.364	4.116	4.775	4.724	4.715	4.762	4.480	4.111	4.278	4.482	3.804	1.625	0.378	-1.397	-2.971	-4.273	-6.005	-7.900	-8.910	-26.277
3.430	4.233	4.671	4.640	4.856	4.580	4.112	3.756	3.536	3.648	3.158	0.533	-1.217	-3.206	-4.992	-7.344	-8.837	-10.757	-12.406	-33.987
3.582	4.232	4.750	4.674	4.489	4.251	3.643	3.467	3.373	3.342	3.036	-0.638	-3.186	-5.179	-7.871	-10.533	-13.752	-14.795	-17.141	-46.000
3.742	4.374	4.880	4.680	4.769	4.268	3.438	3.053	3.204	2.829	2.098	-1.563	-4.253	-7.054	-9.941	-13.451	-16.649	-19.556	-21.690	-54.847
3.862	4.636	5.026	4.755	4.232	4.346	3.375	2.557	2.281	2.362	1.791	-2.623	-6.219	-9.787	-12.776	-16.573	-20.772	-24.261	-26.881	-64.146
3.954	4.783	4.850	4.976	4.587	4.032	3.160	2.312	1.885	1.722	0.753	-3.846	-8.272	-12.142	-15.430	-19.608	-24.686	-29.382	-32.890	-70.783
4.014	4.758	5.121	4.906	3.935	3.549	2.403	2.197	1.390	1.065	0.287	-5.782	-10.023	-15.302	-18.497	-23.827	-28.690	-34.315	-38.540	-84.601
4.042	4.783	4.873	4.622	3.953	3.045	2.073	1.584	0.475	0.086	-0.490	-7.175	-12.391	-17.664	-21.937	-27.239	-33.775	-38.709	-42.448	-92.432
4.334	5.111	5.051	4.679	3.792	2.978	1.920	1.206	0.056	-0.533	-1.106	-7.655	-14.229	-19.497	-23.680	-29.701	-37.624	-43.033	-46.541	-105.645
4.300	5.205	5.024	4.722	3.638	3.051	1.742	0.974	-0.644	-1.317	-2.018	-9.477	-16.788	-22.043	-27.055	-34.301	-42.256	-48.005	-53.700	-124.039
4.384	5.229	5.274	4.627	3.715	2.601	1.043	0.473	-0.997	-2.321	-2.862	-11.457	-18.718	-25.474	-29.504	-37.616	-44.904	-50.125	-60.827	-131.029
4.392	5.112	5.269	4.588	3.179	2.087	0.367	0.014	-2.177	-2.796	-3.903	-13.976	-21.119	-29.109	-33.151	-42.303	-49.900	-57.751	-66.032	-140.844
4.488	5.308	5.117	4.478	2.858	1.596	0.111	-0.390	-2.622	-3.966	-5.513	-15.699	-23.350	-31.229	-35.956	-46.756	-57.589	-64.564	-73.014	-151.974

fixed vec[rec] fixcv_enct(3)(20,20)
1.269	1.233	1.272	1.310	1.299	1.256	1.295	1.304	1.297	1.305	1.301	1.282	1.292	1.260	1.255	1.286	1.261	1.269	1.266	1.278
1.293	1.254	1.269	1.316	1.304	1.284	1.312	1.287	1.272	1.286	1.282	1.268	1.265	1.247	1.248	1.249	1.247	1.242	1.247	1.239
1.301	1.250	1.269	1.287	1.306	1.305	1.339	1.314	1.290	1.291	1.311	1.265	1.237	1.248	1.274	1.249	1.243	1.235	1.231	1.277
1.292	1.278	1.302	1.272	1.261	1.284	1.286	1.267	1.292	1.271	1.273	1.278	1.264	1.267	1.287	1.272	1.253	1.228	1.216	1.276
1.274	1.278	1.262	1.261	1.232	1.248	1.284	1.263	1.239	1.254	1.258	1.261	1.256	1.281	1.298	1.259	1.269	1.257	1.260	1.249
1.269	1.255	1.271	1.262	1.261	1.275	1.316	1.290	1.281	1.284	1.290	1.252	1.273	1.284	1.328	1.290	1.279	1.261	1.260	1.259
1.298	1.285	1.305	1.275	1.286	1.284	1.305	1.262	1.261	1.265	1.264	1.249	1.268	1.262	1.311	1.299	1.295	1.286	1.269	1.243
1.307	1.247	1.276	1.279	1.300	1.290	1.309	1.275	1.254	1.258	1.270	1.259	1.258	1.266	1.290	1.306	1.302	1.284	1.275	1.259
1.310	1.280	1.297	1.279	1.288	1.295	1.302	1.274	1.266	1.269	1.283	1.240	1.263	1.255	1.277	1.283	1.259	1.281	1.254	1.276
1.300	1.289	1.277	1.263	1.272	1.283	1.272	1.291	1.255	1.249	1.275	1.227	1.259	1.251	1.236	1.265	1.238	1.249	1.222	1.264
1.300	1.286	1.300	1.256	1.244	1.252	1.278	1.293	1.269	1.265	1.278	1.223	1.250	1.255	1.246	1.261	1.269	1.254	1.239	1.286
1.291	1.299	1.303	1.264	1.279	1.273	1.301	1.315	1.274	1.258	1.270	1.216	1.241	1.252	1.252	1.272	1.271	1.267	1.222	1.290
1.290	1.299	1.307	1.259	1.287	1.286	1.302	1.302	1.262	1.267	1.294	1.229	1.252	1.261	1.259	1.290	1.302	1.304	1.249	1.305
1.294	1.287	1.300	1.268	1.301	1.280	1.286	1.304	1.265	1.238	1.282	1.257	1.261	1.278	1.261	1.278	1.280	1.272	1.259	1.317
1.300	1.297	1.300	1.262	1.304	1.281	1.282	1.262	1.249	1.231	1.255	1.247	1.239	1.276	1.264	1.261	1.258	1.279	1.245	1.314
1.296	1.300	1.314	1.295	1.312	1.285	1.278	1.264	1.297	1.263	1.275	1.265	1.286	1.291	1.291	1.295	1.294	1.298	1.264	1.295
1.300	1.295	1.315	1.283	1.312	1.300	1.296	1.267	1.303	1.280	1.274	1.258	1.297	1.307	1.291	1.293	1.287	1.294	1.270	1.303
1.283	1.274	1.298	1.274	1.304	1.281	1.275	1.274	1.292	1.271	1.279	1.256	1.297	1.306	1.286	1.290	1.285	1.273	1.257	1.299
1.291	1.283	1.292	1.261	1.294	1.268	1.270	1.280	1.267	1.263	1.291	1.251	1.291	1.295	1.292	1.288	1.300	1.279	1.268	1.300
1.278	1.269	1.288	1.254	1.273	1.288	1.288	1.272	1.256	1.256	1.262	1.249	1.272	1.279	1.312	1.306	1.302	1.287	1.284	1.305
																			
1.622	1.570	1.598	1.685	1.681	1.606	1.663	1.672	1.664	1.682	1.658	1.652	1.679	1.662	1.628	1.662	1.655	1.638	1.631	1.633
1.673	1.604	1.659	1.677	1.680	1.645	1.666	1.680	1.684	1.681	1.679	1.671	1.642	1.602	1.581	1.637	1.626	1.587	1.583	1.613
1.665	1.629	1.666	1.684	1.675	1.666	1.681	1.693	1.659	1.694	1.686	1.638	1.595	1.653	1.588	1.639	1.608	1.613	1.600	1.626
1.647	1.697	1.684	1.700	1.672	1.616	1.667	1.641	1.650	1.633	1.640	1.654	1.599	1.644	1.641	1.639	1.596	1.662	1.625	1.620
1.645	1.632	1.650	1.667	1.632	1.608	1.644	1.612	1.629	1.644	1.625	1.607	1.617	1.672	1.653	1.632	1.602	1.615	1.654	1.608
1.610	1.630	1.632	1.632	1.648	1.630	1.678	1.637	1.662	1.694	1.634	1.599	1.613	1.670	1.680	1.652	1.635	1.642	1.634	1.640
1.649	1.624	1.632	1.646	1.650	1.638	1.649	1.646	1.627	1.655	1.623	1.604	1.641	1.641	1.651	1.656	1.638	1.667	1.664	1.601
1.660	1.642	1.656	1.634	1.648	1.645	1.675	1.646	1.648	1.632	1.617	1.611	1.627	1.604	1.616	1.645	1.655	1.654	1.673	1.628
1.646	1.653	1.655	1.623	1.652	1.653	1.668	1.646	1.643	1.628	1.619	1.624	1.591	1.601	1.621	1.624	1.642	1.660	1.656	1.647
1.631	1.671	1.636	1.618	1.628	1.625	1.631	1.638	1.612	1.635	1.606	1.585	1.589	1.579	1.606	1.612	1.624	1.641	1.631	1.655
1.634	1.654	1.662	1.598	1.626	1.650	1.652	1.646	1.637	1.636	1.617	1.621	1.600	1.578	1.615	1.624	1.616	1.659	1.633	1.651
1.638	1.647	1.650	1.602	1.626	1.631	1.675	1.676	1.635	1.636	1.608	1.600	1.624	1.612	1.607	1.616	1.607	1.646	1.593	1.653
1.642	1.640	1.673	1.610	1.621	1.640	1.670	1.678	1.667	1.664	1.659	1.603	1.654	1.624	1.610	1.644	1.622	1.643	1.604	1.660
1.663	1.641	1.654	1.598	1.620	1.628	1.634	1.652	1.644	1.657	1.655	1.567	1.610	1.639	1.635	1.646	1.646	1.658	1.608	1.672
1.676	1.630	1.641	1.619	1.647	1.633	1.654	1.637	1.619	1.618	1.618	1.562	1.630	1.635	1.643	1.625	1.636	1.640	1.590	1.660
1.682	1.661	1.648	1.658	1.660	1.645	1.669	1.638	1.623	1.640	1.656	1.603	1.642	1.672	1.673	1.655	1.662	1.657	1.613	1.681
1.684	1.677	1.672	1.626	1.648	1.661	1.685	1.650	1.619	1.613	1.650	1.633	1.639	1.655	1.649	1.638	1.683	1.657	1.613	1.660
1.628	1.644	1.672	1.609	1.663	1.646	1.684	1.645	1.618	1.604	1.636	1.648	1.642	1.660	1.662	1.673	1.696	1.652	1.636	1.658
1.629	1.656	1.683	1.597	1.650	1.649	1.670	1.645	1.614	1.601	1.644	1.631	1.635	1.660	1.668	1.662	1.669	1.682	1.623	1.681
1.640	1.668	1.651	1.605	1.640	1.637	1.671	1.641	1.616	1.602	1.634	1.623	1.649	1.669	1.669	1.653	1.696	1.674	1.624	1.681
																			
2.321	2.184	2.222	2.275	2.334	2.320	2.336	2.338	2.318	2.415	2.379	2.324	2.331	2.321	2.340	2.361	2.342	2.314	2.341	2.286
2.323	2.349	2.212	2.305	2.387	2.292	2.364	2.354	2.271	2.359	2.297	2.315	2.317	2.263	2.268	2.356	2.291	2.314	2.324	2.269
2.260	2.313	2.413	2.374	2.327	2.224	2.318	2.361	2.298	2.389	2.327	2.345	2.339	2.291	2.285	2.366	2.327	2.280	2.304	2.335
2.266	2.310	2.370	2.330	2.255	2.281	2.369	2.328	2.332	2.361	2.396	2.286	2.321	2.289	2.343	2.364	2.343	2.286	2.334	2.331
2.255	2.309	2.324	2.346	2.230	2.229	2.276	2.327	2.317	2.333	2.311	2.264	2.268	2.346	2.345	2.371	2.329	2.297	2.321	2.310
2.300	2.398	2.306	2.316	2.238	2.258	2.346	2.321	2.309	2.293	2.293	2.230	2.288	2.309	2.337	2.349	2.303	2.315	2.265	2.282
2.315	2.403	2.269	2.238	2.266	2.322	2.339	2.328	2.322	2.317	2.298	2.295	2.359	2.353	2.399	2.426	2.333	2.394	2.302	2.380
2.276	2.372	2.282	2.242	2.278	2.296	2.275	2.298	2.324	2.347	2.319	2.257	2.271	2.325	2.427	2.347	2.324	2.344	2.336	2.373
2.297	2.344	2.269	2.207	2.254	2.299	2.290	2.324	2.325	2.343	2.340	2.271	2.316	2.282	2.378	2.335	2.327	2.360	2.372	2.370
2.287	2.263	2.263	2.216	2.231	2.266	2.273	2.291	2.300	2.314	2.323	2.356	2.306	2.277	2.356	2.297	2.298	2.317	2.299	2.280
2.318	2.299	2.287	2.235	2.298	2.339	2.286	2.339	2.295	2.355	2.346	2.334	2.315	2.315	2.348	2.369	2.346	2.343	2.317	2.329
2.293	2.338	2.269	2.300	2.236	2.343	2.341	2.329	2.298	2.345	2.356	2.338	2.310	2.308	2.360	2.312	2.351	2.312	2.316	2.332
2.309	2.336	2.243	2.307	2.286	2.327	2.334	2.322	2.280	2.342	2.304	2.327	2.257	2.283	2.360	2.351	2.347	2.371	2.301	2.330
2.278	2.312	2.289	2.278	2.277	2.310	2.283	2.326	2.318	2.340	2.317	2.295	2.310	2.283	2.347	2.340	2.292	2.303	2.319	2.333
2.303	2.269	2.257	2.213	2.275	2.281	2.317	2.351	2.277	2.339	2.333	2.307	2.275	2.313	2.360	2.302	2.292	2.302	2.344	2.334
2.319	2.352	2.319	2.252	2.276	2.288	2.338	2.332	2.312	2.290	2.319	2.328	2.277	2.288	2.306	2.346	2.331	2.345	2.318	2.371
2.314	2.378	2.304	2.284	2.309	2.284	2.311	2.375	2.348	2.318	2.312	2.373	2.268	2.299	2.311	2.335	2.316	2.355	2.317	2.357
2.361	2.334	2.352	2.315	2.319	2.294	2.293	2.344	2.349	2.302	2.297	2.317	2.288	2.314	2.345	2.340	2.337	2.395	2.341	2.322
2.312	2.309	2.366	2.330	2.285	2.285	2.288	2.339	2.323	2.296	2.277	2.303	2.336	2.312	2.374	2.366	2.328	2.363	2.419	2.361
2.327	2.366	2.336	2.334	2.297	2.352	2.331	2.369	2.298	2.309	2.315	2.318	2.298	2.303	2.327	2.330	2.317	2.331	2.360	2.311

fixed vec[rec] fixcv_encf(3)(20,20)
0.322	0.449	0.670	0.795	0.942	1.074	1.190	1.305	1.378	1.421	1.462	1.786	2.086	2.319	2.552	2.738	2.933	3.070	3.235	4.442
0.500	0.701	1.003	1.247	1.473	1.621	1.782	1.931	2.047	2.128	2.193	2.703	3.060	3.436	3.835	4.219	4.457	4.654	4.905	6.898
0.663	0.882	1.322	1.607	1.860	2.065	2.286	2.417	2.642	2.779	2.954	3.533	3.990	4.475	4.976	5.340	5.604	5.995	6.217	8.693
0.761	1.066	1.525	1.797	2.075	2.375	2.601	2.797	2.962	3.124	3.342	4.054	4.707	5.305	5.894	6.331	6.604	7.049	7.327	10.706
0.851	1.201	1.671	2.020	2.359	2.641	2.930	3.125	3.354	3.535	3.721	4.595	5.423	6.123	6.757	7.108	7.524	7.927	8.384	12.001
0.922	1.301	1.882	2.277	2.640	3.015	3.322	3.504	3.802	3.967	4.214	5.050	5.988	6.732	7.421	7.789	8.236	8.846	9.484	13.212
1.030	1.440	2.078	2.531	2.906	3.230	3.517	3.771	4.084	4.277	4.520	5.500	6.422	7.153	7.983	8.495	9.165	9.553	10.028	14.345
1.106	1.559	2.208	2.719	3.125	3.515	3.822	4.089	4.360	4.564	4.854	5.904	6.828	7.758	8.555	9.139	9.880	10.396	10.998	15.396
1.189	1.659	2.411	2.893	3.308	3.767	4.093	4.370	4.600	4.937	5.226	6.322	7.256	8.143	8.977	9.723	10.243	10.846	11.544	16.737
1.252	1.754	2.486	3.018	3.450	3.915	4.282	4.585	4.830	5.157	5.527	6.526	7.610	8.531	9.301	10.064	10.646	11.392	11.990	17.206
1.316	1.868	2.616	3.178	3.567	4.053	4.565	4.927	5.099	5.461	5.726	6.939	8.050	9.090	9.873	10.872	11.486	12.167	12.820	18.417
1.389	1.973	2.750	3.310	3.819	4.340	4.840	5.173	5.387	5.681	5.970	7.301	8.440	9.397	10.424	11.258	11.877	12.795	13.197	19.281
1.443	2.049	2.880	3.413	4.034	4.506	4.968	5.376	5.640	5.979	6.267	7.568	8.868	9.902	10.977	11.978	12.762	13.362	14.111	20.668
1.497	2.099	2.962	3.590	4.248	4.623	5.107	5.569	5.878	6.213	6.602	7.864	9.370	10.461	11.306	12.249	13.087	13.962	14.730	21.511
1.544	2.203	3.097	3.711	4.367	4.829	5.201	5.631	5.992	6.283	6.716	8.100	9.526	10.839	11.763	12.539	13.435	14.611	15.106	22.094
1.618	2.278	3.275	3.980	4.623	5.046	5.435	5.920	6.373	6.654	7.078	8.607	10.216	11.401	12.362	13.234	14.063	15.073	15.543	22.889
1.642	2.361	3.397	4.019	4.758	5.244	5.681	6.192	6.590	6.933	7.125	8.862	10.517	11.871	12.761	13.781	14.553	15.617	16.223	23.610
1.700	2.380	3.421	4.137	4.841	5.376	5.858	6.339	6.840	7.150	7.494	9.099	10.953	12.179	13.097	14.208	14.842	16.029	16.577	24.234
1.742	2.449	3.487	4.225	4.919	5.459	6.042	6.540	6.928	7.294	7.582	9.345	11.158	12.502	13.686	14.535	15.696	16.469	17.227	24.990
1.777	2.520	3.586	4.344	5.021	5.675	6.207	6.665	7.085	7.489	7.787	9.597	11.351	12.736	14.124	15.133	16.214	17.034	17.920	26.120
																			
0.489	0.698	1.037	1.231	1.477	1.622	1.792	1.881	2.089	2.187	2.272	2.804	3.118	3.563	3.977	4.255	4.446	4.715	4.984	7.009
0.729	1.002	1.486	1.790	2.049	2.307	2.545	2.733	2.948	3.071	3.249	3.963	4.557	4.993	5.571	5.797	6.394	6.706	6.966	9.843
0.901	1.227	1.799	2.234	2.618	2.925	3.240	3.477	3.742	3.808	4.042	4.790	5.584	6.316	6.805	7.277	7.763	8.254	8.728	12.615
1.015	1.433	2.119	2.527	2.886	3.284	3.561	3.856	4.161	4.341	4.429	5.464	6.392	6.958	7.857	8.585	9.103	9.535	10.044	14.630
1.135	1.582	2.295	2.807	3.185	3.520	3.999	4.267	4.524	4.808	5.044	6.100	7.124	8.171	8.781	9.442	9.944	10.415	11.151	16.007
1.251	1.738	2.506	3.051	3.509	3.980	4.409	4.690	5.033	5.393	5.658	6.797	7.945	8.950	9.789	10.499	11.188	11.583	12.317	17.617
1.395	1.910	2.729	3.303	3.872	4.290	4.847	5.098	5.393	5.861	6.048	7.297	8.396	9.543	10.470	11.230	11.900	12.527	13.508	19.239
1.497	2.032	2.905	3.548	4.223	4.696	5.121	5.421	5.800	6.178	6.495	7.818	9.026	10.064	11.156	12.022	12.819	13.814	14.636	20.584
1.552	2.205	3.145	3.744	4.511	5.013	5.433	5.937	6.243	6.620	6.829	8.314	9.582	10.545	11.772	12.599	13.703	14.483	15.401	22.077
1.607	2.312	3.238	3.896	4.602	5.213	5.679	6.099	6.412	6.778	7.084	8.682	9.986	11.069	12.266	13.354	14.257	15.123	16.201	23.007
1.709	2.405	3.398	4.049	4.814	5.501	5.958	6.507	6.858	7.125	7.456	9.217	10.409	11.778	13.246	14.100	15.184	16.401	17.081	24.601
1.775	2.534	3.614	4.273	5.058	5.715	6.274	6.705	7.126	7.494	7.868	9.670	10.975	12.266	13.753	14.897	15.757	16.765	17.684	25.853
1.876	2.628	3.771	4.479	5.338	5.987	6.536	7.107	7.388	7.888	8.444	10.029	11.717	13.282	14.547	15.733	16.785	17.811	18.353	26.565
1.987	2.729	3.862	4.598	5.478	6.145	6.710	7.283	7.677	8.119	8.656	10.452	12.132	13.670	15.099	16.018	17.400	18.283	18.921	27.987
2.039	2.791	4.055	4.873	5.695	6.265	6.779	7.365	7.763	8.315	8.875	10.749	12.582	14.121	15.635	16.528	17.887	18.831	19.299	28.547
2.125	2.914	4.204	5.060	5.879	6.532	7.122	7.757	8.169	8.698	9.306	11.174	13.222	14.812	16.226	17.405	18.831	19.796	19.919	29.394
2.178	3.064	4.380	5.244	6.079	6.719	7.310	7.973	8.345	8.896	9.615	11.535	13.526	15.082	16.681	17.898	19.181	20.398	21.036	30.226
2.218	3.102	4.490	5.373	6.158	6.866	7.602	8.110	8.646	9.068	9.783	12.118	13.973	15.676	17.012	18.286	19.825	20.829	21.616	31.178
2.288	3.221	4.551	5.486	6.400	7.114	7.775	8.324	8.872	9.316	10.101	12.366	14.263	15.976	17.685	18.881	20.490	21.729	22.207	32.003
2.307	3.337	4.725	5.578	6.577	7.381	7.999	8.550	9.085	9.753	10.382	12.598	14.693	16.164	18.280	19.592	20.966	22.180	22.884	33.130
																			
0.908	1.300	1.941	2.421	2.813	3.069	3.307	3.550	3.821	4.016	4.196	5.445	6.059	6.537	7.129	7.910	8.470	8.951	9.420	12.111
1.154	1.669	2.411	2.987	3.579	3.870	4.233	4.535	4.777	4.922	5.082	6.236	7.169	8.275	8.574	9.971	10.767	11.466	12.211	16.493
1.469	1.965	2.935	3.492	4.060	4.531	4.914	5.360	5.603	6.027	6.309	7.769	9.098	10.536	11.058	12.314	12.913	13.569	14.365	19.846
1.610	2.236	3.168	3.978	4.658	5.256	5.613	6.232	6.570	7.085	7.333	8.863	10.268	11.723	13.053	13.929	14.982	15.774	15.984	22.361
1.758	2.474	3.344	4.190	5.013	5.431	5.961	6.493	6.891	7.536	7.675	9.671	11.242	12.544	13.794	15.222	15.965	17.047	17.757	24.897
1.982	2.672	3.783	4.697	5.423	6.191	6.633	7.287	7.653	8.278	8.791	10.250	12.024	13.389	15.567	16.025	17.208	18.573	19.019	28.221
2.120	2.893	4.068	4.995	5.937	6.812	7.392	7.747	8.163	8.943	9.495	11.235	13.602	15.349	17.187	17.942	19.307	20.224	20.373	29.608
2.218	2.983	4.348	5.394	6.321	7.294	7.906	8.329	8.730	9.012	9.632	12.132	13.904	15.597	17.513	19.115	19.868	21.187	22.172	31.593
2.301	3.163	4.623	5.551	6.673	7.330	8.101	8.879	9.003	9.876	10.113	12.625	14.714	16.345	18.464	19.917	21.353	22.935	23.454	32.494
2.392	3.267	4.753	5.839	6.753	7.617	8.328	9.160	9.706	10.487	10.794	13.424	15.384	17.124	19.621	20.310	21.758	23.933	24.878	34.674
2.545	3.478	4.889	6.012	7.279	8.266	9.032	9.904	10.190	11.103	11.557	14.151	16.403	17.844	20.306	22.005	22.815	24.891	25.915	36.344
2.672	3.628	5.080	6.215	7.441	8.536	9.283	10.258	10.797	11.182	11.696	14.415	17.013	19.507	20.987	22.670	23.460	25.371	26.603	38.042
2.744	3.911	5.493	6.496	7.671	8.652	9.550	10.409	11.287	11.809	12.291	15.025	17.510	19.991	21.400	23.559	24.827	26.226	27.535	39.303
2.830	4.053	5.653	6.883	8.016	9.042	9.994	10.800	11.713	12.188	12.518	15.426	18.205	20.642	21.899	24.708	25.284	27.473	28.474	39.827
2.951	4.136	5.697	6.938	8.201	9.365	10.613	11.191	12.466	12.685	13.095	16.074	18.523	21.420	22.762	24.851	26.314	27.643	29.613	40.600
3.079	4.356	6.039	7.214	8.566	9.546	10.779	11.524	12.774	13.003	13.385	16.709	19.238	21.613	23.622	25.253	27.299	28.803	31.094	43.024
3.161	4.455	6.390	7.478	8.706	10.136	11.164	11.760	13.011	13.300	14.118	17.548	19.948	23.021	24.444	26.241	28.167	29.617	32.237	44.777
3.245	4.580	6.539	7.707	8.851	10.238	11.310	12.098	13.271	13.452	14.141	17.978	20.302	23.217	25.509	27.263	28.586	30.646	32.681	45.847
3.215	4.644	6.513	7.858	9.086	10.233	11.282	12.332	13.562	13.794	14.731	18.368	20.954	23.660	25.809	28.369	29.001	30.909	33.151	47.480
3.419	4.851	6.817	8.324	9.550	10.689	11.803	12.769	13.774	14.005	14.886	19.048	21.324	24.264	26.148	28.530	30.507	31.560	34.417	49.622

dim me(2) mae(2) rmse(2) mse(2)
dim fc(2) error(2) error2(2)
clear fc error error2

ENTER(varying) res_ser
ENTER(varying) unres_ser

**** determining sample over which all data are available
inquire(regressorlist) stser>>start endser>>end
# LHS res_ser unres_ser
comp stshared = stser
comp endshared = endser

if fcst <= stshared .or. fcst > endshared
{
dis '!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
dis '!!!! Inappropriate argument for "fcst"                            !!!!!'
dis '!!!! Please make sure that "fcst" falls within the sample-period. !!!!!' 
dis '!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
dis 'Procedure Terminated'
branch done
}

*******
*******Step 1. generating forecast for the restricted model
*******

*** NOTE: to speed computations, model estimates are obtained using the
*** recursive updating algorithm described in such sources as
*** Judge, Griffiths, Hill, Lutkepohl, and Lee (1985, p. 174), rather
*** than by re-estimating the regression at each forecast date

linreg(noprint) LHS stshared fcst-1
# res_ser
comp sig1 = %xx
comp R = %nobs

do time = fcst,endshared
 **** updating, adding an observation (and storing forecast)
 comp xvec = %eqnxvector(0,time)
 comp proj1 = %dot(xvec,%beta)
 comp fc(1)(time) = proj1
 if scheme<>1
  {
   **** adding an observation
   comp midterm = %qform(sig1,xvec)
   comp xs = sig1*xvec
   comp sig = sig1 - (1./(midterm+1.))*%outerxx(xs)
   comp %beta = %beta + (sig*xvec)*(LHS(time)-proj1)
   comp sig1 = sig
      **** dropping an observation if scheme=rolling
   if scheme==2
    {
     comp xvec = %eqnxvector(0,time-R)
     comp proj1 = %dot(xvec,%beta)
     comp midterm = %qform(sig1,xvec)
     comp xs = sig1*xvec
     comp sig = sig1 + (1./(1.-midterm))*%outerxx(xs)
     comp %beta = %beta - (sig*xvec)*((LHS(time-R))-proj1)
     comp sig1 = sig
    }
    
  }
end do time

*******
*******Step 2. generating forecast for the unrestricted model
*******
linreg(noprint) LHS stshared fcst-1
# res_ser unres_ser
comp sig1 = %xx
comp R = %nobs

**** checking that models are appropriately nested, by checking for singularity in (x'x)^-1
if %det(%xx)==0.
 {
  dis ' '
  dis '!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
  dis '!!!! failure to properly specify nested models                    !!!!!'
  dis '!!!! make sure card 1 includes variables in null model            !!!!!' 
  dis '!!!! and card 2 includes the additional vars  of the altern model !!!!!' 
  dis '!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
  dis 'Procedure Terminated'
  branch done
 }

do time = fcst,endshared
 **** updating, adding an observation (and storing forecast)
 comp xvec = %eqnxvector(0,time)
 comp proj1 = %dot(xvec,%beta)
 comp fc(2)(time) = proj1
 if scheme<>1
  {
   **** adding an observation
   comp midterm = %qform(sig1,xvec)
   comp xs = sig1*xvec
   comp sig = sig1 - (1./(midterm+1.))*%outerxx(xs)
   comp %beta = %beta + (sig*xvec)*(LHS(time)-proj1)
   comp sig1 = sig
      **** dropping an observation if scheme=rolling
   if scheme==2
    {
     comp xvec = %eqnxvector(0,time-R)
     comp proj1 = %dot(xvec,%beta)
     comp midterm = %qform(sig1,xvec)
     comp xs = sig1*xvec
     comp sig = sig1 + (1./(1.-midterm))*%outerxx(xs)
     comp %beta = %beta - (sig*xvec)*((LHS(time-R))-proj1)
     comp sig1 = sig
    }    
  }
end do time
*******
*******Step 3. Computing summary stats and MSE-F, MSE-t, ENC-F, and ENC-t statistics.
*******
do i =1,2
   set error(i) fcst endshared =  LHS{0} - fc(i){0} 
   set error2(i) fcst endshared = error(i){0}**2
   set abserror fcst endshared = abs(error(i){0})
   stats(noprint) error(i) fcst endshared
   comp me(i) = %mean
   stats(noprint) abserror fcst endshared
   comp mae(i) = %mean
   stats(noprint) error2(i) fcst endshared
   comp mse(i) = %mean
   comp rmse(i) = sqrt(mse(i))
   if i == 1
   {
   dis ' '
   dis '***Restricted model'   
   dis 'In-sample' %Datelabel(stshared) '-' %Datelabel(fcst-1) 
   dis 'Out-sample' %Datelabel(fcst) '-' %Datelabel(endshared) 
   }
   if i == 2
   {
   dis '***Unrestricted model'   
   dis 'In-sample' %Datelabel(stshared) '-' %Datelabel(fcst-1) 
   dis 'Out-sample' %Datelabel(fcst) '-' %Datelabel(endshared) 
   }
   dis @15 'ME' @27 'MAE' @39 'RMSE' @51 'MSE' 
   dis @11 '--------' @24 '--------' @37 '--------' @49 '--------' 
   dis @7 ######.##### me(i) ######.##### mae(i) ######.##### rmse(i) ######.##### mse(i)
end do i

set dser fcst endshared = error2(1) - error2(2)
set cser fcst endshared = error2(1) - error(1){0}*error(2){0}

**** MSE-F and MSE-t
linreg(noprint) dser fcst endshared
# constant
comp msef = %nobs*(%beta(1)/mse(2))
comp mset = %tstats(1)

**** ENC-F and ENC-t
linreg(noprint) cser fcst endshared
# constant
comp encf = %nobs*(%beta(1)/mse(2))
comp enct = %tstats(1)

if %defined(fc_usr)
 {
  dim fc_usr(2)
  do i=1,2
   set fc_usr(i) = fc(i)
  end do i
 }

if %defined(error_usr)
 {
  dim error_usr(2)
  do i=1,2
   set error_usr(i) = error(i) 
  end do i
 }

*******
*******Step 4. pulling appropriate critical-values
*******

comp cvpis = ||.1,.2,.4,.6,.8,1,1.2,1.4,1.6,1.8,2,3,4,5,6,7,8,9,10,20|| ;* pi vals for which crit vals are available

comp maxnsig = 3           ;* number of significance levels for which crit vals are available in the data files
                            * order:  (1) .10, (2) .05, (3) .01
comp maxnpi = %rows(cvpis) ;* number of pi vals for which crit vals are available
comp maxnk = 20            ;* maximum number of d.o.f. (k2) for which critical values are available

****** figuring out pi value
comp P = endshared - fcst + 1 ;* Out-of-sample observations
comp R = (fcst-1) - stshared + 1  ;* In-sample observations
comp pival = float(P)/float(R)
set assignser 1 maxnpi = abs(cvpis(t) - pival)
extremum(noprint) assignser 1 maxnpi
comp cmpos1 = %minent
if pival < cvpis(cmpos1)
 comp cmpos2 = %minent - 1
else 
 comp cmpos2 = %minent + 1
 
dis '*****'
dis 'Actual pi val = ' @35 pival   
dis ':nearest CM pi val =' @35 cvpis(cmpos1)
dis ':second nearest CM pi val =' @35 cvpis(cmpos2) 
dis '*****'

comp ndf = %rlcount(unres_ser)
dis 'number of restrictions in test =' ndf                             
                            
***** set-up here: one matrix for each significance level; dimension of matrix is max(k2) by # pi values
***** all matrices are stored in a text file, with one file for each forecast scheme
dim cv_msef(maxnsig) cv_mset(maxnsig) cv_encf(maxnsig) cv_enct(maxnsig)
do i = 1,maxnsig
  dim cv_msef(i)(maxnk,maxnpi)
  dim cv_mset(i)(maxnk,maxnpi)
  dim cv_encf(i)(maxnk,maxnpi)
  dim cv_enct(i)(maxnk,maxnpi)
end do i

do i = 1,maxnsig
 if scheme==1
  {
   comp cv_mset(i) = fixcv_mset(i)
   comp cv_msef(i) = fixcv_msef(i)
   comp cv_enct(i) = fixcv_enct(i)
   comp cv_encf(i) = fixcv_encf(i)
  }
 else if scheme==2
  {
   comp cv_mset(i) = rollcv_mset(i)
   comp cv_msef(i) = rollcv_msef(i)
   comp cv_enct(i) = rollcv_enct(i)
   comp cv_encf(i) = rollcv_encf(i)
  }
 else
  {
   comp cv_mset(i) = reccv_mset(i)
   comp cv_msef(i) = reccv_msef(i)
   comp cv_enct(i) = reccv_enct(i)
   comp cv_encf(i) = reccv_encf(i)
  }
end do i

******************** pulling relevant critical value for each test
******************** using linear interpolation between the critical values associated
******************** with the pi values surrounding P/R
comp pidist = cvpis(cmpos1)-cvpis(cmpos2)
comp msefcv10 = cv_msef(1)(ndf,cmpos1)*(pival-cvpis(cmpos2))/(pidist)+cv_msef(1)(ndf,cmpos2)*(cvpis(cmpos1)-pival)/(pidist)
comp msefcv05 = cv_msef(2)(ndf,cmpos1)*(pival-cvpis(cmpos2))/(pidist)+cv_msef(2)(ndf,cmpos2)*(cvpis(cmpos1)-pival)/(pidist)
comp msefcv01 = cv_msef(3)(ndf,cmpos1)*(pival-cvpis(cmpos2))/(pidist)+cv_msef(3)(ndf,cmpos2)*(cvpis(cmpos1)-pival)/(pidist)

comp msetcv10 = cv_mset(1)(ndf,cmpos1)*(pival-cvpis(cmpos2))/(pidist)+cv_mset(1)(ndf,cmpos2)*(cvpis(cmpos1)-pival)/(pidist)
comp msetcv05 = cv_mset(2)(ndf,cmpos1)*(pival-cvpis(cmpos2))/(pidist)+cv_mset(2)(ndf,cmpos2)*(cvpis(cmpos1)-pival)/(pidist)
comp msetcv01 = cv_mset(3)(ndf,cmpos1)*(pival-cvpis(cmpos2))/(pidist)+cv_mset(3)(ndf,cmpos2)*(cvpis(cmpos1)-pival)/(pidist)

comp encfcv10 = cv_encf(1)(ndf,cmpos1)*(pival-cvpis(cmpos2))/(pidist)+cv_encf(1)(ndf,cmpos2)*(cvpis(cmpos1)-pival)/(pidist)
comp encfcv05 = cv_encf(2)(ndf,cmpos1)*(pival-cvpis(cmpos2))/(pidist)+cv_encf(2)(ndf,cmpos2)*(cvpis(cmpos1)-pival)/(pidist)
comp encfcv01 = cv_encf(3)(ndf,cmpos1)*(pival-cvpis(cmpos2))/(pidist)+cv_encf(3)(ndf,cmpos2)*(cvpis(cmpos1)-pival)/(pidist)
     
comp enctcv10 = cv_enct(1)(ndf,cmpos1)*(pival-cvpis(cmpos2))/(pidist)+cv_enct(1)(ndf,cmpos2)*(cvpis(cmpos1)-pival)/(pidist)
comp enctcv05 = cv_enct(2)(ndf,cmpos1)*(pival-cvpis(cmpos2))/(pidist)+cv_enct(2)(ndf,cmpos2)*(cvpis(cmpos1)-pival)/(pidist)
comp enctcv01 = cv_enct(3)(ndf,cmpos1)*(pival-cvpis(cmpos2))/(pidist)+cv_enct(3)(ndf,cmpos2)*(cvpis(cmpos1)-pival)/(pidist)


*******
*******Step 5. Displaying the results.
*******

dis '**************************** test results ****************************'
dis ' '
dis
dis 'test' @8 'test value' @22 'asymptotic critical values (1-sided)'
dis  @25 '10%' @34 '5%' @42 '1%'
dis ' '
dis 'MSF-F' ###.### @10 msef @22 msefcv10 msefcv05 msefcv01
dis 'MSF-t' ###.### @10 mset @22 msetcv10 msetcv05 msetcv01
dis 'ENC-F' ###.### @10 encf @22 encfcv10 encfcv05 encfcv01
dis 'ENC-t' ###.### @10 enct @22 enctcv10 enctcv05 enctcv01

if print==1
 {
  dis @28 'restricted' @46 'unrestricted'
  dis @28 'model fc' @46 'model fc'
  print fcst endshared LHS fc(1) fc(2)
 }

:done

end proc
