CG method implementation

[Hiroyuki Kasai]

Dear all,

Reading the code of CG solver (conjugategradient.m), I found the code below;

            Pnewgrad = newgrad;

            newgradPnewgrad = inner_product(newgrad, Pnewgrad);  

From seeing the entire code, I cannot understand that why "Pnewgrad" is needed.

It would be appreciated to tell me the purpose of this "Pnewgrad"  variable.

Regards,

Hiro

[BM]

Hello Hiro,

Which version are you using? In the current version, this line is not there. Could you verify it again?

Cheers,

BM

[Hiroyuki Kasai]

Hello Bamdev,

This is completely my mistake. I am sorry for that.

Now, I understood that when no preconditioner provided, Pnewgrad is the same as newgrad.

Is this correct?

Regards,

Hiro

[Nicolas Boumal]

Hello Hiro,

That is right: if you specify no preconditioner, than the preconditioner is the identity map, so that Pnewgrad = newgrad.

Best,

Nicolas