[petsc-users] AMG preconditioners

[Manav Bhatia]

Thanks! 

I am intrigued by this. 

I am solving a system of equations where the natural block size is 2x2. This is obtained by rewriting a complex linear system of equations in real numbers. 

I use MatSetValuesBlocked to set the values of each 2x2 block in my matrix. I am assuming that this method works only for blocked matrix. So, to make the code generic, I will need to use the explicit row/column numberings along with a MatSetValues method. 

-Manav


> On Mar 9, 2016, at 8:59 AM, Barry Smith <bsmith at mcs.anl.gov> wrote:
> 
> 
>> On Mar 9, 2016, at 8:43 AM, Manav Bhatia <bhatiamanav at gmail.com> wrote:
>> 
>> Hi, 
>> 
>>   Out of GAMG, ML and HYPRE, which are expected to be compatible with SEQBAIJ and MPIBAIJ matrices? 
> 
>  Not necessarily.  Note that you can create matrices with MatCreate() and control at run time the matrix version. You don't need to hardwire BAIJ or AIJ (in fact we recommend not hardwiring it). 
> 
>  Barry
> 
> 
>> 
>> Regards,
>> Manav
>> 
>