Hi,
I am using DMDAVecGetArrayF90 with a pointer in petsc 3.4.3. Is there anyway of
finding out what the start and finish indices of the pointer coming back from
this routine are per processor?
I have declared:
DM :: da
Vec :: X
PetscScalar,pointer :: xv(:,:)=>null()
DMDAVecGetArrayF90(da,X,xv,ierr)
Do j = 1, m
Do i = 1, n
xv(i-1.j-1) = rhs2d(i,j)
enddo
enddo
DMDAVecRestoreArrayF90(da,X,xv,ierr)
When we have xv(0:n-1,0:m-1) from DMDAVecGetArrayF90 it works fine but other
processors may give back eg; xv(0:n-1,1:m-1) or something else and it will fail
as the second index
now goes below 0 to 1. Any ideas of what I can do next?
Thanks,
Carol
Dr Carol Brickley
BSc,PhD,ARCS,DIC,MBCS
Senior Software Engineer
Applied Computer Science
DS+T,
AWE
Aldermaston
Reading
Berkshire
RG7 4PR
Direct: 0118 9855035
___________________________________________________
____________________________
The information in this email and in any attachment(s) is
commercial in confidence. If you are not the named addressee(s)
or
if you receive this email in error then any distribution, copying or
use of this communication or the information in it is strictly
prohibited. Please notify us immediately by email at
admin.internet(at)awe.co.uk, and then delete this message from
your computer. While attachments are virus checked, AWE plc
does not accept any liability in respect of any virus which is not
detected.
AWE Plc
Registered in England and Wales
Registration No 02763902
AWE, Aldermaston, Reading, RG7 4PR
