Oleg, let me give you a more detailed explanation.
1. SOBOL QRNG is implemented so that if params array is incorrectly filled the default implementation with dimension 1 will be used.
2. soboluserdirnums.f example shows how to initialize SOBOL QRNG with a complete direction numbers table. If you'd like you provide the generator with both initial direction numbers and primitive polynomials youll need to slightly change params array filling like this:
params(1) = dm1
params(2) = VSL_USER_QRNG_INITIAL_VALUES
params(3) = IOR(VSL_USER_INIT_DIRECTION_NUMBERS, &
& VSL_USER_PRIMITIVE_POLYMS)
params(3) = IOR(params(3), VSL_QRNG_OVERRIDE_1ST_DIM_INIT)
do i = 1, dm1
params(i+3) = poly(i)
end do
params(4 + dm1) = maxdeg
k = 5 + dm1
do i=1,dm1
do j=1,maxdeg
params(k) = vinit(i,j)
k = k + 1
end do
end do
nparams = 4 + dm1 + dm1*maxdeg
