include 'mod_numerique.f'

c******************************************************************

	program test_colloc
	
c	test pour colloc

c	Auteur: P. Morel, Departement J.D. Cassini, O.C.A., Observatoire de Nice
c	CESAM95
	
c-----------------------------------------------------------------

	use mod_kind
	use mod_numerique
	
	real (kind=dp), allocatable, dimension(:) :: x, xt, xx, derxx, lderxx,
	1	xl 
	
	integer :: m, r, n, knot, i, rp1, ne, ncoll, nlim, dim, ord
		
c-------------------------------------------------------------------

2000	format(1x,8es10.3)

	n=20 ; m=2 ; r=1 ; ne=3
	ncoll=(n-1)*m ; dim=ncoll+r ; nlim=ne*r ; ord=m+r ; rp1=r+1
	
	allocate(x(n),xt(dim+ord),xl(nlim),xx(ncoll),derxx(ord*ncoll*rp1),
	1	lderxx(rp1*ord*nlim))	
	do i=1,n
	 x(i)=real(i-1)/(n-1)
	enddo
	
	xl=x(1)		!limites en x(1)

        call noedif(x,n,m,r,xt,knot)
         
	call colloc(m,r,n,x,xt,knot,xx,derxx,nlim,lderxx,xl)
	print*,'xx'
	write(*,2000)xx
	print*,'derxx'
	write(*,2000)derxx
	print*,'lderxx'
	write(*,2000)lderxx
	
	stop
	
	end program test_colloc