 
 =#=#=#=#=#=#=#=#=#=#=kppSolve=#=#=#=#=#=#=#=#=#=#=
 
-__device__ void ros_Solve(double *Ghimj, double *K, int &Nsol, int istage,  int ros_S)
+__device__ void ros_Solve(double * __restrict__ Ghimj, double * __restrict__ K, int &Nsol, const int istage, const int ros_S)
 {
 
     #pragma unroll 4 

}
__device__ double ros_ErrorNorm(double *var, double * __restrict__ varNew, double * __restrict__ varErr, 
                                const double * __restrict__ absTol, const double * __restrict__ relTol,
                                int vectorTol )
{
    double err, scale, varMax;


                                     //  cuda global mem buffers              
                                     d_absTol, d_relTol,   
                                     // extra - vector lenght and processor
                                     VL_GLO);

    GPU_DEBUG();

