kscience/kmath
9
0
Fork 1
Code Issues Pull Requests 6 Packages Projects Releases 16 Wiki Activity

Gauss-Legendre working test

Browse Source
This commit is contained in:
Alexander Nozik 2021-04-15 19:39:46 +03:00
parent 93bc371622
commit 19ec6a57a4
2 changed files with 12 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

18
kmath-functions/src/commonMain/kotlin/space/kscience/kmath/integration/GaussIntegratorRuleFactory.kt
View File

@ -2,7 +2,6 @@ package space.kscience.kmath.integration
import space.kscience.kmath.operations.DoubleField import space.kscience.kmath.operations.DoubleField
import space.kscience.kmath.operations.Field import space.kscience.kmath.operations.Field
import space.kscience.kmath.operations.Ring
import space.kscience.kmath.structures.* import space.kscience.kmath.structures.*
import kotlin.jvm.Synchronized import kotlin.jvm.Synchronized
import kotlin.math.ulp import kotlin.math.ulp
@ -27,14 +26,19 @@ public fun <T : Comparable<T>> GaussIntegratorRuleFactory<T>.build(
range: ClosedRange<T>, range: ClosedRange<T>,
): Pair<Buffer<T>, Buffer<T>> { ): Pair<Buffer<T>, Buffer<T>> {
val normalized = build(numPoints) val normalized = build(numPoints)
val points = with(algebra) { with(algebra) {
val length = range.endInclusive - range.start val length = range.endInclusive - range.start
normalized.first.map(bufferFactory) {
range.start + length / 2 + length * it/2
}
}
return points to normalized.second val points = normalized.first.map(bufferFactory) {
range.start + length / 2 + length * it / 2
}
val weights = normalized.second.map(bufferFactory) {
it * length / 2
}
return points to weights
}
} }

2
kmath-functions/src/commonTest/kotlin/space/kscience/kmath/integration/GaussIntegralTest.kt
View File

@ -23,7 +23,7 @@ class GaussIntegralTest {
0.0 0.0
} }
} }
assertEquals(20.0, res.value!!, 0.1) assertEquals(20.0, res.value!!, 0.5)
} }