anchor.diff:      Testval = gen(Alg),
       case Refval =:= Testval of
           true -> ok;
  +        false when Refval =:= not_implemented ->
  +            exit({not_implemented,Alg});
           false ->
   	    io:format("Failed: ~p~n",[Alg]),
   	    io:format("Length ~p ~p~n",[length(Refval), length(Testval)]),
  
positive.code:      17936751184378118743, 4224632875737239207, 15888641556987476199, 9586888813112229805,
       9476861567287505094, 14909536929239540332, 17996844556292992842, 2699310519182298856];
  
  reference_jump_val(exs64) -> [not_implemented].
  
negative.code: %% Check that the algorithms generate sound values.
  
  basic_stats_uniform_1(Config) when is_list(Config) ->
      ct:timetrap({minutes,15}), %% valgrind needs a lot of time
      [basic_uniform_1(?LOOP, rand:seed_s(Alg), 0.0, array:new([{default, 0}]))
       || Alg <- algs()],
      ok.
  
  basic_stats_uniform_2(Config) when is_list(Config) ->
      ct:timetrap({minutes,15}), %% valgrind needs a lot of time
      [basic_uniform_2(?LOOP, rand:seed_s(Alg), 0, array:new([{default, 0}]))
       || Alg <- algs()],
      ok.
  
  basic_stats_standard_normal(Config) when is_list(Config) ->
      ct:timetrap({minutes,6}), %% valgrind needs a lot of time
      io:format("Testing standard normal~n",[]),
      IntendedMean = 0,
      IntendedVariance = 1,
      [basic_normal_1(?LOOP, IntendedMean, IntendedVariance,
                      rand:seed_s(Alg), 0, 0)
       || Alg <- algs()],
      ok.
  
  basic_stats_normal(Config) when is_list(Config) ->
      IntendedMeans = [-1.0e6, -50, -math:pi(), -math:exp(-1),
                       0.12345678, math:exp(1), 100, 1.0e6],
      IntendedVariances = [1.0e-6, math:exp(-1), 1, math:pi(), 1.0e6],
      IntendedMeanVariancePairs =
          [{Mean, Variance} || Mean <- IntendedMeans,
                               Variance <- IntendedVariances],
  
      ct:timetrap({minutes, 6 * length(IntendedMeanVariancePairs)}), %% valgrind needs a lot of time
      lists:foreach(
        fun ({IntendedMean, IntendedVariance}) ->
                ct:pal(
                  "Testing normal(~.2f, ~.2f)~n",
                  [float(IntendedMean), float(IntendedVariance)]),
                [basic_normal_1(?LOOP, IntendedMean, IntendedVariance,
                                rand:seed_s(Alg), 0, 0)
                 || Alg <- algs()]
        end,
        IntendedMeanVariancePairs).
  
  basic_uniform_1(N, S0, Sum, A0) when N > 0 ->
      {X,S} = rand:uniform_s(S0),
  