Commit 43021c02 authored by Anastasiia's avatar Anastasiia

new input intervals to avoid divbyzero for some benchmarks

parent 7fbe5cbb
......@@ -6,27 +6,6 @@ import Real._
object Turbine {
// def turbine1(v: Real): Real = {
// require(-4.5 <= v && v <= -0.3)
//
// 3 + 2/(v*v) - 0.125*(3-2*v)*(v*v*v*v)/(1-v) - 4.5
//
// }
//
// def turbine2(v: Real, w: Real, r: Real): Real = {
// require(-4.5 <= v && v <= -0.3 && 0.4 <= w && w <= 0.9 && 3.8 <= r && r <= 7.8)
//
// 6*v - 0.5 * v * (v*v*v*v) / (1-v) - 2.5
//
// }
//
// def turbine3(v: Real, w: Real, r: Real): Real = {
// require(-4.5 <= v && v <= -0.3 && 0.4 <= w && w <= 0.9 && 3.8 <= r && r <= 7.8)
//
// 3 - 2/(v*v) - 0.125 * (1+2*v) * (v*v*v*v) / (1-v) - 0.5
//
// }
def turbine1(v: Real, w: Real, r: Real): Real = {
require(-4.5 <= v && v <= -0.3 && 0.4 <= w && w <= 0.9 && 3.8 <= r && r <= 7.8)
......
import daisy.lang._
import Real._
object Sqrt {
def sqroot(x: Real): Real = {
require(x >= 0.0 && x < 1) // FPTaylor uses [ 0; 1 ]
1.0 + 0.5*x - 0.125*x*x + 0.0625*x*x*x - 0.0390625*x*x*x*x
}
}
\ No newline at end of file
import daisy.lang._
import Real._
object Turbine {
def turbine1(v: Real, w: Real, r: Real): Real = {
require(-4.5 <= v && v <= -0.3 && 0.4 <= w && w <= 0.9 && 3.8 <= r && r <= 7.8)
3 + 2/(r*r) - 0.125*(3-2*v)*(w*w*r*r)/(1-v) - 4.5
}
def turbine2(v: Real, w: Real, r: Real): Real = {
require(-4.5 <= v && v <= -3.3 && -0.4 <= w && w <= -0.1 && 3.8 <= r && r <= 7.8)
6*v - 0.5 * v * (w*w*r*r) / (1-v) - 2.5
}
def turbine3(v: Real, w: Real, r: Real): Real = {
require(-4.5 <= v && v <= -0.3 && 0.4 <= w && w <= 0.9 && 3.8 <= r && r <= 7.8)
3 - 2/(r*r) - 0.125 * (1+2*v) * (w*w*r*r) / (1-v) - 0.5
}
}
\ No newline at end of file
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment