Commit e976945b authored by Zolnierczuk, Piotr's avatar Zolnierczuk, Piotr
Browse files

small bug in testpolynom

parent 5b2ff03e

Makefile.version

+1 −1
Original line number Diff line number Diff line
@@ -3,7 +3,7 @@ export PROJECT=drspine 


export VERSION_MAJOR=2

export VERSION_MINOR=0

export VERSION_RELEASE=0-a7

export VERSION_RELEASE=0-a8



export PROJLIB=lib$(PROJECT).a

export PROJARCH=$(PROJECT)-$(VERSION_MAJOR).$(VERSION_MINOR)

tests/testpolynom.f90

+2 −2
Original line number Diff line number Diff line
@@ -49,7 +49,7 @@ subroutine test_polyfit(x, y, n, aexp) 
    write(output_unit,'("polyfit n=",i0)') n

    do i=n,0,-1

    delta = abs((a(i)-aexp(i))/aexp(i))

    write(cmsg,'("a[",i2,"]",g14.8," ",g14.8, " d="g14.8)') i,a(i),aexp(i),delta

    write(cmsg,'("a[",i2,"]",g14.8," ",g14.8," d=",g14.8)') i,a(i),aexp(i),delta

    write(output_unit,*) trim(cmsg)

    call assert(delta<EPS,'polyfit '//trim(cmsg))

    end do
@@ -59,7 +59,7 @@ subroutine test_polyfit(x, y, n, aexp) 
    do i=1, NDATA

    yf = polyeval(x(i),a,n)

    delta = abs((y(i)-yf)/yf)

    write(cmsg,'("y[",i2,"]",g14.8," ",g14.8, " d="g14.8)') i, y(i),yf, delta

    write(cmsg,'("y[",i2,"]",g14.8," ",g14.8," d=",g14.8)') i, y(i),yf, delta

    write(output_unit,*) trim(cmsg)

    call assert(abs(y(i)-yf)<EPS,'polyeval '//trim(cmsg))

    end do