PROGRAM TESTV1
C
C	  PROGRAM TO TEST FOLLOWING ROUTINES ON MATLIB
C	  VADDD, VCLRD, VCOMD, VMAGD, VMOVD, VSCLD, VSUBD
C	  (ALL ARE DOUBLE PRECISION ROUTINES)
C
	REAL*8 A(5),B(5),C(5),P,Q
	DATA P,Q/ 5.,2./
C
	DO 10 I=1,5
	A(I) = I
	B(I) = 2*I
10	CONTINUE
	TYPE *,' DOUBLE PRECISON VECTOR TEST'
	TYPE 20,A,B,P,Q
20	FORMAT(' A=',5F5.1,/' B=',5F5.1,/' P='F5.1,' Q='F5.1)
	CALL VADDD(C,A,B,4)
	TYPE *,' VADDD ',C
	CALL VCLRD(C,4)
	TYPE *,' VCLRD ',C
	CALL VCOMD(C,P,A,Q,B,4)
	TYPE *,' VCOMD ',C
	C(1) = VMAGD(A,4)
	TYPE *,' VMAGD ',C(1)
	CALL VMOVD(C,A,4)
	TYPE *,' VMOVD ',C
	CALL VSCLD(C,P,A,4)
	TYPE *,' VSCLD ',C
	CALL VSUBD(C,A,B,4)
	TYPE *,' VSUBD ',C
	DO 30 I=1,5
 30	C(I) = 0.
	TYPE *,' 3-VECTOR TEST'
	CALL V3ADDD(C,A(3),B(3))
	TYPE *,' V3ADDD ',C
	CALL V3CLRD(C)
	TYPE *,' V3CLRD ',C
	CALL V3COMD(C,P,A(3),Q,B(3))
	TYPE *,' V3COMD ',C
	C(1) = V3MAGD(A(3))
	TYPE *,' V3MAGD ',C(1)
	CALL V3SCLD(C,P,A(3))
	TYPE *,' V3SCLD ',C
	CALL V3SUBD(C,A(3),B(3))
	TYPE *,' V3SUBD ',C
	STOP
	END