A bug in dcg_get?

[Zbigniew_L_]

I wrote a small code calling the RCI CG routine with the ILU0 preconditioner. The code is attached. The code returns an incorrect vector of all zeros as a solution. I have noticed that the vector tmp contains the correct solution but it is not transfered to the vector x (= solution).

Zbigniew

[Gennady_F_Intel]

is that linux or windows?

[Alexander_K_Intel2]

Hi,

I see that you call dcg_get to return correct value of solution, but dcg_get return only current number of iteration

Output Parameters

itercount

INTEGER. Returns the current iteration number.

Thanks,

Alex

[Zbigniew_L_]

It is Linux (64-bit).

[Zbigniew_L_]

Here is my output:

tmp:
tmp[0]=-0.0909091, tmp[1]=-0.181818, tmp[2]=-0.272727, tmp[3]=-0.363636, tmp[4]=-0.454545, tmp[5]=-0.545455, tmp[6]=-0.636364, tmp[7]=-0.727273, tmp[8]=-0.818182, tmp[9]=-0.909091, 
# of iterations = 0
solution x:
x[0]=0, x[1]=0, x[2]=0, x[3]=0, x[4]=0, x[5]=0, x[6]=0, x[7]=0, x[8]=0, x[9]=0, 

It is evident that the solution (x) is not correct.

Zbigniew

[Zbigniew_L_]

I should have asked the question "A bug in dcg?" because it is dcg which returns an incorrect solution (all zeros).

Also, why the number of iterations is 0?

Zbigniew