CHARACTER*40 INFILE,FILENAME
	CHARACTER*80 STRINGA,STRINGB
	character*3 type,id
	LOGICAL EFLG
	real*8 M1(3,4),M2(3,4),M3(3,4)
	character*40 string
	INTEGER*4 MM,NN,I,J

1200	FORMAT(A)


c 	read operator into m2
6	WRITE(6,*) 'ENTER NAME OF rotation-trans op (.o):'
	READ(5,1200) INFILE
	IF (INFILE(:5).EQ.'     ') STOP
	  INQUIRE(FILE=INFILE, NAME=FILENAME,EXIST=EFLG)
	  IF (.NOT.EFLG) GOTO 6
        OPEN (UNIT=2,FILE=INFILE,STATUS='OLD',READONLY)
c	READ(2,*) '.SPACE_GROUP_operators    r   12 (3f12.7)'
	READ(2,'(A40)') string
	do 60 j=1,4
60	read(2,*)(M2(I,J),I=1,3)
	CLOSE(UNIT=2)
	write(6,*)
	MM=3
	NN=4
c	CALL MATPRN (M2,MM,NN)
	IATOM=1

c 	read vector into m1
5	WRITE(6,*)'ENTER coordinates of point to be operated on:'
50	read(5,'(A)',end=1160) STRINGA
	I=INDEX(STRINGA,' ')
	J=I+INDEX(STRINGA(I+1:),' ')
	ID=STRINGA(I:J)
	STRINGB=STRINGA(J+1:)
c	write(6,*)'i=',i,' j=',j,' id=',id
c	write(6,*)stringb

	DECODE (80,52,STRINGB) M1(1,1)
c	J=INDEX(STRINGB,' ')
	stringb=stringb(1+index(stringb,' '):)
	DECODE (80,52,STRINGB) M1(2,1)
	stringb=stringb(1+index(stringb,' '):)
	DECODE (80,52,STRINGB) M1(3,1)
	stringb=stringb(1+index(stringb,' '):)
	type=stringb(1:3)
	stringb=stringb(1+index(stringb,' '):)
	DECODE (80,52,STRINGB) occ
	stringb=stringb(1+index(stringb,' '):)
	DECODE (80,52,STRINGB) B

c	DECODE (80,52,STRINGB) (M1(I,1),I=1,3),type,occ,b
52	format(3f9.6,A3,2f9.6)
C50	read(5,*,end=1160) (m1(i,1),I=1,3)
c	write(6,*) ID, (m1(i,1),I=1,3),type,occ,b
c	CALL MATPRN (M1,MM,NN)

c	Multiply m1 matrix and transl by m2 matrix (first 3 cols), put im M3

	do 100 i=1,3
	do 100 j=1,1
	m3(i,j)=0
	do 100 k=1,3
	m3(i,j)=m3(i,j)+m2(i,k)*m1(k,j) 
C	WRITE(6,*)I,J,K,M2(I,K),M1(K,J),M3(I,J)
100	CONTINUE
c	now add the translation of the second operator:
	do 110 i=1,3
110	m3(i,1)=m3(i,1)+m2(i,4)
	WRITE(6,51) iatom,type,ID,iatom,(M3(I,1),i=1,3),occ,b
51 	format (' ATOM  ',I5,1x,A3,1x,A3,' A',i5,4x,3f8.3,2f6.2)
9	FORMAT(' ',3F14.6)
	IATOM=IATOM+1
	goto 50
1150	continue
1160	continue
	END

	SUBROUTINE MATPRN(ARRAY,M,N)
	REAL*8 ARRAY(3,4)
	WRITE(6,*)'M, N=',M,N
	DO 10 I=1,M
10	WRITE (6,50) (ARRAY(I,J),J=1,N)
50	format(8f10.5)
  	RETURN
        END