diffsl 0.11.9

A compiler for a domain-specific language for ordinary differential equations (ODE).
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42

// RUN: %eopt --enzyme %s | FileCheck %s

module {
  func.func @square(%x : f64) -> (f64, f64) {
    %y = arith.mulf %x, %x : f64
    return %y, %x : f64, f64
  }
  func.func @dsq(%x : f64, %dx : f64) -> (f64, f64) {
    %r:2 = enzyme.fwddiff @square(%x, %dx) { activity=[#enzyme<activity enzyme_dup>], ret_activity=[#enzyme<activity enzyme_dupnoneed>, #enzyme<activity enzyme_dupnoneed>] } : (f64, f64) -> (f64, f64)
    return %r#0, %r#1 : f64, f64
  }
  func.func @dsq2(%x : f64, %dx : f64) -> (f64, f64, f64, f64) {
    %r:4 = enzyme.fwddiff @square(%x, %dx) { activity=[#enzyme<activity enzyme_dup>], ret_activity=[#enzyme<activity enzyme_dup>, #enzyme<activity enzyme_dup>] } : (f64, f64) -> (f64, f64, f64, f64)
    return %r#0, %r#1, %r#2, %r#3 : f64, f64, f64, f64
  }
}

// CHECK:   func.func @dsq(%[[arg0:.+]]: f64, %[[arg1:.+]]: f64) -> (f64, f64) {
// CHECK-NEXT:     %[[i0:[0-9]+]]:2 = call @fwddiffesquare(%[[arg0]], %[[arg1]]) : (f64, f64) -> (f64, f64)
// CHECK-NEXT:     return %[[i0]]#0, %[[i0]]#1 : f64, f64
// CHECK-NEXT:   }

// CHECK:   func.func @dsq2(%[[arg0:.+]]: f64, %[[arg1:.+]]: f64) -> (f64, f64, f64, f64) {
// CHECK-NEXT:     %[[i0:[0-9]+]]:4 = call @fwddiffesquare_0(%[[arg0]], %[[arg1]]) : (f64, f64) -> (f64, f64, f64, f64)
// CHECK-NEXT:     return %[[i0]]#0, %[[i0]]#1, %[[i0]]#2, %[[i0]]#3 : f64, f64, f64, f64
// CHECK-NEXT:   }

// CHECK:   func.func private @fwddiffesquare(%[[arg0:.+]]: f64, %[[arg1:.+]]: f64) -> (f64, f64) {
// CHECK-NEXT:     %[[i0:.+]] = arith.mulf %[[arg1]], %[[arg0]] : f64
// CHECK-NEXT:     %[[i1:.+]] = arith.mulf %[[arg1]], %[[arg0]] : f64
// CHECK-NEXT:     %[[i2:.+]] = arith.addf %[[i0]], %[[i1]] : f64
// CHECK-NEXT:     %[[i3:.+]] = arith.mulf %[[arg0]], %[[arg0]] : f64
// CHECK-NEXT:     return %[[i2]], %[[arg1]] : f64, f64
// CHECK-NEXT:   }

// CHECK:   func.func private @fwddiffesquare_0(%[[arg0:.+]]: f64, %[[arg1:.+]]: f64) -> (f64, f64, f64, f64) {
// CHECK-NEXT:     %[[i0:.+]] = arith.mulf %[[arg1]], %[[arg0]] : f64
// CHECK-NEXT:     %[[i1:.+]] = arith.mulf %[[arg1]], %[[arg0]] : f64
// CHECK-NEXT:     %[[i2:.+]] = arith.addf %[[i0]], %[[i1]] : f64
// CHECK-NEXT:     %[[i3:.+]] = arith.mulf %[[arg0]], %[[arg0]] : f64
// CHECK-NEXT:     return %[[i3]], %[[i2]], %[[arg0]], %[[arg1]] : f64, f64, f64, f64
// CHECK-NEXT:   }