c
c  ____________________________
c
	subroutine symm_pre_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)+r(ia,ic,index)*eee(ic,ib)
c  premultiply by the symmetry operator (c1)
c  
1940	continue
c
	return
	end
c
c
c  -----------------------------
c