09-08-2010 01:21 AM
In order to switch-off iterative refinement in DSS, please addMKL_DSS_REFINEMENT_OFF constant to the opt parameter before calling dss_solve_rela or dss_solve_complex routine.
opt = MKL_DSS_REFINEMENT_OFF;
error = dss_solve_real( handle, opt, rhs, nRhs, solValues );
opt = MKL_DSS_DEFAULTS; // this is needed if you use opt variable later in the code for other dss calls
09-08-2010 05:06 AM
thanks for reply.
I am getting one problem, in order to solve equation of Ax=B type i am using DSS which is running fine. But when compared with Gauss Siedel method for same thing then results does not match.
09-08-2010 11:31 PM
i plotted the results along a line in the given domain. I got an error of 10% when compared with results obtained from gauss siedel solver.
I also wanted to know one more thing, if in a given coefficient 'A' (of Ax=b) i want to apply Neumann boundary condition at one end Then what necessary changes should i make in my coefficient matrix.
09-09-2010 09:27 PM
If you have zero Neumann boundary condition then nothing change both initial matrix and rhs, if non zero Neumann boundary condition only right hand side change (This is true if you construct you matrix using FEM, FVM or some method like it, in other case we need to see initial differential equation and way of descritization to answer your question.).
With best regards,