Skip to content
Snippets Groups Projects
Commit bcd8aeaa authored by Krueger Jasmin's avatar Krueger Jasmin
Browse files

Examples Hensel validatin a posteriori

parent 6fcbf37e
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
Showing
with 63 additions and 9 deletions
examplesHensel/example1.jl
+ 4
− 1
View file @ bcd8aeaa
......@@ -77,6 +77,7 @@ hensel_test(fr, gr, hr, fq, gq, hq, 7, R)
####################################
#hensel in intervals and validation#
####################################
R = 0.3999999999068677425384521484375
gq = y_q -(2*x_q^3 -5*x_q^2 + 8*x_q -5)
hq = y_q - (x_q^3 -3* x_q^2 + 3*x_q -1)
fq = y_q^2 + 6*y_q +5 -10*x_q
......@@ -113,3 +114,5 @@ error_measure1(gistar, R/2)
error_measure1(histar, R/2)
error_measure1(gvstar, R/2)
error_measure1(hvstar, R/2)
error_measure1(grv, R/2)
error_measure1(hrv, R/2)
\ No newline at end of file
examplesHensel/example10.jl
+ 7
− 2
View file @ bcd8aeaa
......@@ -52,11 +52,12 @@ hr = convert_polynomial(hq,[RXY,RX],RDF)
fi = convert_polynomial(fq, [RBXY,RBX], RR)
fr = convert_polynomial(fq,[RXY,RX],RDF)
deg = 128
deg = 64#128
@time (gqstar, hqstar) = henseltruncate(deg,fq,gq,hq)
@time (gistar, histar) = henseltruncate(deg,fi,gi,hi)
@time (gvstar, hvstar) = henseltruncate(deg,fi,gi,hi,validation = true)#@time (grstar, hrstar) = henseltruncate(128,fr,gr,hr)
@time (grstar, hrstar) = henseltruncate(deg,fr,gr,hr)
@time (grv,hrv) = hensel_validate(deg,fi,gi,hi)
norm1error(grstar,gqstar)
norm1error(hrstar,hqstar)
......@@ -64,6 +65,8 @@ norm1error(gistar)
norm1error(histar)
norm1error(gvstar)
norm1error(hvstar)
norm1error(grv)
norm1error(hrv)
range_of_coefficients(fq)
......@@ -73,3 +76,5 @@ error_measure1(gistar, R/2)
error_measure1(histar, R/2)
error_measure1(gvstar, R/2)
error_measure1(hvstar, R/2)
error_measure1(grv, R/2)
error_measure1(hrv, R/2)
\ No newline at end of file
examplesHensel/example11.jl
+ 5
− 0
View file @ bcd8aeaa
......@@ -53,6 +53,7 @@ deg = 128
@time (gistar, histar) = henseltruncate(deg,fi,gi,hi)
@time (gvstar, hvstar) = henseltruncate(deg,fi,gi,hi,validation = true)#@time (grstar, hrstar) = henseltruncate(128,fr,gr,hr)
@time (grstar, hrstar) = henseltruncate(deg,fr,gr,hr)
@time (grv,hrv) = hensel_validate(deg,fi,gi,hi)
norm1error(grstar,gqstar)
norm1error(hrstar,hqstar)
......@@ -60,6 +61,8 @@ norm1error(gistar)
norm1error(histar)
norm1error(gvstar)
norm1error(hvstar)
norm1error(grv)
norm1error(hrv)
range_of_coefficients(fq)
......@@ -70,3 +73,5 @@ error_measure1(gistar, R/2)
error_measure1(histar, R/2)
error_measure1(gvstar, R/2)
error_measure1(hvstar, R/2)
error_measure1(grv, R/2)
error_measure1(hrv, R/2)
\ No newline at end of file
examplesHensel/example12.jl
+ 6
− 1
View file @ bcd8aeaa
......@@ -57,6 +57,7 @@ deg = 128
@time (gistar, histar) = henseltruncate(deg,fi,gi,hi)
@time (gvstar, hvstar) = henseltruncate(deg,fi,gi,hi,validation = true)#@time (grstar, hrstar) = henseltruncate(128,fr,gr,hr)
@time (grstar, hrstar) = henseltruncate(deg,fr,gr,hr)
@time (grv,hrv) = hensel_validate(deg,fi,gi,hi)
norm1error(grstar,gqstar)
norm1error(hrstar,hqstar)
......@@ -64,6 +65,8 @@ norm1error(gistar)
norm1error(histar)
norm1error(gvstar)
norm1error(hvstar)
norm1error(grv)
norm1error(hrv)
range_of_coefficients(fq)
......@@ -73,3 +76,5 @@ error_measure1(gistar, R/2)
error_measure1(histar, R/2)
error_measure1(gvstar, R/2)
error_measure1(hvstar, R/2)
error_measure1(grv, R/2)
error_measure1(hrv, R/2)
\ No newline at end of file
examplesHensel/example13.jl
+ 5
− 0
View file @ bcd8aeaa
......@@ -54,6 +54,7 @@ deg = 128
@time (gistar, histar) = henseltruncate(deg,fi,gi,hi)
@time (gvstar, hvstar) = henseltruncate(deg,fi,gi,hi,validation = true)#@time (grstar, hrstar) = henseltruncate(128,fr,gr,hr)
@time (grstar, hrstar) = henseltruncate(deg,fr,gr,hr)
@time (grv,hrv) = hensel_validate(deg,fi,gi,hi);
norm1error(grstar,gqstar)
norm1error(hrstar,hqstar)
......@@ -61,6 +62,8 @@ norm1error(gistar)
norm1error(histar)
norm1error(gvstar)
norm1error(hvstar)
norm1error(grv)
norm1error(hrv)
range_of_coefficients(fq)
......@@ -70,3 +73,5 @@ error_measure1(gistar, R/2)
error_measure1(histar, R/2)
error_measure1(gvstar, R/2)
error_measure1(hvstar, R/2)
error_measure1(grv, R/2)
error_measure1(hrv, R/2)
\ No newline at end of file
examplesHensel/example3.jl
+ 3
− 0
View file @ bcd8aeaa
......@@ -104,6 +104,7 @@ deg = 32
@time (gistar, histar) = henseltruncate(deg,fi,gi,hi)
@time (gvstar, hvstar) = henseltruncate(deg,fi,gi,hi,validation = true)#@time (grstar, hrstar) = henseltruncate(128,fr,gr,hr)
@time (grstar, hrstar) = henseltruncate(deg,fr,gr,hr)
@time (grv,hrv) = hensel_validate(deg,fi,gi,hi, imax = 5);
norm1error(grstar,gqstar)
norm1error(hrstar,hqstar)
......@@ -111,6 +112,8 @@ norm1error(gistar)
norm1error(histar)
norm1error(gvstar)
norm1error(hvstar)
norm1error(grv)
norm1error(hrv)
range_of_coefficients(Fq)
......
examplesHensel/example4.jl
+ 3
− 0
View file @ bcd8aeaa
......@@ -97,6 +97,7 @@ deg = 128
@time (gistar, histar) = henseltruncate(deg,fi,gi,hi)
@time (gvstar, hvstar) = henseltruncate(deg,fi,gi,hi,validation = true)#@time (grstar, hrstar) = henseltruncate(128,fr,gr,hr)
@time (grstar, hrstar) = henseltruncate(deg,fr,gr,hr)
@time (grv,hrv) = hensel_validate(deg,fi,gi,hi);
norm1error(grstar,gqstar)
norm1error(hrstar,hqstar)
......@@ -104,6 +105,8 @@ norm1error(gistar)
norm1error(histar)
norm1error(gvstar)
norm1error(hvstar)
norm1error(grv)
norm1error(hrv)
range_of_coefficients(fq)
range_of_coefficients(gq)
......
examplesHensel/example5.jl
+ 6
− 1
View file @ bcd8aeaa
......@@ -93,6 +93,7 @@ deg = 128
@time (gistar, histar) = henseltruncate(deg,fi,gi,hi)
@time (gvstar, hvstar) = henseltruncate(deg,fi,gi,hi,validation = true)#@time (grstar, hrstar) = henseltruncate(128,fr,gr,hr)
@time (grstar, hrstar) = henseltruncate(deg,fr,gr,hr)
@time (grv,hrv) = hensel_validate(deg,fi,gi,hi);
norm1error(grstar,gqstar)
norm1error(hrstar,hqstar)
......@@ -100,6 +101,8 @@ norm1error(gistar)
norm1error(histar)
norm1error(gvstar)
norm1error(hvstar)
norm1error(grv)
norm1error(hrv)
range_of_coefficients(fq)
range_of_coefficients(gq)
......@@ -111,3 +114,5 @@ error_measure1(gistar, R/2)
error_measure1(histar, R/2)
error_measure1(gvstar, R/2)
error_measure1(hvstar, R/2)
error_measure1(grv, R/2)
error_measure1(hrv, R/2)
\ No newline at end of file
examplesHensel/example6.jl
+ 6
− 1
View file @ bcd8aeaa
......@@ -97,6 +97,7 @@ deg = 128
@time (gistar, histar) = henseltruncate(deg,fi,gi,hi)
@time (gvstar, hvstar) = henseltruncate(deg,fi,gi,hi,validation = true)#@time (grstar, hrstar) = henseltruncate(128,fr,gr,hr)
@time (grstar, hrstar) = henseltruncate(deg,fr,gr,hr)
@time (grv,hrv) = hensel_validate(deg,fi,gi,hi);
norm1error(grstar,gqstar)
norm1error(hrstar,hqstar)
......@@ -104,6 +105,8 @@ norm1error(gistar)
norm1error(histar)
norm1error(gvstar)
norm1error(hvstar)
norm1error(grv)
norm1error(hrv)
range_of_coefficients(fq)
range_of_coefficients(gq)
......@@ -115,3 +118,5 @@ error_measure1(gistar, R/2)
error_measure1(histar, R/2)
error_measure1(gvstar, R/2)
error_measure1(hvstar, R/2)
error_measure1(grv, R/2)
error_measure1(hrv, R/2)
examplesHensel/example7.jl
+ 6
− 1
View file @ bcd8aeaa
......@@ -68,6 +68,7 @@ deg = 128
@time (gistar, histar) = henseltruncate(deg,fi,gi,hi)
@time (gvstar, hvstar) = henseltruncate(deg,fi,gi,hi,validation = true)#@time (grstar, hrstar) = henseltruncate(128,fr,gr,hr)
@time (grstar, hrstar) = henseltruncate(deg,fr,gr,hr)
@time (grv,hrv) = hensel_validate(deg,fi,gi,hi);
norm1error(grstar,gqstar)
norm1error(hrstar,hqstar)
......@@ -75,6 +76,8 @@ norm1error(gistar)
norm1error(histar)
norm1error(gvstar)
norm1error(hvstar)
norm1error(grv)
norm1error(hrv)
range_of_coefficients(fq)
......@@ -84,3 +87,5 @@ error_measure1(gistar, R/2)
error_measure1(histar, R/2)
error_measure1(gvstar, R/2)
error_measure1(hvstar, R/2)
error_measure1(grv, R/2)
error_measure1(hrv, R/2)
\ No newline at end of file
examplesHensel/example8.jl
+ 6
− 1
View file @ bcd8aeaa
......@@ -100,6 +100,7 @@ deg = 128
@time (gistar, histar) = henseltruncate(deg,fi,gi,hi)
@time (gvstar, hvstar) = henseltruncate(deg,fi,gi,hi,validation = true)#@time (grstar, hrstar) = henseltruncate(128,fr,gr,hr)
@time (grstar, hrstar) = henseltruncate(deg,fr,gr,hr)
@time (grv,hrv) = hensel_validate(deg,fi,gi,hi);
norm1error(grstar,gqstar)
norm1error(hrstar,hqstar)
......@@ -107,6 +108,8 @@ norm1error(gistar)
norm1error(histar)
norm1error(gvstar)
norm1error(hvstar)
norm1error(grv)
norm1error(hrv)
range_of_coefficients(fq)
......@@ -116,3 +119,5 @@ error_measure1(gistar, R/2)
error_measure1(histar, R/2)
error_measure1(gvstar, R/2)
error_measure1(hvstar, R/2)
error_measure1(grv, R/2)
error_measure1(hrv, R/2)
\ No newline at end of file
examplesHensel/example9.jl
+ 6
− 1
View file @ bcd8aeaa
......@@ -59,6 +59,7 @@ deg = 128
@time (gistar, histar) = henseltruncate(deg,fi,gi,hi)
@time (gvstar, hvstar) = henseltruncate(deg,fi,gi,hi,validation = true)#@time (grstar, hrstar) = henseltruncate(128,fr,gr,hr)
@time (grstar, hrstar) = henseltruncate(deg,fr,gr,hr)
@time (grv,hrv) = hensel_validate(deg,fi,gi,hi)
norm1error(grstar,gqstar)
norm1error(hrstar,hqstar)
......@@ -66,6 +67,8 @@ norm1error(gistar)
norm1error(histar)
norm1error(gvstar)
norm1error(hvstar)
norm1error(grv)
norm1error(hrv)
range_of_coefficients(fq)
......@@ -75,3 +78,5 @@ error_measure1(gistar, R/2)
error_measure1(histar, R/2)
error_measure1(gvstar, R/2)
error_measure1(hvstar, R/2)
error_measure1(grv, R/2)
error_measure1(hrv, R/2)
\ No newline at end of file
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment