cgetRotMat4.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cgetrotmat4 (m, ang, rm)
subroutine	cinvrotmat4 (rm, rmn)
subroutine	cmultrotmat4 (a, b, c)
subroutine	capplyrot4 (a, v, vn)

Function/Subroutine Documentation

capplyrot4()

subroutine capplyrot4	(	real*8, dimension(4,4)	a,
		type(fmom)	v,
		type(fmom)	vn
	)

Definition at line 148 of file cgetRotMat4.f.

Referenced by clorep(), and cuonsphere().

!          3-d transformation matrix a is multiplied by
!          a vector v to obtain a new vector vn.
!          vn can  be v.
!
           implicit none
!----           include '../Zptcl.h'
#include  "Zptcl.h"
           type(fmom):: v, vn
           real*8 a(4,4)
           type(fmom):: tmp
!
           real*8 sum
           integer i, j
!
          do  i=1, 3
             sum=0.
             do  j=1, 4
                 sum=sum + a(i, j)*v%p(j)
             enddo    
             tmp%p(i)=sum
          enddo
          vn = tmp

Here is the caller graph for this function:

cgetrotmat4()

subroutine cgetrotmat4	(	integer	m,
		real*8	ang,
		real*8, dimension(4, 4)	rm
	)

Definition at line 79 of file cgetRotMat4.f.

References cos, and d0.

Referenced by clorep(), and cuonsphere().

          implicit none
!
          integer m
          real*8 ang
          real*8 rm(4, 4)
!
          integer i, j, m1, m2
          real*8 c, s
!
          if(m .ge. 1 .and. m .le. 3) then
             do i=1, 4
                do j=1, 4
                   rm(i, j)= 0.d0
                enddo
             enddo   
             rm(4,4)=1.d0
             rm(m, m)=1.d0
             m1=mod(m,3)+1
             m2=mod(m1,3)+1
             c=cos(ang)
             s=sin(ang)
             rm(m1,m1)=c
             rm(m2,m2)=c
             rm(m1,m2)=s
             rm(m2,m1)=-s
          else
             write(*,*) ' invalid m=',m,' to cgetRotMat4 '
             stop  9999
          endif

Here is the caller graph for this function:

cinvrotmat4()

subroutine cinvrotmat4	(	real*8, dimension(4,4)	rm,
		real*8, dimension(4,4)	rmn
	)

Definition at line 110 of file cgetRotMat4.f.

Referenced by clorep().

!             Invert rotation matrix rm and put into rmn.
!             rm should be a roation matrix made by calling
!             cgetRotMat4 (with ang).  rmn can be made by calling
!             cgetRotMat4 with -ang, too.  this one is to avoid
!             computing the cos and sin for reducing time.
!             rmn cannot be the same arrays as rm.
!             rmn is nothing but the transposed matrix of rm.
          implicit none
          real*8 rm(4,4), rmn(4,4)
          integer i, j
          do  i=1,4
             do  j=i+1, 4
                rmn(i,j)=rm(j, i)
                rmn(j,i)=rm(i, j)
             enddo
          enddo
          do  i=1, 4
             rmn(i,i)=rm(i, i)
          enddo

Here is the caller graph for this function:

cmultrotmat4()

subroutine cmultrotmat4	(	real*8, dimension(4,4)	a,
		real*8, dimension(4,4)	b,
		real*8, dimension(4,4)	c
	)

Definition at line 131 of file cgetRotMat4.f.

Referenced by clorep(), and cuonsphere().

!            3-d matrix product c=a*b
!            c cannot be either of a or b.
          implicit none
          real*8 a(4,4), b(4,4), c(4,4)
          integer i, j, k
          real*8 ab
          do  i=1, 4
             do j=1, 4
                ab=0.
                do k=1, 4
                   ab=ab+ a(i,k)*b(k,j)
                enddo
                c(i,j)=ab
             enddo   
          enddo

Here is the caller graph for this function: