2021-03-22 14:32:08 +03:00
|
|
|
# Module kmath-ast
|
2020-12-22 14:02:02 +03:00
|
|
|
|
2021-05-03 00:14:19 +03:00
|
|
|
Performance and visualization extensions to MST API.
|
2020-12-22 14:02:02 +03:00
|
|
|
|
|
|
|
${features}
|
|
|
|
|
|
|
|
${artifact}
|
|
|
|
|
|
|
|
## Dynamic expression code generation
|
|
|
|
|
|
|
|
### On JVM
|
|
|
|
|
2021-05-03 00:14:19 +03:00
|
|
|
`kmath-ast` JVM module supports runtime code generation to eliminate overhead of tree traversal. Code generator builds a
|
|
|
|
special implementation of `Expression<T>` with implemented `invoke` function.
|
2020-12-22 14:02:02 +03:00
|
|
|
|
|
|
|
For example, the following builder:
|
|
|
|
|
|
|
|
```kotlin
|
2021-05-26 16:24:29 +03:00
|
|
|
import space.kscience.kmath.expressions.Symbol.Companion.x
|
2021-05-03 00:14:19 +03:00
|
|
|
import space.kscience.kmath.expressions.*
|
|
|
|
import space.kscience.kmath.operations.*
|
|
|
|
import space.kscience.kmath.asm.*
|
|
|
|
|
2021-05-26 16:24:29 +03:00
|
|
|
MstField { x + 2 }.compileToExpression(DoubleField)
|
2020-12-22 14:02:02 +03:00
|
|
|
```
|
|
|
|
|
2021-04-23 04:49:01 +03:00
|
|
|
... leads to generation of bytecode, which can be decompiled to the following Java class:
|
2020-12-22 14:02:02 +03:00
|
|
|
|
|
|
|
```java
|
2021-02-19 10:52:56 +03:00
|
|
|
package space.kscience.kmath.asm.generated;
|
2020-12-22 14:02:02 +03:00
|
|
|
|
|
|
|
import java.util.Map;
|
2021-05-03 00:14:19 +03:00
|
|
|
|
2020-12-22 14:02:02 +03:00
|
|
|
import kotlin.jvm.functions.Function2;
|
2021-02-18 11:17:28 +03:00
|
|
|
import space.kscience.kmath.asm.internal.MapIntrinsics;
|
|
|
|
import space.kscience.kmath.expressions.Expression;
|
|
|
|
import space.kscience.kmath.expressions.Symbol;
|
2020-12-22 14:02:02 +03:00
|
|
|
|
|
|
|
public final class AsmCompiledExpression_45045_0 implements Expression<Double> {
|
|
|
|
private final Object[] constants;
|
|
|
|
|
2021-02-11 19:35:45 +03:00
|
|
|
public final Double invoke(Map<Symbol, ? extends Double> arguments) {
|
2021-05-03 00:14:19 +03:00
|
|
|
return (Double) ((Function2) this.constants[0]).invoke((Double) MapIntrinsics.getOrFail(arguments, "x"), 2);
|
2020-12-22 14:02:02 +03:00
|
|
|
}
|
|
|
|
|
|
|
|
public AsmCompiledExpression_45045_0(Object[] constants) {
|
|
|
|
this.constants = constants;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
```
|
|
|
|
|
|
|
|
#### Known issues
|
|
|
|
|
2021-05-03 00:14:19 +03:00
|
|
|
- The same classes may be generated and loaded twice, so it is recommended to cache compiled expressions to avoid class
|
|
|
|
loading overhead.
|
2020-12-22 14:02:02 +03:00
|
|
|
- This API is not supported by non-dynamic JVM implementations (like TeaVM and GraalVM) because of using class loaders.
|
|
|
|
|
|
|
|
### On JS
|
|
|
|
|
|
|
|
A similar feature is also available on JS.
|
|
|
|
|
|
|
|
```kotlin
|
2021-05-26 16:24:29 +03:00
|
|
|
import space.kscience.kmath.expressions.Symbol.Companion.x
|
2021-05-03 00:14:19 +03:00
|
|
|
import space.kscience.kmath.expressions.*
|
|
|
|
import space.kscience.kmath.operations.*
|
|
|
|
import space.kscience.kmath.estree.*
|
|
|
|
|
2021-05-26 16:24:29 +03:00
|
|
|
MstField { x + 2 }.compileToExpression(DoubleField)
|
2020-12-22 14:02:02 +03:00
|
|
|
```
|
|
|
|
|
|
|
|
The code above returns expression implemented with such a JS function:
|
|
|
|
|
|
|
|
```js
|
|
|
|
var executable = function (constants, arguments) {
|
2021-05-03 00:14:19 +03:00
|
|
|
return constants[1](constants[0](arguments, "x"), 2);
|
2020-12-22 14:02:02 +03:00
|
|
|
};
|
|
|
|
```
|
|
|
|
|
2021-05-07 15:59:21 +03:00
|
|
|
JS also supports experimental expression optimization with [WebAssembly](https://webassembly.org/) IR generation.
|
2021-05-03 00:14:19 +03:00
|
|
|
Currently, only expressions inside `DoubleField` and `IntRing` are supported.
|
2021-04-20 17:39:45 +03:00
|
|
|
|
|
|
|
```kotlin
|
2021-05-26 16:24:29 +03:00
|
|
|
import space.kscience.kmath.expressions.Symbol.Companion.x
|
2021-05-03 00:14:19 +03:00
|
|
|
import space.kscience.kmath.expressions.*
|
|
|
|
import space.kscience.kmath.operations.*
|
2021-04-20 17:39:45 +03:00
|
|
|
import space.kscience.kmath.wasm.*
|
|
|
|
|
2021-05-26 16:24:29 +03:00
|
|
|
MstField { x + 2 }.compileToExpression(DoubleField)
|
2021-04-20 17:39:45 +03:00
|
|
|
```
|
|
|
|
|
2021-04-23 04:49:01 +03:00
|
|
|
An example of emitted Wasm IR in the form of WAT:
|
2021-04-20 17:39:45 +03:00
|
|
|
|
|
|
|
```lisp
|
|
|
|
(func \$executable (param \$0 f64) (result f64)
|
|
|
|
(f64.add
|
|
|
|
(local.get \$0)
|
|
|
|
(f64.const 2)
|
|
|
|
)
|
|
|
|
)
|
|
|
|
```
|
|
|
|
|
2020-12-22 14:02:02 +03:00
|
|
|
#### Known issues
|
|
|
|
|
2021-04-20 17:39:45 +03:00
|
|
|
- ESTree expression compilation uses `eval` which can be unavailable in several environments.
|
|
|
|
- WebAssembly isn't supported by old versions of browsers (see https://webassembly.org/roadmap/).
|
2021-03-25 19:57:47 +03:00
|
|
|
|
|
|
|
## Rendering expressions
|
|
|
|
|
2021-04-20 17:39:45 +03:00
|
|
|
kmath-ast also includes an extensible engine to display expressions in LaTeX or MathML syntax.
|
2021-03-25 19:57:47 +03:00
|
|
|
|
|
|
|
Example usage:
|
|
|
|
|
|
|
|
```kotlin
|
|
|
|
import space.kscience.kmath.ast.*
|
|
|
|
import space.kscience.kmath.ast.rendering.*
|
2021-05-03 00:14:19 +03:00
|
|
|
import space.kscience.kmath.misc.*
|
2021-03-25 19:57:47 +03:00
|
|
|
|
2021-05-03 00:14:19 +03:00
|
|
|
@OptIn(UnstableKMathAPI::class)
|
2021-03-25 19:57:47 +03:00
|
|
|
public fun main() {
|
2021-05-03 00:14:19 +03:00
|
|
|
val mst = "exp(sqrt(x))-asin(2*x)/(2e10+x^3)/(12)+x^(2/3)".parseMath()
|
2021-03-25 19:57:47 +03:00
|
|
|
val syntax = FeaturedMathRendererWithPostProcess.Default.render(mst)
|
|
|
|
val latex = LatexSyntaxRenderer.renderWithStringBuilder(syntax)
|
|
|
|
println("LaTeX:")
|
|
|
|
println(latex)
|
|
|
|
println()
|
|
|
|
val mathML = MathMLSyntaxRenderer.renderWithStringBuilder(syntax)
|
|
|
|
println("MathML:")
|
|
|
|
println(mathML)
|
|
|
|
}
|
|
|
|
```
|
|
|
|
|
2021-04-20 17:39:45 +03:00
|
|
|
Result LaTeX:
|
2021-03-25 19:57:47 +03:00
|
|
|
|
2021-05-07 15:59:21 +03:00
|
|
|
<div style="background-color:white;">
|
|
|
|
|
2021-05-03 00:14:19 +03:00
|
|
|
![](https://latex.codecogs.com/gif.latex?%5Coperatorname{exp}%5C,%5Cleft(%5Csqrt{x}%5Cright)-%5Cfrac{%5Cfrac{%5Coperatorname{arcsin}%5C,%5Cleft(2%5C,x%5Cright)}{2%5Ctimes10^{10}%2Bx^{3}}}{12}+x^{2/3})
|
2021-05-07 15:59:21 +03:00
|
|
|
</div>
|
2021-03-25 19:57:47 +03:00
|
|
|
|
2021-05-03 00:14:19 +03:00
|
|
|
Result MathML (can be used with MathJax or other renderers):
|
|
|
|
|
|
|
|
<details>
|
2021-03-25 19:57:47 +03:00
|
|
|
|
|
|
|
```html
|
2021-05-03 00:14:19 +03:00
|
|
|
<math xmlns="https://www.w3.org/1998/Math/MathML">
|
|
|
|
<mrow>
|
|
|
|
<mo>exp</mo>
|
|
|
|
<mspace width="0.167em"></mspace>
|
|
|
|
<mfenced open="(" close=")" separators="">
|
|
|
|
<msqrt>
|
|
|
|
<mi>x</mi>
|
|
|
|
</msqrt>
|
|
|
|
</mfenced>
|
|
|
|
<mo>-</mo>
|
|
|
|
<mfrac>
|
|
|
|
<mrow>
|
|
|
|
<mfrac>
|
|
|
|
<mrow>
|
|
|
|
<mo>arcsin</mo>
|
|
|
|
<mspace width="0.167em"></mspace>
|
|
|
|
<mfenced open="(" close=")" separators="">
|
|
|
|
<mn>2</mn>
|
|
|
|
<mspace width="0.167em"></mspace>
|
|
|
|
<mi>x</mi>
|
|
|
|
</mfenced>
|
|
|
|
</mrow>
|
|
|
|
<mrow>
|
|
|
|
<mn>2</mn>
|
|
|
|
<mo>×</mo>
|
|
|
|
<msup>
|
|
|
|
<mrow>
|
|
|
|
<mn>10</mn>
|
|
|
|
</mrow>
|
|
|
|
<mrow>
|
|
|
|
<mn>10</mn>
|
|
|
|
</mrow>
|
|
|
|
</msup>
|
|
|
|
<mo>+</mo>
|
|
|
|
<msup>
|
|
|
|
<mrow>
|
|
|
|
<mi>x</mi>
|
|
|
|
</mrow>
|
|
|
|
<mrow>
|
|
|
|
<mn>3</mn>
|
|
|
|
</mrow>
|
|
|
|
</msup>
|
|
|
|
</mrow>
|
|
|
|
</mfrac>
|
|
|
|
</mrow>
|
|
|
|
<mrow>
|
|
|
|
<mn>12</mn>
|
|
|
|
</mrow>
|
|
|
|
</mfrac>
|
|
|
|
<mo>+</mo>
|
|
|
|
<msup>
|
|
|
|
<mrow>
|
|
|
|
<mi>x</mi>
|
|
|
|
</mrow>
|
|
|
|
<mrow>
|
|
|
|
<mn>2</mn>
|
|
|
|
<mo>/</mo>
|
|
|
|
<mn>3</mn>
|
|
|
|
</mrow>
|
|
|
|
</msup>
|
|
|
|
</mrow>
|
|
|
|
</math>
|
2021-03-25 19:57:47 +03:00
|
|
|
```
|
|
|
|
|
2021-05-03 00:14:19 +03:00
|
|
|
</details>
|
|
|
|
|
2021-04-20 17:39:45 +03:00
|
|
|
It is also possible to create custom algorithms of render, and even add support of other markup languages
|
2021-03-25 19:57:47 +03:00
|
|
|
(see API reference).
|