Re: [petsc-users] MATSETVALUES: Fortran problem

[Matthew Knepley]

On Fri, Mar 15, 2024 at 9:55 AM Frank Bramkamp  wrote:

> Dear PETSc Team, I am using the latest petsc version 3. 20. 5. I would
> like to create a matrix using MatCreateSeqAIJ To insert values, I use
> MatSetValues. It seems that the Fortran interface/stubs are missing for
> MatsetValues, as the linker does
> ZjQcmQRYFpfptBannerStart
> This Message Is From an External Sender
> This message came from outside your organization.
>
> ZjQcmQRYFpfptBannerEnd
>
> Dear PETSc Team,
>
> I am using the latest petsc version 3.20.5.
>
>
> I would like to create a matrix using
> MatCreateSeqAIJ
>
> To insert values, I use MatSetValues.
> It seems that the Fortran interface/stubs are missing for MatsetValues, as 
> the linker does not find any subroutine with that name.
> MatSetValueLocal seems to be fine.
>
>
Here is a Fortran example calling MatSetValues():


https://urldefense.us/v3/__https://gitlab.com/petsc/petsc/-/blob/main/src/ksp/ksp/tutorials/ex1f.F90?ref_type=heads__;!!G_uCfscf7eWS!ag7gfsg4qE_beGvp2_VPt5PyN1bZcFr8xkVlbukEpUnqYHx_awYDepWWZfT-7rtVF5-lHC2GXj1ETtSIaOLc$
 


> Typically I am using a blocked matrix format (BAIJ), which works fine in 
> fortran.
> Soon we want to try PETSC on GPUs, using the format MATAIJCUSPARSE, since 
> there seems not to be a blocked format available in PETSC for GPUs so far.
>
> You can use the blocked input API, like MatSetValuesBlocked(), with all
the storage formats. There is no explicit blocking in the MATAIJCUSPARSE format
because Nvidia handles the optimization differently.

  Thanks,

 Matt

Therefore I first want to try the pointwise format MatCreateSeqAIJ
format on a CPU, before using the GPU format.
>
> I think that CUDA also supports a block format now ?! Maybe that would be 
> also useful to have one day.
>
>
> Greetings, Frank Bramkamp
>
>
>
>
>
>
>
>

-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which their
experiments lead.
-- Norbert Wiener

https://urldefense.us/v3/__https://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!ag7gfsg4qE_beGvp2_VPt5PyN1bZcFr8xkVlbukEpUnqYHx_awYDepWWZfT-7rtVF5-lHC2GXj1ETryTqsT3$
  

Re: [petsc-users] MATSETVALUES: Fortran problem

[Barry Smith]

> On Mar 15, 2024, at 9:53 AM, Frank Bramkamp  wrote:
> 
> This Message Is From an External Sender
> This message came from outside your organization.
> Dear PETSc Team,
> 
> I am using the latest petsc version 3.20.5.
> 
> 
> I would like to create a matrix using
> MatCreateSeqAIJ
> 
> To insert values, I use MatSetValues.
> It seems that the Fortran interface/stubs are missing for MatsetValues, as 
> the linker does not find any subroutine with that name.
> MatSetValueLocal seems to be fine.

   Please send the exact error message (cut and paste), there are definitely 
Fortran stubs for this function but it could be you exact parameter input does 
not have a stub yet.

   Barry

> 
> 
> Typically I am using a blocked matrix format (BAIJ), which works fine in 
> fortran.
> Soon we want to try PETSC on GPUs, using the format MATAIJCUSPARSE, since 
> there seems not to be a blocked format available in PETSC for GPUs so far.
> Therefore I first want to try the pointwise format MatCreateSeqAIJ format on 
> a CPU, before using the GPU format.
> 
> I think that CUDA also supports a block format now ?! Maybe that would be 
> also useful to have one day.
> 
> 
> Greetings, Frank Bramkamp
> 
> 
> 
> 
> 
> 
> 

[petsc-users] MATSETVALUES: Fortran problem

[Frank Bramkamp]

 Dear PETSc Team, I am using the latest petsc version 3. 20. 5. I would like to create a matrix using MatCreateSeqAIJ To insert values, I use MatSetValues. It seems that the Fortran interface/stubs are missing for MatsetValues, as the linker does




ZjQcmQRYFpfptBannerStart




  

  
	This Message Is From an External Sender
  
  
This message came from outside your organization.
  



 
  


ZjQcmQRYFpfptBannerEnd




Dear PETSc Team,

I am using the latest petsc version 3.20.5.


I would like to create a matrix using
MatCreateSeqAIJ

To insert values, I use MatSetValues.
It seems that the Fortran interface/stubs are missing for MatsetValues, as the linker does not find any subroutine with that name.
MatSetValueLocal seems to be fine.


Typically I am using a blocked matrix format (BAIJ), which works fine in fortran.
Soon we want to try PETSC on GPUs, using the format MATAIJCUSPARSE, since there seems not to be a blocked format available in PETSC for GPUs so far.
Therefore I first want to try the pointwise format MatCreateSeqAIJ format on a CPU, before using the GPU format.

I think that CUDA also supports a block format now ?! Maybe that would be also useful to have one day.


Greetings, Frank Bramkamp