user: 下面的汇编指令对应的C源代码是什么:

```asm
<Init_Qmat_WAG>:
  endbr64
  movapd D693(%rip),%xmm0
  mov    %rdi,%rcx
  lea    0x8(%rdi),%r8
  mov    D100(%rip),%rax
  lea    0xa0(%rdi),%r9
  add    $0xa8,%rcx
  movups %xmm0,0x140(%rdi)
  movapd D863(%rip),%xmm0
  mov    %rax,0xa0(%rdi)
  mov    D621(%rip),%rax
  movups %xmm0,0x1e0(%rdi)
  movapd D558(%rip),%xmm0
  mov    %rax,0x1f0(%rdi)
  mov    D954(%rip),%rax
  movups %xmm0,0x280(%rdi)
  movapd D905(%rip),%xmm0
  mov    %rax,0x340(%rdi)
  mov    D1071(%rip),%rax
  movups %xmm0,0x290(%rdi)
  movapd D706(%rip),%xmm0
  mov    %rax,0x490(%rdi)
  movups %xmm0,0x320(%rdi)
  movapd D246(%rip),%xmm0
  movups %xmm0,0x330(%rdi)
  movapd D794(%rip),%xmm0
  movups %xmm0,0x3c0(%rdi)
  movapd D557(%rip),%xmm0
  movups %xmm0,0x3d0(%rdi)
  movapd D190(%rip),%xmm0
  movups %xmm0,0x3e0(%rdi)
  movapd D797(%rip),%xmm0
  movups %xmm0,0x460(%rdi)
  movapd D480(%rip),%xmm0
  movups %xmm0,0x470(%rdi)
  movapd D4(%rip),%xmm0
  movups %xmm0,0x480(%rdi)
  movapd D637(%rip),%xmm0
  mov    D164(%rip),%rax
  movups %xmm0,0x500(%rdi)
  movapd D695(%rip),%xmm0
  mov    %rax,0x5e0(%rdi)
  movups %xmm0,0x510(%rdi)
  movapd D38(%rip),%xmm0
  movups %xmm0,0x520(%rdi)
  movapd D945(%rip),%xmm0
  movups %xmm0,0x530(%rdi)
  movapd D729(%rip),%xmm0
  movups %xmm0,0x5a0(%rdi)
  movapd D737(%rip),%xmm0
  movups %xmm0,0x5b0(%rdi)
  movapd D177(%rip),%xmm0
  movups %xmm0,0x5c0(%rdi)
  movapd D800(%rip),%xmm0
  movups %xmm0,0x5d0(%rdi)
  movapd D518(%rip),%xmm0
  movups %xmm0,0x640(%rdi)
  movapd D24(%rip),%xmm0
  movups %xmm0,0x650(%rdi)
  movapd D421(%rip),%xmm0
  movups %xmm0,0x660(%rdi)
  movapd D782(%rip),%xmm0
  movups %xmm0,0x670(%rdi)
  movapd D790(%rip),%xmm0
  movups %xmm0,0x680(%rdi)
  movapd D172(%rip),%xmm0
  movups %xmm0,0x6e0(%rdi)
  movapd D968(%rip),%xmm0
  movups %xmm0,0x6f0(%rdi)
  movapd D799(%rip),%xmm0
  mov    D366(%rip),%rax
  movups %xmm0,0x700(%rdi)
  movapd D676(%rip),%xmm0
  mov    %rax,0x730(%rdi)
  movups %xmm0,0x710(%rdi)
  movapd D1035(%rip),%xmm0
  movups %xmm0,0x720(%rdi)
  movapd D269(%rip),%xmm0
  movups %xmm0,0x780(%rdi)
  movapd D1047(%rip),%xmm0
  movups %xmm0,0x790(%rdi)
  movapd D417(%rip),%xmm0
  movups %xmm0,0x7a0(%rdi)
  movapd D580(%rip),%xmm0
  movups %xmm0,0x7b0(%rdi)
  movapd D95(%rip),%xmm0
  movups %xmm0,0x7c0(%rdi)
  movapd D809(%rip),%xmm0
  movups %xmm0,0x7d0(%rdi)
  movapd D859(%rip),%xmm0
  movups %xmm0,0x820(%rdi)
  movapd D118(%rip),%xmm0
  movups %xmm0,0x830(%rdi)
  movapd D618(%rip),%xmm0
  movups %xmm0,0x840(%rdi)
  movapd D973(%rip),%xmm0
  movups %xmm0,0x850(%rdi)
  movapd D504(%rip),%xmm0
  movups %xmm0,0x860(%rdi)
  movapd D582(%rip),%xmm0
  movups %xmm0,0x870(%rdi)
  mov    D36(%rip),%rax
  movapd D895(%rip),%xmm0
  movups %xmm0,0x8c0(%rdi)
  movapd D856(%rip),%xmm0
  mov    %rax,0x880(%rdi)
  mov    D315(%rip),%rax
  movups %xmm0,0x8d0(%rdi)
  movapd D1049(%rip),%xmm0
  mov    %rax,0x9d0(%rdi)
  movups %xmm0,0x8e0(%rdi)
  movapd D354(%rip),%xmm0
  movups %xmm0,0x8f0(%rdi)
  movapd D84(%rip),%xmm0
  movups %xmm0,0x900(%rdi)
  movapd D699(%rip),%xmm0
  movups %xmm0,0x910(%rdi)
  movapd D892(%rip),%xmm0
  movups %xmm0,0x920(%rdi)
  movapd D910(%rip),%xmm0
  movups %xmm0,0x960(%rdi)
  movapd D1017(%rip),%xmm0
  movups %xmm0,0x970(%rdi)
  movapd D334(%rip),%xmm0
  movups %xmm0,0x980(%rdi)
  movapd D964(%rip),%xmm0
  movups %xmm0,0x990(%rdi)
  movapd D128(%rip),%xmm0
  movups %xmm0,0x9a0(%rdi)
  movapd D829(%rip),%xmm0
  movups %xmm0,0x9b0(%rdi)
  movapd D808(%rip),%xmm0
  movups %xmm0,0x9c0(%rdi)
  movapd D41(%rip),%xmm0
  movups %xmm0,0xa00(%rdi)
  movapd D14(%rip),%xmm0
  movups %xmm0,0xa10(%rdi)
  movapd D509(%rip),%xmm0
  movups %xmm0,0xa20(%rdi)
  movapd D312(%rip),%xmm0
  movups %xmm0,0xa30(%rdi)
  movapd D530(%rip),%xmm0
  movups %xmm0,0xa40(%rdi)
  movapd D862(%rip),%xmm0
  movups %xmm0,0xa50(%rdi)
  movapd D949(%rip),%xmm0
  movups %xmm0,0xa60(%rdi)
  movapd D891(%rip),%xmm0
  movups %xmm0,0xa70(%rdi)
  movapd D841(%rip),%xmm0
  movups %xmm0,0xaa0(%rdi)
  movapd D383(%rip),%xmm0
  movups %xmm0,0xab0(%rdi)
  movapd D115(%rip),%xmm0
  movups %xmm0,0xac0(%rdi)
  movapd D929(%rip),%xmm0
  movups %xmm0,0xad0(%rdi)
  movapd D830(%rip),%xmm0
  movups %xmm0,0xae0(%rdi)
  movapd D181(%rip),%xmm0
  movups %xmm0,0xaf0(%rdi)
  movapd D251(%rip),%xmm0
  movups %xmm0,0xb00(%rdi)
  movapd D355(%rip),%xmm0
  movups %xmm0,0xb10(%rdi)
  mov    D437(%rip),%rax
  movapd D391(%rip),%xmm0
  movups %xmm0,0xb40(%rdi)
  movapd D199(%rip),%xmm0
  mov    %rax,0xb20(%rdi)
  movups %xmm0,0xb50(%rdi)
  movapd D998(%rip),%xmm0
  movups %xmm0,0xb60(%rdi)
  movapd D74(%rip),%xmm0
  movups %xmm0,0xb70(%rdi)
  movapd D272(%rip),%xmm0
  movups %xmm0,0xb80(%rdi)
  movapd D771(%rip),%xmm0
  movups %xmm0,0xb90(%rdi)
  movapd D363(%rip),%xmm0
  movups %xmm0,0xba0(%rdi)
  movapd D754(%rip),%xmm0
  movups %xmm0,0xbb0(%rdi)
  movapd D53(%rip),%xmm0
  movups %xmm0,0xbc0(%rdi)
  movapd D293(%rip),%xmm0
  movups %xmm0,0xbe0(%rdi)
  movapd D171(%rip),%xmm0
  movups %xmm0,0xbf0(%rdi)
  movapd D529(%rip),%xmm0
  movups %xmm0,0xc00(%rdi)
  movapd D478(%rip),%xmm0
  movups %xmm0,0xc10(%rdi)
  movapd D68(%rip),%xmm0
  movups %xmm0,0xc20(%rdi)
  movapd D955(%rip),%xmm0
  movups %xmm0,0xc30(%rdi)
  movapd D789(%rip),%xmm0
  mov    D2(%rip),%rax
  movups %xmm0,0xc40(%rdi)
  movapd D375(%rip),%xmm0
  mov    %rax,0xc70(%rdi)
  movups %xmm0,0xc50(%rdi)
  movapd D757(%rip),%xmm0
  movups %xmm0,0xc60(%rdi)
  mov    $0x1,%edi
L16:
  mov    %r8,%rdx
  mov    %r9,%rax
  cs nopw 0x0(%rax,%rax,1)
L8:
  movsd  (%rax),%xmm0
  add    $0x8,%rax
  add    $0xa0,%rdx
  movsd  %xmm0,-0xa0(%rdx)
  cmp    %rax,%rcx
  jne    L8
  add    $0x1,%edi
  add    $0xa0,%r9
  add    $0x8,%r8
  add    $0xa8,%rcx
  cmp    $0x14,%edi
  jne    L16
  movapd D842(%rip),%xmm0
  mov    $0x1,%eax
  movups %xmm0,(%rsi)
  movapd D227(%rip),%xmm0
  movups %xmm0,0x10(%rsi)
  movapd D287(%rip),%xmm0
  movups %xmm0,0x20(%rsi)
  movapd D629(%rip),%xmm0
  movups %xmm0,0x30(%rsi)
  movapd D195(%rip),%xmm0
  movups %xmm0,0x40(%rsi)
  movapd D505(%rip),%xmm0
  movups %xmm0,0x50(%rsi)
  movapd D965(%rip),%xmm0
  movups %xmm0,0x60(%rsi)
  movapd D289(%rip),%xmm0
  movups %xmm0,0x70(%rsi)
  movapd D589(%rip),%xmm0
  movups %xmm0,0x80(%rsi)
  movapd D3(%rip),%xmm0
  movups %xmm0,0x90(%rsi)
  ret
```

tool_calls: name: read_data
arguments: data_label: D693
type_infered: f64[2]
name: read_data
arguments: data_label: D100
type_infered: f64
name: read_data
arguments: data_label: D863
type_infered: f64[2]
name: read_data
arguments: data_label: D621
type_infered: f64
name: read_data
arguments: data_label: D558
type_infered: f64[2]
name: read_data
arguments: data_label: D954
type_infered: f64
name: read_data
arguments: data_label: D905
type_infered: f64[2]
name: read_data
arguments: data_label: D1071
type_infered: f64
name: read_data
arguments: data_label: D706
type_infered: f64[2]
name: read_data
arguments: data_label: D246
type_infered: f64[2]
name: read_data
arguments: data_label: D794
type_infered: f64[2]
name: read_data
arguments: data_label: D557
type_infered: f64[2]
name: read_data
arguments: data_label: D190
type_infered: f64[2]
name: read_data
arguments: data_label: D797
type_infered: f64[2]
name: read_data
arguments: data_label: D480
type_infered: f64[2]
name: read_data
arguments: data_label: D4
type_infered: f64[2]
name: read_data
arguments: data_label: D637
type_infered: f64[2]
name: read_data
arguments: data_label: D164
type_infered: f64
name: read_data
arguments: data_label: D695
type_infered: f64[2]
name: read_data
arguments: data_label: D38
type_infered: f64[2]
name: read_data
arguments: data_label: D945
type_infered: f64[2]
name: read_data
arguments: data_label: D729
type_infered: f64[2]
name: read_data
arguments: data_label: D737
type_infered: f64[2]
name: read_data
arguments: data_label: D177
type_infered: f64[2]
name: read_data
arguments: data_label: D800
type_infered: f64[2]
name: read_data
arguments: data_label: D518
type_infered: f64[2]
name: read_data
arguments: data_label: D24
type_infered: f64[2]
name: read_data
arguments: data_label: D421
type_infered: f64[2]
name: read_data
arguments: data_label: D782
type_infered: f64[2]
name: read_data
arguments: data_label: D790
type_infered: f64[2]
name: read_data
arguments: data_label: D172
type_infered: f64[2]
name: read_data
arguments: data_label: D968
type_infered: f64[2]
name: read_data
arguments: data_label: D799
type_infered: f64[2]
name: read_data
arguments: data_label: D366
type_infered: f64
name: read_data
arguments: data_label: D676
type_infered: f64[2]
name: read_data
arguments: data_label: D1035
type_infered: f64[2]
name: read_data
arguments: data_label: D269
type_infered: f64[2]
name: read_data
arguments: data_label: D1047
type_infered: f64[2]
name: read_data
arguments: data_label: D417
type_infered: f64[2]
name: read_data
arguments: data_label: D580
type_infered: f64[2]
name: read_data
arguments: data_label: D95
type_infered: f64[2]
name: read_data
arguments: data_label: D809
type_infered: f64[2]
name: read_data
arguments: data_label: D859
type_infered: f64[2]
name: read_data
arguments: data_label: D118
type_infered: f64[2]
name: read_data
arguments: data_label: D618
type_infered: f64[2]
name: read_data
arguments: data_label: D973
type_infered: f64[2]
name: read_data
arguments: data_label: D504
type_infered: f64[2]
name: read_data
arguments: data_label: D582
type_infered: f64[2]
name: read_data
arguments: data_label: D36
type_infered: f64
name: read_data
arguments: data_label: D895
type_infered: f64[2]
name: read_data
arguments: data_label: D856
type_infered: f64[2]
name: read_data
arguments: data_label: D315
type_infered: f64
name: read_data
arguments: data_label: D1049
type_infered: f64[2]
name: read_data
arguments: data_label: D354
type_infered: f64[2]
name: read_data
arguments: data_label: D84
type_infered: f64[2]
name: read_data
arguments: data_label: D699
type_infered: f64[2]
name: read_data
arguments: data_label: D892
type_infered: f64[2]
name: read_data
arguments: data_label: D910
type_infered: f64[2]
name: read_data
arguments: data_label: D1017
type_infered: f64[2]
name: read_data
arguments: data_label: D334
type_infered: f64[2]
name: read_data
arguments: data_label: D964
type_infered: f64[2]
name: read_data
arguments: data_label: D128
type_infered: f64[2]
name: read_data
arguments: data_label: D829
type_infered: f64[2]
name: read_data
arguments: data_label: D808
type_infered: f64[2]
name: read_data
arguments: data_label: D41
type_infered: f64[2]
name: read_data
arguments: data_label: D14
type_infered: f64[2]
name: read_data
arguments: data_label: D509
type_infered: f64[2]
name: read_data
arguments: data_label: D312
type_infered: f64[2]
name: read_data
arguments: data_label: D530
type_infered: f64[2]
name: read_data
arguments: data_label: D862
type_infered: f64[2]
name: read_data
arguments: data_label: D949
type_infered: f64[2]
name: read_data
arguments: data_label: D891
type_infered: f64[2]
name: read_data
arguments: data_label: D841
type_infered: f64[2]
name: read_data
arguments: data_label: D383
type_infered: f64[2]
name: read_data
arguments: data_label: D115
type_infered: f64[2]
name: read_data
arguments: data_label: D929
type_infered: f64[2]
name: read_data
arguments: data_label: D830
type_infered: f64[2]
name: read_data
arguments: data_label: D181
type_infered: f64[2]
name: read_data
arguments: data_label: D251
type_infered: f64[2]
name: read_data
arguments: data_label: D355
type_infered: f64[2]
name: read_data
arguments: data_label: D437
type_infered: f64
name: read_data
arguments: data_label: D391
type_infered: f64[2]
name: read_data
arguments: data_label: D199
type_infered: f64[2]
name: read_data
arguments: data_label: D998
type_infered: f64[2]
name: read_data
arguments: data_label: D74
type_infered: f64[2]
name: read_data
arguments: data_label: D272
type_infered: f64[2]
name: read_data
arguments: data_label: D771
type_infered: f64[2]
name: read_data
arguments: data_label: D363
type_infered: f64[2]
name: read_data
arguments: data_label: D754
type_infered: f64[2]
name: read_data
arguments: data_label: D53
type_infered: f64[2]
name: read_data
arguments: data_label: D293
type_infered: f64[2]
name: read_data
arguments: data_label: D171
type_infered: f64[2]
name: read_data
arguments: data_label: D529
type_infered: f64[2]
name: read_data
arguments: data_label: D478
type_infered: f64[2]
name: read_data
arguments: data_label: D68
type_infered: f64[2]
name: read_data
arguments: data_label: D955
type_infered: f64[2]
name: read_data
arguments: data_label: D789
type_infered: f64[2]
name: read_data
arguments: data_label: D2
type_infered: f64
name: read_data
arguments: data_label: D375
type_infered: f64[2]
name: read_data
arguments: data_label: D757
type_infered: f64[2]
name: read_data
arguments: data_label: D842
type_infered: f64[2]
name: read_data
arguments: data_label: D227
type_infered: f64[2]
name: read_data
arguments: data_label: D287
type_infered: f64[2]
name: read_data
arguments: data_label: D629
type_infered: f64[2]
name: read_data
arguments: data_label: D195
type_infered: f64[2]
name: read_data
arguments: data_label: D505
type_infered: f64[2]
name: read_data
arguments: data_label: D965
type_infered: f64[2]
name: read_data
arguments: data_label: D289
type_infered: f64[2]
name: read_data
arguments: data_label: D589
type_infered: f64[2]
name: read_data
arguments: data_label: D3
type_infered: f64[2]

tool: {"D693": [50.9848,63.5346]}

tool: D100: 55.1571

tool: {"D863": [73.8998,14.7304]}

tool: D621: 542.942

tool: {"D558": [102.704,52.8191]}

tool: D954: 9.88179

tool: {"D905": [26.5256,3.02949]}

tool: D1071: 56.7717

tool: {"D706": [90.8598,303.55]}

tool: {"D246": [154.364,61.6783]}

tool: {"D794": [158.285,43.9157]}

tool: {"D557": [94.7198,617.416]}

tool: {"D190": [2.1352,546.947]}

tool: {"D797": [141.672,58.4665]}

tool: {"D480": [112.556,86.5584]}

tool: {"D4": [30.6674,33.0052]}

tool: {"D637": [31.6954,213.715]}

tool: D164: 13.819

tool: {"D695": [395.629,93.0676]}

tool: {"D38": [24.8972,429.411]}

tool: {"D945": [57.0025,24.941]}

tool: {"D729": [19.3335,18.6979]}

tool: {"D737": [55.4236,3.9437]}

tool: {"D177": [17.0135,11.3917]}

tool: {"D800": [12.7395,3.04501]}

tool: {"D518": [39.7915,49.7671]}

tool: {"D24": [13.1528,8.48047]}

tool: {"D421": [38.4287,86.9489]}

tool: {"D782": [15.4263,6.13037]}

tool: {"D790": [49.9462,317.097]}

tool: {"D172": [90.6265,535.142]}

tool: {"D968": [301.201,47.9855]}

tool: {"D799": [7.40339,389.49]}

tool: D366: 25.7555

tool: {"D676": [258.443,37.3558]}

tool: {"D1035": [89.0432,32.3832]}

tool: {"D269": [89.3496,68.3162]}

tool: {"D1047": [19.8221,10.3754]}

tool: {"D417": [39.0482,154.526]}

tool: {"D580": [31.5124,17.41]}

tool: {"D95": [40.4141,425.746]}

tool: {"D809": [485.402,93.4276]}

tool: {"D859": [21.0494,10.2711]}

tool: {"D118": [9.61621,4.67304]}

tool: {"D618": [39.802,9.99208]}

tool: {"D973": [8.11339,4.9931]}

tool: {"D504": [67.9371,105.947]}

tool: {"D582": [211.517,8.8836]}

tool: D36: 119.063

tool: {"D895": [143.855,67.9489]}

tool: {"D856": [19.5081,42.3984]}

tool: D315: 161.328

tool: {"D1049": [10.9404,93.3372]}

tool: {"D354": [68.2355,24.357]}

tool: {"D84": [69.6198,9.99288]}

tool: {"D699": [41.5844,55.6896]}

tool: {"D892": [17.1329,16.1444]}

tool: {"D910": [337.079,122.419]}

tool: {"D1017": [397.423,107.176]}

tool: {"D334": [140.766,102.887]}

tool: {"D964": [70.4939,134.182]}

tool: {"D128": [74.0169,31.944]}

tool: {"D829": [34.4739,96.713]}

tool: {"D808": [49.3905,54.5931]}

tool: {"D41": [212.111,55.4413]}

tool: {"D14": [203.006,37.4866]}

tool: {"D509": [51.2984,85.7928]}

tool: {"D312": [82.2765,22.5833]}

tool: {"D530": [47.3307,145.816]}

tool: {"D862": [32.6622,138.698]}

tool: {"D949": [151.612,17.1903]}

tool: {"D891": [79.5384,437.802]}

tool: {"D841": [11.3133,116.392]}

tool: {"D383": [7.19167,12.9767]}

tool: {"D115": [71.707,21.5737]}

tool: {"D929": [15.6557,33.6983]}

tool: {"D830": [26.2569,21.2483]}

tool: {"D181": [66.5309,13.7505]}

tool: {"D251": [51.5706,152.964]}

tool: {"D355": [13.9405,52.3742]}

tool: D437: 11.0864

tool: {"D391": [24.0735,38.1533]}

tool: {"D199": [108.6,32.5711]}

tool: {"D998": [54.3833,22.771]}

tool: {"D74": [19.6303,10.3604]}

tool: {"D272": [387.344,42.017]}

tool: {"D771": [39.8618,13.3264]}

tool: {"D363": [42.8437,645.428]}

tool: {"D754": [21.6046,78.6993]}

tool: {"D53": [29.1148,248.539]}

tool: {"D293": [200.601,25.1849]}

tool: {"D171": [19.6246,15.2335]}

tool: {"D529": [100.214,30.1281]}

tool: {"D478": [58.8731,18.7247]}

tool: {"D68": [11.8358,782.13]}

tool: {"D955": [180.034,30.5434]}

tool: {"D789": [205.845,64.9892]}

tool: D2: 31.473

tool: {"D375": [31.4887,23.2739]}

tool: {"D757": [138.823,36.5369]}

tool: {"D842": [0.0866279,0.043972]}

tool: {"D227": [0.0390894,0.0570451]}

tool: {"D287": [0.0193078,0.0367281]}

tool: {"D629": [0.0580589,0.0832518]}

tool: {"D195": [0.0244313,0.048466]}

tool: {"D505": [0.086209,0.0620286]}

tool: {"D965": [0.0195027,0.0384319]}

tool: {"D289": [0.0457631,0.0695179]}

tool: {"D589": [0.0610127,0.0143859]}

tool: {"D3": [0.0352742,0.0708956]}

assistant: ```c
int Init_Qmat_WAG(double *daa, double *pi) {
  int i, j, naa;
  naa = 20;
  daa[1 * 20 + 0] = 55.15710;
  daa[2 * 20 + 0] = 50.98480;
  daa[2 * 20 + 1] = 63.53460;
  daa[3 * 20 + 0] = 73.89980;
  daa[3 * 20 + 1] = 14.73040;
  daa[3 * 20 + 2] = 542.94200;
  daa[4 * 20 + 0] = 102.70400;
  daa[4 * 20 + 1] = 52.81910;
  daa[4 * 20 + 2] = 26.52560;
  daa[4 * 20 + 3] = 3.02949;
  daa[5 * 20 + 0] = 90.85980;
  daa[5 * 20 + 1] = 303.55000;
  daa[5 * 20 + 2] = 154.36400;
  daa[5 * 20 + 3] = 61.67830;
  daa[5 * 20 + 4] = 9.88179;
  daa[6 * 20 + 0] = 158.28500;
  daa[6 * 20 + 1] = 43.91570;
  daa[6 * 20 + 2] = 94.71980;
  daa[6 * 20 + 3] = 617.41600;
  daa[6 * 20 + 4] = 2.13520;
  daa[6 * 20 + 5] = 546.94700;
  daa[7 * 20 + 0] = 141.67200;
  daa[7 * 20 + 1] = 58.46650;
  daa[7 * 20 + 2] = 112.55600;
  daa[7 * 20 + 3] = 86.55840;
  daa[7 * 20 + 4] = 30.66740;
  daa[7 * 20 + 5] = 33.00520;
  daa[7 * 20 + 6] = 56.77170;
  daa[8 * 20 + 0] = 31.69540;
  daa[8 * 20 + 1] = 213.71500;
  daa[8 * 20 + 2] = 395.62900;
  daa[8 * 20 + 3] = 93.06760;
  daa[8 * 20 + 4] = 24.89720;
  daa[8 * 20 + 5] = 429.41100;
  daa[8 * 20 + 6] = 57.00250;
  daa[8 * 20 + 7] = 24.94100;
  daa[9 * 20 + 0] = 19.33350;
  daa[9 * 20 + 1] = 18.69790;
  daa[9 * 20 + 2] = 55.42360;
  daa[9 * 20 + 3] = 3.94370;
  daa[9 * 20 + 4] = 17.01350;
  daa[9 * 20 + 5] = 11.39170;
  daa[9 * 20 + 6] = 12.73950;
  daa[9 * 20 + 7] = 3.04501;
  daa[9 * 20 + 8] = 13.81900;
  daa[10 * 20 + 0] = 39.79150;
  daa[10 * 20 + 1] = 49.76710;
  daa[10 * 20 + 2] = 13.15280;
  daa[10 * 20 + 3] = 8.48047;
  daa[10 * 20 + 4] = 38.42870;
  daa[10 * 20 + 5] = 86.94890;
  daa[10 * 20 + 6] = 15.42630;
  daa[10 * 20 + 7] = 6.13037;
  daa[10 * 20 + 8] = 49.94620;
  daa[10 * 20 + 9] = 317.09700;
  daa[11 * 20 + 0] = 90.62650;
  daa[11 * 20 + 1] = 535.14200;
  daa[11 * 20 + 2] = 301.20100;
  daa[11 * 20 + 3] = 47.98550;
  daa[11 * 20 + 4] = 7.40339;
  daa[11 * 20 + 5] = 389.49000;
  daa[11 * 20 + 6] = 258.44300;
  daa[11 * 20 + 7] = 37.35580;
  daa[11 * 20 + 8] = 89.04320;
  daa[11 * 20 + 9] = 32.38320;
  daa[11 * 20 + 10] = 25.75550;
  daa[12 * 20 + 0] = 89.34960;
  daa[12 * 20 + 1] = 68.31620;
  daa[12 * 20 + 2] = 19.82210;
  daa[12 * 20 + 3] = 10.37540;
  daa[12 * 20 + 4] = 39.04820;
  daa[12 * 20 + 5] = 154.52600;
  daa[12 * 20 + 6] = 31.51240;
  daa[12 * 20 + 7] = 17.41000;
  daa[12 * 20 + 8] = 40.41410;
  daa[12 * 20 + 9] = 425.74600;
  daa[12 * 20 + 10] = 485.40200;
  daa[12 * 20 + 11] = 93.42760;
  daa[13 * 20 + 0] = 21.04940;
  daa[13 * 20 + 1] = 10.27110;
  daa[13 * 20 + 2] = 9.61621;
  daa[13 * 20 + 3] = 4.67304;
  daa[13 * 20 + 4] = 39.80200;
  daa[13 * 20 + 5] = 9.99208;
  daa[13 * 20 + 6] = 8.11339;
  daa[13 * 20 + 7] = 4.99310;
  daa[13 * 20 + 8] = 67.93710;
  daa[13 * 20 + 9] = 105.94700;
  daa[13 * 20 + 10] = 211.51700;
  daa[13 * 20 + 11] = 8.88360;
  daa[13 * 20 + 12] = 119.06300;
  daa[14 * 20 + 0] = 143.85500;
  daa[14 * 20 + 1] = 67.94890;
  daa[14 * 20 + 2] = 19.50810;
  daa[14 * 20 + 3] = 42.39840;
  daa[14 * 20 + 4] = 10.94040;
  daa[14 * 20 + 5] = 93.33720;
  daa[14 * 20 + 6] = 68.23550;
  daa[14 * 20 + 7] = 24.35700;
  daa[14 * 20 + 8] = 69.61980;
  daa[14 * 20 + 9] = 9.99288;
  daa[14 * 20 + 10] = 41.58440;
  daa[14 * 20 + 11] = 55.68960;
  daa[14 * 20 + 12] = 17.13290;
  daa[14 * 20 + 13] = 16.14440;
  daa[15 * 20 + 0] = 337.07900;
  daa[15 * 20 + 1] = 122.41900;
  daa[15 * 20 + 2] = 397.42300;
  daa[15 * 20 + 3] = 107.17600;
  daa[15 * 20 + 4] = 140.76600;
  daa[15 * 20 + 5] = 102.88700;
  daa[15 * 20 + 6] = 70.49390;
  daa[15 * 20 + 7] = 134.18200;
  daa[15 * 20 + 8] = 74.01690;
  daa[15 * 20 + 9] = 31.94400;
  daa[15 * 20 + 10] = 34.47390;
  daa[15 * 20 + 11] = 96.71300;
  daa[15 * 20 + 12] = 49.39050;
  daa[15 * 20 + 13] = 54.59310;
  daa[15 * 20 + 14] = 161.32800;
  daa[16 * 20 + 0] = 212.11100;
  daa[16 * 20 + 1] = 55.44130;
  daa[16 * 20 + 2] = 203.00600;
  daa[16 * 20 + 3] = 37.48660;
  daa[16 * 20 + 4] = 51.29840;
  daa[16 * 20 + 5] = 85.79280;
  daa[16 * 20 + 6] = 82.27650;
  daa[16 * 20 + 7] = 22.58330;
  daa[16 * 20 + 8] = 47.33070;
  daa[16 * 20 + 9] = 145.81600;
  daa[16 * 20 + 10] = 32.66220;
  daa[16 * 20 + 11] = 138.69800;
  daa[16 * 20 + 12] = 151.61200;
  daa[16 * 20 + 13] = 17.19030;
  daa[16 * 20 + 14] = 79.53840;
  daa[16 * 20 + 15] = 437.80200;
  daa[17 * 20 + 0] = 11.31330;
  daa[17 * 20 + 1] = 116.39200;
  daa[17 * 20 + 2] = 7.19167;
  daa[17 * 20 + 3] = 12.97670;
  daa[17 * 20 + 4] = 71.70700;
  daa[17 * 20 + 5] = 21.57370;
  daa[17 * 20 + 6] = 15.65570;
  daa[17 * 20 + 7] = 33.69830;
  daa[17 * 20 + 8] = 26.25690;
  daa[17 * 20 + 9] = 21.24830;
  daa[17 * 20 + 10] = 66.53090;
  daa[17 * 20 + 11] = 13.75050;
  daa[17 * 20 + 12] = 51.57060;
  daa[17 * 20 + 13] = 152.96400;
  daa[17 * 20 + 14] = 13.94050;
  daa[17 * 20 + 15] = 52.37420;
  daa[17 * 20 + 16] = 11.08640;
  daa[18 * 20 + 0] = 24.07350;
  daa[18 * 20 + 1] = 38.15330;
  daa[18 * 20 + 2] = 108.60000;
  daa[18 * 20 + 3] = 32.57110;
  daa[18 * 20 + 4] = 54.38330;
  daa[18 * 20 + 5] = 22.77100;
  daa[18 * 20 + 6] = 19.63030;
  daa[18 * 20 + 7] = 10.36040;
  daa[18 * 20 + 8] = 387.34400;
  daa[18 * 20 + 9] = 42.01700;
  daa[18 * 20 + 10] = 39.86180;
  daa[18 * 20 + 11] = 13.32640;
  daa[18 * 20 + 12] = 42.84370;
  daa[18 * 20 + 13] = 645.42800;
  daa[18 * 20 + 14] = 21.60460;
  daa[18 * 20 + 15] = 78.69930;
  daa[18 * 20 + 16] = 29.11480;
  daa[18 * 20 + 17] = 248.53900;
  daa[19 * 20 + 0] = 200.60100;
  daa[19 * 20 + 1] = 25.18490;
  daa[19 * 20 + 2] = 19.62460;
  daa[19 * 20 + 3] = 15.23350;
  daa[19 * 20 + 4] = 100.21400;
  daa[19 * 20 + 5] = 30.12810;
  daa[19 * 20 + 6] = 58.87310;
  daa[19 * 20 + 7] = 18.72470;
  daa[19 * 20 + 8] = 11.83580;
  daa[19 * 20 + 9] = 782.13000;
  daa[19 * 20 + 10] = 180.03400;
  daa[19 * 20 + 11] = 30.54340;
  daa[19 * 20 + 12] = 205.84500;
  daa[19 * 20 + 13] = 64.98920;
  daa[19 * 20 + 14] = 31.48870;
  daa[19 * 20 + 15] = 23.27390;
  daa[19 * 20 + 16] = 138.82300;
  daa[19 * 20 + 17] = 36.53690;
  daa[19 * 20 + 18] = 31.47300;
  for (i = 0; i < naa; i++)
    for (j = 0; j < i; j++)
      daa[j * naa + i] = daa[i * naa + j];
  pi[0] = 0.0866279;
  pi[1] = 0.043972;
  pi[2] = 0.0390894;
  pi[3] = 0.0570451;
  pi[4] = 0.0193078;
  pi[5] = 0.0367281;
  pi[6] = 0.0580589;
  pi[7] = 0.0832518;
  pi[8] = 0.0244313;
  pi[9] = 0.048466;
  pi[10] = 0.086209;
  pi[11] = 0.0620286;
  pi[12] = 0.0195027;
  pi[13] = 0.0384319;
  pi[14] = 0.0457631;
  pi[15] = 0.0695179;
  pi[16] = 0.0610127;
  pi[17] = 0.0143859;
  pi[18] = 0.0352742;
  pi[19] = 0.0708956;
  return 1;
}
```