cgetRotMat3.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	cgetrotmat3 (m, cosa, sina, rm)
subroutine	cinvrotmat3 (rm, rmn)
subroutine	cmultrotmat3 (a, b, c)
subroutine	capplyrot3 (a, v, vn)

Function/Subroutine Documentation

capplyrot3()

subroutine capplyrot3	(	real*8, dimension(3,3)	a,
		real*8, dimension(3)	v,
		real*8, dimension(3)	vn
	)

Definition at line 145 of file cgetRotMat3.f.

!          3-d transformation matrix a is multiplied by
!          a vector v to obtain a new vector vn.
!          vn can  be v.
!
           implicit none
           real*8 a(3,3), v(3), vn(3)
           
!
           real*8 sum
           integer i, j
!
          do  i=1, 3
             sum=0.
             do  j=1, 3
                 sum=sum + a(i, j)*v(j)
             enddo    
             vn(i)=sum
          enddo

cgetrotmat3()

subroutine cgetrotmat3	(	integer	m,
		real*8	cosa,
		real*8	sina,
		real*8, dimension(3, 3)	rm
	)

Definition at line 80 of file cgetRotMat3.f.

References cerrormsg(), and d0.

          implicit none
!
          integer m
          real*8 cosa, sina
          real*8 rm(3, 3)
!
          integer i, j, m1, m2
          character*70 msg
!
          if(m .ge. 1 .and. m .le. 3) then
             do i=1, 3
                do j=1, 3
                   rm(i, j)= 0.d0
                enddo
             enddo   
             rm(m, m)=1.d0
             m1=mod(m,3)+1
             m2=mod(m1,3)+1
             rm(m1,m1)=cosa
             rm(m2,m2)=cosa
             rm(m1,m2)=sina
             rm(m2,m1)=-sina
          else
             write(msg, *) ' invalid m=',m,' to cgetRotMat3 '
             call cerrormsg(msg,0)
          endif

Here is the call graph for this function:

cinvrotmat3()

subroutine cinvrotmat3	(	real*8, dimension(3,3)	rm,
		real*8, dimension(3,3)	rmn
	)

Definition at line 108 of file cgetRotMat3.f.

!             Invert rotation matrix rm and put into rmn.
!             rm should be a roation matrix made by calling
!             cgetRotMat3 (with cosa, sina).  rmn can be made by calling
!             cgetRotMat3 with cosa, -sina, too.  
!             rmn cannot be the same arrays as rm.
!             rmn is nothing but the transposed matrix of rm.
          implicit none
          real*8 rm(3,3), rmn(3,3)
          integer i, j
          do  i=1,3
             do  j=i+1, 3
                rmn(i,j)=rm(j, i)
                rmn(j,i)=rm(i, j)
             enddo
          enddo
          do  i=1, 3
             rmn(i,i)=rm(i, i)
          enddo

cmultrotmat3()

subroutine cmultrotmat3	(	real*8, dimension(3,3)	a,
		real*8, dimension(3,3)	b,
		real*8, dimension(3,3)	c
	)

Definition at line 128 of file cgetRotMat3.f.

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