c
c  ____________________________
c
	subroutine symm_apply(eee,f,index)
	real r(3,3,48)
	real eee(3,3),f(3,3)
	integer index,ia,ib,ic
c  matrix multiply, for symmetry operations
c
	do 1940, ia=1,3
	do 1940, ib=1,3
	f(ia,ib)=0.
	do 1940, ic=1,3
	f(ia,ib)=f(ia,ib)+eee(ia,ic)*r(ib,ic,index)
c  postmultiply by transpose of the symmetry operator (c1)
c  
1940	continue
c
	return
	end
c
c  ____________________________
c