Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Python skeleton codes #3

Open
wants to merge 7 commits into
base: master
Choose a base branch
from
Open

Python skeleton codes #3

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
7 commits
Select commit Hold shift + click to select a range
7b6c4e5
Python skeleton codes
Jan 29, 2016
1e93704
Adding a README.md
atableman Jan 29, 2016
e3f7fb3
Readme
Jan 29, 2016
851d039
just-in-time compilation for Python skeleton codes
Feb 5, 2016
8877121
updating files
Feb 5, 2016
5ca4bd2
readme change
Feb 8, 2016
30ff3d5
just some changes
Feb 8, 2016
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 2 additions & 2 deletions README.md
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
# Particle-In-Cell Skeleton Codes
# Particle-In-Cell Skeleton Codes
Succinct PIC Implementation

__All source code in this repository is confidential.__
Expand Down Expand Up @@ -88,7 +88,7 @@ Here's a table showing the names of the codes and their parallelism (ES means el
| -------- | --------------------- | ----------- | --------- | ----------- |
|__Basic Serial__|
||pic1|ES|C,Fortran| None-Serial|
||pic2|ES|C,Fortran| None-Serial|
||pic2|ES|C,Fortran,Python| None-Serial|
||pic3|ES|C,Fortran| None-Serial|
||bpic1|EM|C,Fortran| None-Serial|
||bpic2|EM|C,Fortran| None-Serial|
Expand Down
Binary file added serial/pic2/cpic2
View file
Open in desktop
Binary file not shown.
Binary file added serial/pic2/cpic2.o
View file
Open in desktop
Binary file not shown.
Binary file added serial/pic2/cpush2.o
View file
Open in desktop
Binary file not shown.
Binary file added serial/pic2/dtimer.o
View file
Open in desktop
Binary file not shown.
Binary file added serial/pic2/pic2_py/cpush2.so
View file
Open in desktop
Binary file not shown.
135 changes: 135 additions & 0 deletions serial/pic2/pic2_py/cpush2_h.pyf
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,135 @@
! -*- f90 -*-
! Note: the context of this file is case sensitive.

python module cpush2 ! in
interface ! in :cpush2
subroutine cdistr2(part,vtx,vty,vdx,vdy,npx,npy,idimp,nop,nx,ny,ipbc) ! in :cpush2:cpush2_hpy.f90
real dimension(*),intent(inout) :: part
real intent(in) :: vtx
real intent(in) :: vty
real intent(in) :: vdx
real intent(in) :: vdy
integer intent(in) :: npx
integer intent(in) :: npy
integer intent(in) :: idimp
integer intent(in) :: nop
integer intent(in) :: nx
integer intent(in) :: ny
integer intent(in) :: ipbc
end subroutine cdistr2
end interface
interface ! in :cpush2
subroutine cgpush2l(part,fxy,qbm,dt,ek,idimp,nop,nx,ny,nxv,nyv,ipbc) ! in :cpush2:cpush2_hpy.f90
real dimension(*),intent(inout) :: part
real dimension(*),intent(in) :: fxy
real intent(in) :: qbm
real intent(in) :: dt
real intent(inout) :: ek
integer intent(in) :: idimp
integer intent(in) :: nop
integer intent(in) :: nx
integer intent(in) :: ny
integer intent(in) :: nxv
integer intent(in) :: nyv
integer intent(in) :: ipbc
end subroutine cgpush2l
end interface
interface ! in :cpush2
subroutine cgpost2l(part,q,qm,nop,idimp,nxv,nyv) ! in :cpush2:cpush2_hpy.f90
real dimension(*),intent(in) :: part
real dimension(*),intent(inout) :: q
real intent(in) :: qm
integer intent(in) :: nop
integer intent(in) :: idimp
integer intent(in) :: nxv
integer intent(in) :: nyv
end subroutine cgpost2l
end interface
interface ! in :cpush2
subroutine cdsortp2yl(parta,partb,npic,idimp,nop,ny1) ! in :cpush2:cpush2_hpy.f90
real dimension(*),intent(in) :: parta
real dimension(*),intent(inout) :: partb
integer dimension(*),intent(inout) :: npic
integer intent(in) :: idimp
integer intent(in) :: nop
integer intent(in) :: ny1
end subroutine cdsortp2yl
end interface
interface ! in :cpush2
subroutine ccguard2l(fxy,nx,ny,nxe,nye) ! in :cpush2:cpush2_hpy.f90
real dimension(*),intent(inout) :: fxy
integer intent(in) :: nx
integer intent(in) :: ny
integer intent(in) :: nxe
integer intent(in) :: nye
end subroutine ccguard2l
end interface
interface ! in :cpush2
subroutine caguard2l(q,nx,ny,nxe,nye) ! in :cpush2:cpush2_hpy.f90
real dimension(*),intent(inout) :: q
integer intent(in) :: nx
integer intent(in) :: ny
integer intent(in) :: nxe
integer intent(in) :: nye
end subroutine caguard2l
end interface
interface ! in :cpush2
subroutine cpois22(q,fxy,isign,ffc,ax,ay,affp,we,nx,ny,nxvh,nyv,nxhd,nyhd) ! in :cpush2:cpush2_hpy.f90
real dimension(*),intent(in) :: q
real dimension(*),intent(inout) :: fxy
integer intent(in) :: isign
complex dimension(*),intent(inout) :: ffc
real intent(in) :: ax
real intent(in) :: ay
real intent(in) :: affp
real intent(inout) :: we
integer intent(in) :: nx
integer intent(in) :: ny
integer intent(in) :: nxvh
integer intent(in) :: nyv
integer intent(in) :: nxhd
integer intent(in) :: nyhd
end subroutine cpois22
end interface
interface ! in :cpush2
subroutine cwfft2rinit(mixup,sct,indx,indy,nxhyd,nxyhd) ! in :cpush2:cpush2_hpy.f90
integer dimension(*),intent(inout) :: mixup
complex dimension(*),intent(inout) :: sct
integer intent(in) :: indx
integer intent(in) :: indy
integer intent(in) :: nxhyd
integer intent(in) :: nxyhd
end subroutine cwfft2rinit
end interface
interface ! in :cpush2
subroutine cwfft2rx(f,isign,mixup,sct,indx,indy,nxhd,nyd,nxhyd,nxyhd) ! in :cpush2:cpush2_hpy.f90
real dimension(*),intent(inout) :: f
integer :: isign
integer dimension(*),intent(inout) :: mixup
complex dimension(*),intent(inout) :: sct
integer :: indx
integer :: indy
integer :: nxhd
integer :: nyd
integer :: nxhyd
integer :: nxyhd
end subroutine cwfft2rx
end interface
interface ! in :cpush2
subroutine cwfft2r2(f,isign,mixup,sct,indx,indy,nxhd,nyd,nxhyd,nxyhd) ! in :cpush2:cpush2_hpy.f90
real dimension(*),intent(inout) :: f
integer intent(in) :: isign
integer dimension(*),intent(inout) :: mixup
complex dimension(*),intent(inout) :: sct
integer intent(in) :: indx
integer intent(in) :: indy
integer intent(in) :: nxhd
integer intent(in) :: nyd
integer intent(in) :: nxhyd
integer intent(in) :: nxyhd
end subroutine cwfft2r2
end interface
end python module cpush2

! This file was auto-generated with f2py (version:2).
! See http://cens.ioc.ee/projects/f2py2e/
Binary file added serial/pic2/pic2_py/dtimer.so
View file
Open in desktop
Binary file not shown.
15 changes: 15 additions & 0 deletions serial/pic2/pic2_py/dtimer_h.pyf
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,15 @@
! -*- f90 -*-
! Note: the context of this file is case sensitive.

python module dtimer ! in
interface ! in :dtimer
subroutine dtimer(dtime,itime,icntrl) ! in :dtimer:dtimer_hpy.f90
double precision intent(inout) :: dtime
integer dimension(4),intent(inout) :: itime
integer intent(in) :: icntrl
end subroutine dtimer
end interface
end python module dtimer

! This file was auto-generated with f2py (version:2).
! See http://cens.ioc.ee/projects/f2py2e/
Binary file added serial/pic2/pic2_py/fpush2.so
View file
Open in desktop
Binary file not shown.
209 changes: 209 additions & 0 deletions serial/pic2/pic2_py/push2_h.pyf
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,209 @@
! -*- f90 -*-
! Note: the context of this file is case sensitive.

python module fpush2 ! in
interface ! in :fpush2
subroutine distr2(part,vtx,vty,vdx,vdy,npx,npy,idimp,nop,nx,ny,ipbc) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(inout) :: part
real intent(in) :: vtx
real intent(in) :: vty
real intent(in) :: vdx
real intent(in) :: vdy
integer intent(in) :: npx
integer intent(in) :: npy
integer intent(in) :: idimp
integer intent(in) :: nop
integer intent(in) :: nx
integer intent(in) :: ny
integer intent(in) :: ipbc
end subroutine distr2
end interface
interface ! in :fpush2
subroutine gpush2l(part,fxy,qbm,dt,ek,idimp,nop,nx,ny,nxv,nyv,ipbc) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(inout) :: part
real dimension(*),intent(in) :: fxy
real intent(in) :: qbm
real intent(in) :: dt
real intent(inout) :: ek
integer intent(in) :: idimp
integer intent(in) :: nop
integer intent(in) :: nx
integer intent(in) :: ny
integer intent(in) :: nxv
integer intent(in) :: nyv
integer intent(in) :: ipbc
end subroutine gpush2l
end interface
interface ! in :fpush2
subroutine gpost2l(part,q,qm,nop,idimp,nxv,nyv) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(in) :: part
real dimension(*),intent(inout) :: q
real intent(in) :: qm
integer intent(in) :: nop
integer intent(in) :: idimp
integer intent(in) :: nxv
integer intent(in) :: nyv
end subroutine gpost2l
end interface
interface ! in :fpush2
subroutine dsortp2yl(parta,partb,npic,idimp,nop,ny1) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(in) :: parta
real dimension(*),intent(inout) :: partb
integer dimension(*),intent(inout) :: npic
integer intent(in) :: idimp
integer intent(in) :: nop
integer intent(in) :: ny1
end subroutine dsortp2yl
end interface
interface ! in :fpush2
subroutine cguard2l(fxy,nx,ny,nxe,nye) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(inout) :: fxy
integer intent(in) :: nx
integer intent(in) :: ny
integer intent(in) :: nxe
integer intent(in) :: nye
end subroutine cguard2l
end interface
interface ! in :fpush2
subroutine aguard2l(q,nx,ny,nxe,nye) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(inout) :: q
integer intent(in) :: nx
integer intent(in) :: ny
integer intent(in) :: nxe
integer intent(in) :: nye
end subroutine aguard2l
end interface
interface ! in :fpush2
subroutine pois22(q,fxy,isign,ffc,ax,ay,affp,we,nx,ny,nxvh,nyv,nxhd,nyhd) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(in) :: q
real dimension(*),intent(inout) :: fxy
integer intent(in) :: isign
complex dimension(*),intent(inout) :: ffc
real intent(in) :: ax
real intent(in) :: ay
real intent(in) :: affp
real intent(inout) :: we
integer intent(in) :: nx
integer intent(in) :: ny
integer intent(in) :: nxvh
integer intent(in) :: nyv
integer intent(in) :: nxhd
integer intent(in) :: nyhd
end subroutine pois22
end interface
interface ! in :fpush2
subroutine wfft2rinit(mixup,sct,indx,indy,nxhyd,nxyhd) ! in :fpush2:push2_hpy.f90
integer dimension(*),intent(inout) :: mixup
complex dimension(*),intent(inout) :: sct
integer intent(in) :: indx
integer intent(in) :: indy
integer intent(in) :: nxhyd
integer intent(in) :: nxyhd
end subroutine wfft2rinit
end interface
interface ! in :fpush2
subroutine wfft2rx(f,isign,mixup,sct,indx,indy,nxhd,nyd,nxhyd,nxyhd) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(inout) :: f
integer :: isign
integer dimension(*),intent(inout) :: mixup
complex dimension(*),intent(inout) :: sct
integer :: indx
integer :: indy
integer :: nxhd
integer :: nyd
integer :: nxhyd
integer :: nxyhd
end subroutine wfft2rx
end interface
interface ! in :fpush2
subroutine wfft2r2(f,isign,mixup,sct,indx,indy,nxhd,nyd,nxhyd,nxyhd) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(inout) :: f
integer intent(in) :: isign
integer dimension(*),intent(inout) :: mixup
complex dimension(*),intent(inout) :: sct
integer intent(in) :: indx
integer intent(in) :: indy
integer intent(in) :: nxhd
integer intent(in) :: nyd
integer intent(in) :: nxhyd
integer intent(in) :: nxyhd
end subroutine wfft2r2
end interface
interface ! in :fpush2
subroutine fft2rxx(f,isign,mixup,sct,indx,indy,nyi,nyp,nxhd,nyd,nxhyd,nxyhd) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(inout) :: f
integer intent(in) :: isign
integer dimension(*),intent(inout) :: mixup
complex dimension(*),intent(inout) :: sct
integer intent(in) :: indx
integer intent(in) :: indy
integer intent(in) :: nyi
integer intent(in) :: nyp
integer intent(in) :: nxhd
integer intent(in) :: nyd
integer intent(in) :: nxhyd
integer intent(in) :: nxyhd
end subroutine fft2rxx
end interface
interface ! in :fpush2
subroutine fft2rxy(f,isign,mixup,sct,indx,indy,nxi,nxp,nxhd,nyd,nxhyd,nxyhd) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(inout) :: f
integer intent(in) :: isign
integer dimension(*),intent(inout) :: mixup
complex dimension(*),intent(inout) :: sct
integer intent(in) :: indx
integer intent(in) :: indy
integer intent(in) :: nxi
integer intent(in) :: nxp
integer intent(in) :: nxhd
integer intent(in) :: nyd
integer intent(in) :: nxhyd
integer intent(in) :: nxyhd
end subroutine fft2rxy
end interface
interface ! in :fpush2
subroutine fft2r2x(f,isign,mixup,sct,indx,indy,nyi,nyp,nxhd,nyd,nxhyd,nxyhd) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(inout) :: f
integer intent(in) :: isign
integer dimension(*),intent(inout) :: mixup
complex dimension(*),intent(inout) :: sct
integer intent(in) :: indx
integer intent(in) :: indy
integer intent(in) :: nyi
integer intent(in) :: nyp
integer intent(in) :: nxhd
integer intent(in) :: nyd
integer intent(in) :: nxhyd
integer intent(in) :: nxyhd
end subroutine fft2r2x
end interface
interface ! in :fpush2
subroutine fft2r2y(f,isign,mixup,sct,indx,indy,nxi,nxp,nxhd,nyd,nxhyd,nxyhd) ! in :fpush2:push2_hpy.f90
real dimension(*),intent(inout) :: f
integer intent(in) :: isign
integer dimension(*),intent(inout) :: mixup
complex dimension(*),intent(inout) :: sct
integer intent(in) :: indx
integer intent(in) :: indy
integer intent(in) :: nxi
integer intent(in) :: nxp
integer intent(in) :: nxhd
integer intent(in) :: nyd
integer intent(in) :: nxhyd
integer intent(in) :: nxyhd
end subroutine fft2r2y
end interface
interface ! in :fpush2
function ranorm() ! in :fpush2:push2_hpy.f90
double precision :: ranorm
end function ranorm
end interface
interface ! in :fpush2
function randum() ! in :fpush2:push2_hpy.f90
double precision :: randum
end function randum
end interface
end python module fpush2

! This file was auto-generated with f2py (version:2).
! See http://cens.ioc.ee/projects/f2py2e/
Loading