Ziele Untersuchung
mit Columbo Integrität von
Datenbanken Interaktion und
Portierbarkeit Ergonomie der
Schnittstellen

Angebot Produkte Projekt Beratung

Mittel Analytik Modellierung Sprachen Algebra Logik Hardware Denken Kreativität

Zusammenhänge Gesellschaft Wirtschaft Branche Firma


Quellcodebibliothek Statistik Leitseite products/Sources/formale Sprachen/PVS/vectors/ (Beweissystem der NASA Version 6.0.9^©) Datei vom 28.9.2014 mit Größe 8 kB

SSL sqrt_approx.pvs Interaktion und
PortierbarkeitPVS

sqrt_approx : theory
begin
%-----------------------------------------------------------------------------
%
%  Includes sqrt_newton, sqrt_approx
%
%  This theory was originally developed by Michel Levy (LSR IMAG, France) after
%  Exercise Page 66 N°8 of Reasoned Programming by Krysia Broda,
%  Susan Eisenbach, Hessam Koshnevisan, and Steve Vickers, 1994. Prentice Hall.
%
%  It was modified and integrated to sqrt_approx in the NASA PVS Library
%  by C. Munoz (NIA) on March 28 2005. The definition of sqrt_bisect was
%  moved to the Interval package.
%
%-----------------------------------------------------------------------------

  IMPORTING sq,sqrt,log_nat

  x,y           : VAR nnreal
  n             : VAR nat
  z             : VAR posreal
  eps,eps1,eps2 : VAR posreal

  sq_sq_aux : lemma
    sq((z+x/z)/2)-x = sq((sq(z)-x)/(2*z))

  sqrt_newton_aux(x,(y | x < sq(y)),eps) :
    recursive {z | x < sq(z) AND sq(z)-x < eps} =
    if sq(y)-x < eps then y
    else sqrt_newton_aux(x,((y+x/y)/2),eps)
    endif
    measure if sq(y)-x < eps then 0 else 1+log_nat((sq(y)-x)/eps,4)`1 endif

  sqrt_newton_aux_NOT_increasing: LEMMA
    EXISTS (x1,x2:nnreal,(y1:nnreal|x1<sq(y1)),(y2:nnreal|x2<sq(y2)),eps):
      x1<x2 AND y1<y2 IMPLIES
        sqrt_newton_aux(x1,y1,eps) >= sqrt_newton_aux(x2,y2,eps)

  sqrt_newton(x,eps) : {z | x < sq(z) AND sq(z)-x < eps} =
    sqrt_newton_aux(x,x+1,eps)

  sqrt_newton_is_NOT_increasing: LEMMA
    sqrt_newton(1.61,0.99) < sqrt_newton(1.6,0.99)

  sq_sqrt_newton         : LEMMA sq(sqrt_newton(x,eps)) > x

  sqrt_newton_gt         : LEMMA sqrt_newton(x,eps) > sqrt(x)

  sqrt_newton_lt         : LEMMA x /= 0 IMPLIES x/sqrt_newton(x,eps) < sqrt(x)

  sqrt_newton_le         : LEMMA x/sqrt_newton(x,eps) <= sqrt(x)

  sqrt_approx(x,n) : {z | abs(sq(z)-x) < 10^(-n)/2} =
    LET eps = 10^(-n),
        N = 2*10^n,
k = floor(x*N),
xx = k/N
    IN sqrt_newton(xx,10^(-n)/2)

  sqrt_approx_increasing: LEMMA x<=y IMPLIES sqrt_approx(x,n)<=sqrt_approx(y,n)

  sqrt_ub_old(x,n)    : posreal = sqrt_approx(x,2*n)+10^(-n)

  sqrt_lb_old(x,n)    : nnreal = max(sqrt_approx(x,2*n)-10^(-n),0)

  sqrt_old_bounds      : LEMMA
sqrt_lb_old(x,n) <= sqrt(x) AND
        sqrt(x) < sqrt_ub_old(x,n)

  sqrt_approx_fast(x,n) : {z | abs(sq(z)-x) < 10^(-n)/2} =
    IF x<=64 THEN LET sa = sqrt_approx(x,n) IN min(sa,8) ELSE
      LET twoexp = log_nat(x,4)`1,
          newx = x/(4^twoexp),
   exp2err = floor(twoexp/3)+1,
   newn = n + 2*exp2err
      IN LET sa = (2^twoexp)*min(sqrt_approx(newx,newn),2) IN max(sa,8)
    ENDIF

  sqrt_approx_fast_increasing: LEMMA
      x<=y IMPLIES sqrt_approx_fast(x,n)<=sqrt_approx_fast(y,n)

  sqrt_ub(x,n)    : posreal = sqrt_approx_fast(x,2*n)+10^(-n)

  sqrt_lb(x,n)    : nnreal = max(sqrt_approx_fast(x,2*n)-10^(-n),0)

  sqrt_bounds      : LEMMA
sqrt_lb(x,n) <= sqrt(x) AND
        sqrt(x) < sqrt_ub(x,n)

end sqrt_approx

quality97%

¤ Diese beiden folgenden Angebotsgruppen bietet das Unternehmen0.1Angebot Wie Sie bei der Firma Beratungs- und Dienstleistungen beauftragen können ¤

*Eine klare Vorstellung vom Zielzustand

Wurzel

Suchen

Beweissystem der NASA

Beweissystem Isabelle

NIST Cobol Testsuite

Cephes Mathematical Library

Wiener Entwicklungsmethode

Haftungshinweis

Die Informationen auf dieser Webseite wurden nach bestem Wissen sorgfältig zusammengestellt. Es wird jedoch weder Vollständigkeit, noch Richtigkeit, noch Qualität der bereit gestellten Informationen zugesichert.