[petsc-users] Error during code compile with MATGETOWNERSHIPRANGE

Barry Smith bsmith at petsc.dev
Thu Aug 22 09:39:08 CDT 2024


  Hmm, try using a standalone variable

  PetscReal value 
  value = int_impl(k,5)]


>>>>> call MatSetValues(A_mat_uv,[ione],II,[ione],[value],impl_mat_A,INSERT_VALUES,ierr)
>>>>> 

unfortunately Fortan compilers in this situation are not good at telling us exactly what argument is giving it grief.

Barry



> On Aug 22, 2024, at 10:33 AM, TAY Wee Beng <zonexo at gmail.com> wrote:
> 
> 
> 
> On 22/8/2024 10:28 pm, Barry Smith wrote:
>> 
>>   Should be PetscReal as it is a numerical value.
> Ok, I changed it but I still get the same error.
>> 
>> 
>>> On Aug 22, 2024, at 9:55 AM, TAY Wee Beng <zonexo at gmail.com> <mailto:zonexo at gmail.com> wrote:
>>> 
>>> 
>>> 
>>> On 22/8/2024 9:54 pm, Barry Smith wrote:
>>>> 
>>>>   What is int_impl(k,5) defined type?
>>> PetscInt
>>>> 
>>>>> On Aug 22, 2024, at 9:45 AM, TAY Wee Beng <zonexo at gmail.com> <mailto:zonexo at gmail.com> wrote:
>>>>> 
>>>>> Hi Barry,
>>>>> 
>>>>> Do you mean that I change from:
>>>>> 
>>>>> call MatSetValues(A_mat_uv,ione,II,ione,int_impl(k,5),impl_mat_A,INSERT_VALUES,ierr)
>>>>> 
>>>>> to
>>>>> 
>>>>> call MatSetValues(A_mat_uv,[ione],II,[ione],[int_impl(k,5)],impl_mat_A,INSERT_VALUES,ierr)
>>>>> 
>>>>> ?
>>>>> 
>>>>> I did it but the error is still there.
>>>>> 
>>>>> On 22/8/2024 9:40 pm, Barry Smith wrote:
>>>>>> 
>>>>>>    Fortran 90 type checking is very tight; The dimension of the array, or scalar passed as arguments must match the expected dimension (f77 did not do this type checking). Thus the ione argument must be a 1-d array as well as the numerical values so do 
>>>>>> 
>>>>>>> call MatSetValues(A_mat_uv,[ione],II,[ione],[int_impl(k,5)],impl_mat_A,INSERT_VALUES,ierr)
>>>>>>> 
>>>>>> 
>>>>>> See Fortran at https://urldefense.us/v3/__https://petsc.org/main/changes/dev/__;!!G_uCfscf7eWS!b2f2s3CzjxVoPuzYvafOupk882EeGhmMhGIGBixXIP3SRNSBy6JkxfMOXE_aqVmSFgCwaEIhK4WHFyybLs7x4K8$ 
>>>>>> 
>>>>>> I am trying to support the old-fashion F77 model, allowing miss-matches in the array dimensions while still doing proper type checking but it will take some time to simplify the API.
>>>>>> 
>>>>>>    Barry
>>>>>> 
>>>>>> 
>>>>>> 
>>>>>>> On Aug 21, 2024, at 9:44 PM, TAY Wee Beng <zonexo at gmail.com> <mailto:zonexo at gmail.com> wrote:
>>>>>>> 
>>>>>>> Hi Barry,
>>>>>>> 
>>>>>>> I have declared them as integers in Fortran. Is that different from PetscInt and how come it works in debug mode?
>>>>>>> 
>>>>>>> Anyway, I changed them and it solved the problem. However, I have a similar problem in my boundary.F90:
>>>>>>> 
>>>>>>> boundary.F90(6685): error #6285: There is no matching specific subroutine for this generic subroutine call.   [MATSETVALUES]
>>>>>>> call MatSetValues(A_mat_uv,ione,II,ione,int_impl(k,5),impl_mat_A,INSERT_VALUES,ierr)
>>>>>>> -----^
>>>>>>> I changed all to PetscInt and also PetscReal but I still got the error.
>>>>>>> 
>>>>>>> Why is this so now? Any solution?
>>>>>>> 
>>>>>>> Thanks!
>>>>>>> 
>>>>>>> On 22/8/2024 12:03 am, Barry Smith wrote:
>>>>>>>> 
>>>>>>>>   You must declare as
>>>>>>>> 
>>>>>>>>   PetscInt ksta_p,kend_p
>>>>>>>> 
>>>>>>>>   Perhaps they are declared as arrays?
>>>>>>>> 
>>>>>>>> 
>>>>>>>>> On Aug 21, 2024, at 11:19 AM, TAY Wee Beng <zonexo at gmail.com> <mailto:zonexo at gmail.com> wrote:
>>>>>>>>> 
>>>>>>>>> Hi,
>>>>>>>>> 
>>>>>>>>> I am using the latest PETSc thru github. I compiled both the debug and rel ver of PETSc w/o problem.
>>>>>>>>> 
>>>>>>>>> I then use it with my CFD code and the debug ver works.
>>>>>>>>> 
>>>>>>>>> However, I have problems with the rel ver:
>>>>>>>>> 
>>>>>>>>> ftn -o global.o -c -O3 -g -ip -ipo   -fPIC  -save -w -I/home/project/11003851/lib/petsc_210824_intel_rel/include    global.F90
>>>>>>>>> ifort: remark #10448: Intel(R) Fortran Compiler Classic (ifort) is now deprecated and will be discontinued late 2024. Intel recommends that customers transition now to using the LLVM-based Intel(R) Fortran Compiler (ifx) for continued Windows* and Linux* support, new language support, new language features, and optimizations. Use '-diag-disable=10448' to disable this message.
>>>>>>>>> global.F90(444): error #6285: There is no matching specific subroutine for this generic subroutine call.   [MATGETOWNERSHIPRANGE]
>>>>>>>>>         call MatGetOwnershipRange(A_mat,ksta_p,kend_p,ierr)
>>>>>>>>> -------------^
>>>>>>>>> global.F90(720): error #6285: There is no matching specific subroutine for this generic subroutine call.   [MATGETOWNERSHIPRANGE]
>>>>>>>>> call MatGetOwnershipRange(A_mat_uv,ksta_m,kend_m,ierr)
>>>>>>>>> -----^
>>>>>>>>> global.F90(774): error #6285: There is no matching specific subroutine for this generic subroutine call.   [MATGETOWNERSHIPRANGE]
>>>>>>>>> call MatGetOwnershipRange(A_semi_x,ksta_mx,kend_mx,ierr)
>>>>>>>>> -----^
>>>>>>>>> global.F90(776): error #6285: There is no matching specific subroutine for this generic subroutine call.   [MATGETOWNERSHIPRANGE]
>>>>>>>>> call MatGetOwnershipRange(A_semi_y,ksta_my,kend_my,ierr)
>>>>>>>>> -----^
>>>>>>>>> global.F90(949): error #6285: There is no matching specific subroutine for this generic subroutine call.   [MATGETOWNERSHIPRANGE]
>>>>>>>>> call MatGetOwnershipRange(A_semi_x,ksta_mx,kend_mx,ierr)
>>>>>>>>> -----^
>>>>>>>>> global.F90(957): error #6285: There is no matching specific subroutine for this generic subroutine call.   [MATGETOWNERSHIPRANGE]
>>>>>>>>> call MatGetOwnershipRange(A_semi_y,ksta_my,kend_my,ierr)
>>>>>>>>> -----^
>>>>>>>>> compilation aborted for global.F90 (code 1)
>>>>>>>>> 
>>>>>>>>> May I know what's the problem?
>>>>>>>>> 
>>>>>>>>> --
>>>>>>>>> 
>>>>>>>>> Thank you very much.
>>>>>>>>> 
>>>>>>>>> Yours sincerely,
>>>>>>>>> 
>>>>>>>>> ================================================
>>>>>>>>> TAY Wee-Beng 郑伟明 (Zheng Weiming)
>>>>>>>>> ================================================
>>>>>>>>> 
>>>>>>>>> 
>>>>>>>> 
>>>>>>> --
>>>>>>> 
>>>>>>> Thank you very much.
>>>>>>> 
>>>>>>> Yours sincerely,
>>>>>>> 
>>>>>>> ================================================
>>>>>>> TAY Wee-Beng 郑伟明 (Zheng Weiming)
>>>>>>> ================================================
>>>>>>> 
>>>>>>> 
>>>>>> 
>>>>> --
>>>>> 
>>>>> Thank you very much.
>>>>> 
>>>>> Yours sincerely,
>>>>> 
>>>>> ================================================
>>>>> TAY Wee-Beng 郑伟明 (Zheng Weiming)
>>>>> ================================================
>>>>> 
>>>>> 
>>>> 
>>> --
>>> 
>>> Thank you very much.
>>> 
>>> Yours sincerely,
>>> 
>>> ================================================
>>> TAY Wee-Beng 郑伟明 (Zheng Weiming)
>>> ================================================
>>> 
>>> 
>> 
> --
> 
> Thank you very much.
> 
> Yours sincerely,
> 
> ================================================
> TAY Wee-Beng 郑伟明 (Zheng Weiming)
> ================================================
> 
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20240822/37e48066/attachment-0003.html>


More information about the petsc-users mailing list