user: Please show me the C implementation that matches the assembly code listed below:

```asm
<Init_Qmat_MtArt>:
  endbr64
  movapd D22(%rip),%xmm2
  mov    %rdi,%rcx
  lea    0x8(%rdi),%r8
  mov    D300(%rip),%rax
  movsd  D351(%rip),%xmm0
  movapd D351(%rip),%xmm1
  lea    0xa0(%rdi),%r9
  add    $0xa8,%rcx
  movups %xmm2,0x1e0(%rdi)
  movapd D557(%rip),%xmm2
  mov    %rax,0x1f0(%rdi)
  mov    D578(%rip),%rax
  movups %xmm2,0x280(%rdi)
  movapd D5(%rip),%xmm2
  mov    %rax,0x490(%rdi)
  movups %xmm2,0x290(%rdi)
  movapd D1025(%rip),%xmm2
  movsd  %xmm0,0xa0(%rdi)
  movups %xmm2,0x320(%rdi)
  movapd D134(%rip),%xmm2
  movups %xmm1,0x140(%rdi)
  movups %xmm2,0x330(%rdi)
  movapd D463(%rip),%xmm2
  movsd  %xmm0,0x340(%rdi)
  movups %xmm2,0x3d0(%rdi)
  movapd D690(%rip),%xmm2
  movups %xmm1,0x3c0(%rdi)
  movups %xmm2,0x3e0(%rdi)
  movapd D938(%rip),%xmm2
  movups %xmm2,0x460(%rdi)
  movapd D737(%rip),%xmm2
  movups %xmm2,0x470(%rdi)
  movapd D659(%rip),%xmm2
  movups %xmm2,0x480(%rdi)
  movapd D384(%rip),%xmm2
  movups %xmm2,0x500(%rdi)
  movapd D570(%rip),%xmm2
  movsd  %xmm0,0x5e0(%rdi)
  movups %xmm2,0x510(%rdi)
  movapd D847(%rip),%xmm2
  movups %xmm2,0x520(%rdi)
  movapd D949(%rip),%xmm2
  movups %xmm2,0x530(%rdi)
  movapd D931(%rip),%xmm2
  movups %xmm2,0x5a0(%rdi)
  movapd D406(%rip),%xmm2
  movups %xmm2,0x5b0(%rdi)
  movapd D68(%rip),%xmm2
  movups %xmm2,0x5c0(%rdi)
  movapd D446(%rip),%xmm2
  movups %xmm2,0x5d0(%rdi)
  movapd D807(%rip),%xmm2
  movups %xmm2,0x640(%rdi)
  movapd D460(%rip),%xmm2
  movups %xmm2,0x650(%rdi)
  movapd D950(%rip),%xmm2
  movups %xmm2,0x660(%rdi)
  movapd D1035(%rip),%xmm2
  movups %xmm2,0x670(%rdi)
  movapd D1030(%rip),%xmm2
  movups %xmm2,0x680(%rdi)
  movapd D579(%rip),%xmm2
  movups %xmm2,0x6e0(%rdi)
  movapd D1016(%rip),%xmm2
  movups %xmm2,0x6f0(%rdi)
  movapd D144(%rip),%xmm2
  movups %xmm2,0x700(%rdi)
  movapd D442(%rip),%xmm2
  movups %xmm2,0x710(%rdi)
  movapd D206(%rip),%xmm2
  mov    D404(%rip),%rax
  movups %xmm2,0x720(%rdi)
  movapd D187(%rip),%xmm2
  mov    %rax,0x730(%rdi)
  mov    D787(%rip),%rax
  movups %xmm2,0x780(%rdi)
  movapd D978(%rip),%xmm2
  mov    %rax,0x880(%rdi)
  movups %xmm2,0x790(%rdi)
  movapd D915(%rip),%xmm2
  movups %xmm2,0x7a0(%rdi)
  movapd D980(%rip),%xmm2
  movups %xmm2,0x7b0(%rdi)
  movapd D692(%rip),%xmm2
  movups %xmm2,0x7c0(%rdi)
  movapd D278(%rip),%xmm2
  movups %xmm2,0x7d0(%rdi)
  movapd D393(%rip),%xmm2
  movups %xmm2,0x820(%rdi)
  movapd D972(%rip),%xmm2
  movups %xmm2,0x830(%rdi)
  movapd D281(%rip),%xmm2
  movups %xmm2,0x840(%rdi)
  movapd D423(%rip),%xmm2
  movups %xmm2,0x850(%rdi)
  movapd D70(%rip),%xmm2
  movups %xmm2,0x860(%rdi)
  movapd D824(%rip),%xmm2
  movups %xmm2,0x870(%rdi)
  movapd D124(%rip),%xmm2
  movups %xmm2,0x8c0(%rdi)
  movapd D382(%rip),%xmm2
  mov    D331(%rip),%rax
  movups %xmm2,0x8d0(%rdi)
  movapd D743(%rip),%xmm2
  mov    %rax,0x9d0(%rdi)
  movups %xmm2,0x8e0(%rdi)
  movapd D859(%rip),%xmm2
  movups %xmm2,0x8f0(%rdi)
  movapd D775(%rip),%xmm2
  movups %xmm2,0x900(%rdi)
  movapd D933(%rip),%xmm2
  movups %xmm2,0x910(%rdi)
  movapd D250(%rip),%xmm2
  movups %xmm2,0x920(%rdi)
  movapd D957(%rip),%xmm2
  movups %xmm2,0x960(%rdi)
  movapd D569(%rip),%xmm2
  movups %xmm2,0x970(%rdi)
  movapd D508(%rip),%xmm2
  movups %xmm2,0x980(%rdi)
  movapd D504(%rip),%xmm2
  movups %xmm2,0x990(%rdi)
  movapd D802(%rip),%xmm2
  movups %xmm2,0x9a0(%rdi)
  movapd D383(%rip),%xmm2
  movups %xmm2,0x9b0(%rdi)
  movapd D819(%rip),%xmm2
  movups %xmm2,0x9c0(%rdi)
  movapd D432(%rip),%xmm2
  movups %xmm2,0xa00(%rdi)
  movapd D793(%rip),%xmm2
  movups %xmm2,0xa10(%rdi)
  movapd D55(%rip),%xmm2
  movups %xmm1,0xaa0(%rdi)
  movups %xmm2,0xa20(%rdi)
  movapd D568(%rip),%xmm2
  movups %xmm1,0xae0(%rdi)
  movapd D72(%rip),%xmm1
  movups %xmm2,0xa30(%rdi)
  movapd D777(%rip),%xmm2
  movsd  %xmm0,0xb20(%rdi)
  movapd D369(%rip),%xmm0
  movups %xmm2,0xa40(%rdi)
  movapd D246(%rip),%xmm2
  movups %xmm1,0xaf0(%rdi)
  movapd D567(%rip),%xmm1
  movups %xmm2,0xa50(%rdi)
  movapd D157(%rip),%xmm2
  movups %xmm0,0xb40(%rdi)
  movapd D738(%rip),%xmm0
  movups %xmm2,0xa60(%rdi)
  movapd D1031(%rip),%xmm2
  movups %xmm1,0xb00(%rdi)
  movapd D572(%rip),%xmm1
  movups %xmm2,0xa70(%rdi)
  movapd D122(%rip),%xmm2
  movups %xmm0,0xb50(%rdi)
  movapd D96(%rip),%xmm0
  movups %xmm2,0xab0(%rdi)
  movapd D221(%rip),%xmm2
  movups %xmm1,0xb10(%rdi)
  movups %xmm2,0xac0(%rdi)
  movapd D944(%rip),%xmm2
  movups %xmm2,0xad0(%rdi)
  movups %xmm0,0xb60(%rdi)
  movapd D893(%rip),%xmm0
  mov    D243(%rip),%rax
  movups %xmm0,0xb70(%rdi)
  movapd D706(%rip),%xmm0
  mov    %rax,0xc70(%rdi)
  movups %xmm0,0xb80(%rdi)
  movapd D856(%rip),%xmm0
  movups %xmm0,0xb90(%rdi)
  movapd D975(%rip),%xmm0
  movups %xmm0,0xba0(%rdi)
  movapd D1017(%rip),%xmm0
  movups %xmm0,0xbb0(%rdi)
  movapd D547(%rip),%xmm0
  movups %xmm0,0xbc0(%rdi)
  movapd D749(%rip),%xmm0
  movups %xmm0,0xbe0(%rdi)
  movapd D232(%rip),%xmm0
  movups %xmm0,0xbf0(%rdi)
  movapd D1019(%rip),%xmm0
  movups %xmm0,0xc00(%rdi)
  movapd D335(%rip),%xmm0
  movups %xmm0,0xc10(%rdi)
  movapd D1014(%rip),%xmm0
  movups %xmm0,0xc20(%rdi)
  movapd D77(%rip),%xmm0
  movups %xmm0,0xc30(%rdi)
  movapd D73(%rip),%xmm0
  movups %xmm0,0xc40(%rdi)
  movapd D928(%rip),%xmm0
  movups %xmm0,0xc50(%rdi)
  movapd D502(%rip),%xmm0
  movups %xmm0,0xc60(%rdi)
  mov    $0x1,%edi
  nopl   0x0(%rax,%rax,1)
L2:
  mov    %r8,%rdx
  mov    %r9,%rax
  cs nopw 0x0(%rax,%rax,1)
L18:
  movsd  (%rax),%xmm0
  add    $0x8,%rax
  add    $0xa0,%rdx
  movsd  %xmm0,-0xa0(%rdx)
  cmp    %rax,%rcx
  jne    L18
  add    $0x1,%edi
  add    $0xa0,%r9
  add    $0x8,%r8
  add    $0xa8,%rcx
  cmp    $0x14,%edi
  jne    L2
  movapd D626(%rip),%xmm0
  mov    $0x1,%eax
  movups %xmm0,(%rsi)
  movapd D506(%rip),%xmm0
  movups %xmm0,0x10(%rsi)
  movapd D435(%rip),%xmm0
  movups %xmm0,0x20(%rsi)
  movapd D622(%rip),%xmm0
  movups %xmm0,0x30(%rsi)
  movapd D602(%rip),%xmm0
  movups %xmm0,0x40(%rsi)
  movapd D66(%rip),%xmm0
  movups %xmm0,0x50(%rsi)
  movapd D729(%rip),%xmm0
  movups %xmm0,0x60(%rsi)
  movapd D420(%rip),%xmm0
  movups %xmm0,0x70(%rsi)
  movapd D969(%rip),%xmm0
  movups %xmm0,0x80(%rsi)
  movapd D731(%rip),%xmm0
  movups %xmm0,0x90(%rsi)
  ret
```

tool_calls: name: read_data
arguments: tag: D22
type: f64[2]
name: read_data
arguments: tag: D300
type: f64
name: read_data
arguments: tag: D351
type: f64[2]
name: read_data
arguments: tag: D557
type: f64[2]
name: read_data
arguments: tag: D578
type: f64
name: read_data
arguments: tag: D5
type: f64[2]
name: read_data
arguments: tag: D1025
type: f64[2]
name: read_data
arguments: tag: D134
type: f64[2]
name: read_data
arguments: tag: D463
type: f64[2]
name: read_data
arguments: tag: D690
type: f64[2]
name: read_data
arguments: tag: D938
type: f64[2]
name: read_data
arguments: tag: D737
type: f64[2]
name: read_data
arguments: tag: D659
type: f64[2]
name: read_data
arguments: tag: D384
type: f64[2]
name: read_data
arguments: tag: D570
type: f64[2]
name: read_data
arguments: tag: D847
type: f64[2]
name: read_data
arguments: tag: D949
type: f64[2]
name: read_data
arguments: tag: D931
type: f64[2]
name: read_data
arguments: tag: D406
type: f64[2]
name: read_data
arguments: tag: D68
type: f64[2]
name: read_data
arguments: tag: D446
type: f64[2]
name: read_data
arguments: tag: D807
type: f64[2]
name: read_data
arguments: tag: D460
type: f64[2]
name: read_data
arguments: tag: D950
type: f64[2]
name: read_data
arguments: tag: D1035
type: f64[2]
name: read_data
arguments: tag: D1030
type: f64[2]
name: read_data
arguments: tag: D579
type: f64[2]
name: read_data
arguments: tag: D1016
type: f64[2]
name: read_data
arguments: tag: D144
type: f64[2]
name: read_data
arguments: tag: D442
type: f64[2]
name: read_data
arguments: tag: D206
type: f64[2]
name: read_data
arguments: tag: D404
type: f64
name: read_data
arguments: tag: D187
type: f64[2]
name: read_data
arguments: tag: D787
type: f64
name: read_data
arguments: tag: D978
type: f64[2]
name: read_data
arguments: tag: D915
type: f64[2]
name: read_data
arguments: tag: D980
type: f64[2]
name: read_data
arguments: tag: D692
type: f64[2]
name: read_data
arguments: tag: D278
type: f64[2]
name: read_data
arguments: tag: D393
type: f64[2]
name: read_data
arguments: tag: D972
type: f64[2]
name: read_data
arguments: tag: D281
type: f64[2]
name: read_data
arguments: tag: D423
type: f64[2]
name: read_data
arguments: tag: D70
type: f64[2]
name: read_data
arguments: tag: D824
type: f64[2]
name: read_data
arguments: tag: D124
type: f64[2]
name: read_data
arguments: tag: D382
type: f64[2]
name: read_data
arguments: tag: D331
type: f64
name: read_data
arguments: tag: D743
type: f64[2]
name: read_data
arguments: tag: D859
type: f64[2]
name: read_data
arguments: tag: D775
type: f64[2]
name: read_data
arguments: tag: D933
type: f64[2]
name: read_data
arguments: tag: D250
type: f64[2]
name: read_data
arguments: tag: D957
type: f64[2]
name: read_data
arguments: tag: D569
type: f64[2]
name: read_data
arguments: tag: D508
type: f64[2]
name: read_data
arguments: tag: D504
type: f64[2]
name: read_data
arguments: tag: D802
type: f64[2]
name: read_data
arguments: tag: D383
type: f64[2]
name: read_data
arguments: tag: D819
type: f64[2]
name: read_data
arguments: tag: D432
type: f64[2]
name: read_data
arguments: tag: D793
type: f64[2]
name: read_data
arguments: tag: D55
type: f64[2]
name: read_data
arguments: tag: D568
type: f64[2]
name: read_data
arguments: tag: D72
type: f64[2]
name: read_data
arguments: tag: D777
type: f64[2]
name: read_data
arguments: tag: D369
type: f64[2]
name: read_data
arguments: tag: D246
type: f64[2]
name: read_data
arguments: tag: D567
type: f64[2]
name: read_data
arguments: tag: D157
type: f64[2]
name: read_data
arguments: tag: D738
type: f64[2]
name: read_data
arguments: tag: D1031
type: f64[2]
name: read_data
arguments: tag: D572
type: f64[2]
name: read_data
arguments: tag: D122
type: f64[2]
name: read_data
arguments: tag: D96
type: f64[2]
name: read_data
arguments: tag: D221
type: f64[2]
name: read_data
arguments: tag: D944
type: f64[2]
name: read_data
arguments: tag: D893
type: f64[2]
name: read_data
arguments: tag: D243
type: f64
name: read_data
arguments: tag: D706
type: f64[2]
name: read_data
arguments: tag: D856
type: f64[2]
name: read_data
arguments: tag: D975
type: f64[2]
name: read_data
arguments: tag: D1017
type: f64[2]
name: read_data
arguments: tag: D547
type: f64[2]
name: read_data
arguments: tag: D749
type: f64[2]
name: read_data
arguments: tag: D232
type: f64[2]
name: read_data
arguments: tag: D1019
type: f64[2]
name: read_data
arguments: tag: D335
type: f64[2]
name: read_data
arguments: tag: D1014
type: f64[2]
name: read_data
arguments: tag: D77
type: f64[2]
name: read_data
arguments: tag: D73
type: f64[2]
name: read_data
arguments: tag: D928
type: f64[2]
name: read_data
arguments: tag: D502
type: f64[2]
name: read_data
arguments: tag: D626
type: f64[2]
name: read_data
arguments: tag: D506
type: f64[2]
name: read_data
arguments: tag: D435
type: f64[2]
name: read_data
arguments: tag: D622
type: f64[2]
name: read_data
arguments: tag: D602
type: f64[2]
name: read_data
arguments: tag: D66
type: f64[2]
name: read_data
arguments: tag: D729
type: f64[2]
name: read_data
arguments: tag: D420
type: f64[2]
name: read_data
arguments: tag: D969
type: f64[2]
name: read_data
arguments: tag: D731
type: f64[2]

tool: {"D22": [0.6,4.3]}

tool: D300: 500.2

tool: {"D351": [0.2,0.2]}

tool: {"D557": [253.5,35.5]}

tool: D578: 43.9

tool: {"D5": [98.2,10.6]}

tool: {"D1025": [0.2,154.0]}

tool: {"D134": [261.8,0.2]}

tool: {"D463": [183.0,861.8]}

tool: {"D690": [0.2,261.6]}

tool: {"D938": [199.8,0.2]}

tool: {"D737": [120.5,12.5]}

tool: {"D659": [80.5,2.6]}

tool: {"D384": [0.2,41.3]}

tool: {"D570": [179.5,0.2]}

tool: {"D847": [12.4,313.5]}

tool: {"D949": [15.2,0.2]}

tool: {"D931": [25.7,1.8]}

tool: {"D406": [21.3,6.6]}

tool: {"D68": [63.0,10.5]}

tool: {"D446": [6.8,2.7]}

tool: {"D807": [3.7,1.8]}

tool: {"D460": [12.6,1.2]}

tool: {"D950": [78.7,16.3]}

tool: {"D1035": [1.7,1.4]}

tool: {"D1030": [5.5,514.5]}

tool: {"D579": [0.2,208.6]}

tool: {"D1016": [467.3,1.7]}

tool: {"D144": [0.2,349.3]}

tool: {"D442": [106.3,0.2]}

tool: {"D206": [0.2,3.5]}

tool: D404: 3.8

tool: {"D187": [120.6,5.2]}

tool: D787: 321.6

tool: {"D978": [78.8,0.2]}

tool: {"D915": [312.3,67.3]}

tool: {"D980": [0.2,55.7]}

tool: {"D692": [0.2,514.8]}

tool: {"D278": [885.5,105.6]}

tool: {"D393": [13.1,4.7]}

tool: {"D972": [19.7,0.2]}

tool: {"D281": [184.1,0.2]}

tool: {"D423": [0.2,0.8]}

tool: {"D70": [13.8,117.9]}

tool: {"D824": [262.6,10.7]}

tool: {"D124": [49.3,0.2]}

tool: {"D382": [16.5,0.2]}

tool: D331: 86.5

tool: {"D743": [0.2,39.3]}

tool: {"D859": [7.9,0.2]}

tool: {"D775": [0.8,0.2]}

tool: {"D933": [12.2,16.8]}

tool: {"D250": [5.3,14.6]}

tool: {"D957": [673.0,2.7]}

tool: {"D569": [398.4,44.4]}

tool: {"D508": [664.2,52.4]}

tool: {"D504": [31.5,226.0]}

tool: {"D802": [10.6,7.2]}

tool: {"D383": [8.2,144.2]}

tool: {"D819": [111.7,36.1]}

tool: {"D432": [243.9,0.2]}

tool: {"D793": [165.9,0.2]}

tool: {"D55": [182.8,43.7]}

tool: {"D568": [43.4,0.2]}

tool: {"D72": [21.1,16.0]}

tool: {"D777": [18.6,203.7]}

tool: {"D369": [1.2,3.9]}

tool: {"D246": [47.8,69.5]}

tool: {"D567": [70.7,53.7]}

tool: {"D157": [288.6,13.5]}

tool: {"D738": [251.2,0.2]}

tool: {"D1031": [46.8,660.4]}

tool: {"D572": [0.2,2.4]}

tool: {"D122": [7.7,0.2]}

tool: {"D96": [72.0,86.7]}

tool: {"D221": [21.6,6.7]}

tool: {"D944": [11.0,1.9]}

tool: {"D893": [7.7,8.6]}

tool: D243: 1.6

tool: {"D706": [191.4,12.3]}

tool: {"D856": [19.8,117.1]}

tool: {"D975": [70.9,791.6]}

tool: {"D1017": [18.4,30.5]}

tool: {"D547": [46.0,37.7]}

tool: {"D749": [339.9,0.2]}

tool: {"D232": [22.6,0.2]}

tool: {"D1019": [350.4,0.2]}

tool: {"D335": [13.6,2.6]}

tool: {"D1014": [0.2,1854.5]}

tool: {"D77": [84.7,26.1]}

tool: {"D73": [281.3,51.9]}

tool: {"D928": [31.7,60.6]}

tool: {"D502": [544.1,0.2]}

tool: {"D626": [0.054116,0.018227]}

tool: {"D506": [0.039903,0.02016]}

tool: {"D435": [0.009709,0.018781]}

tool: {"D622": [0.024289,0.068183]}

tool: {"D602": [0.024518,0.092639]}

tool: {"D66": [0.148658,0.021718]}

tool: {"D729": [0.061453,0.088668]}

tool: {"D420": [0.041826,0.09103]}

tool: {"D969": [0.049194,0.029786]}

tool: {"D731": [0.039443,0.057701]}

assistant: ```c
int Init_Qmat_MtArt(double *daa, double *pi) {
  int i, j, naa;
  naa = 20;
  daa[1 * 20 + 0] = 0.2;
  daa[2 * 20 + 0] = 0.2;
  daa[2 * 20 + 1] = 0.2;
  daa[3 * 20 + 0] = 0.6;
  daa[3 * 20 + 1] = 4.3;
  daa[3 * 20 + 2] = 500.2;
  daa[4 * 20 + 0] = 253.5;
  daa[4 * 20 + 1] = 35.5;
  daa[4 * 20 + 2] = 98.2;
  daa[4 * 20 + 3] = 10.6;
  daa[5 * 20 + 0] = 0.2;
  daa[5 * 20 + 1] = 154.0;
  daa[5 * 20 + 2] = 261.8;
  daa[5 * 20 + 3] = 0.2;
  daa[5 * 20 + 4] = 0.2;
  daa[6 * 20 + 0] = 0.2;
  daa[6 * 20 + 1] = 0.2;
  daa[6 * 20 + 2] = 183.0;
  daa[6 * 20 + 3] = 861.8;
  daa[6 * 20 + 4] = 0.2;
  daa[6 * 20 + 5] = 261.6;
  daa[7 * 20 + 0] = 199.8;
  daa[7 * 20 + 1] = 0.2;
  daa[7 * 20 + 2] = 120.5;
  daa[7 * 20 + 3] = 12.5;
  daa[7 * 20 + 4] = 80.5;
  daa[7 * 20 + 5] = 2.6;
  daa[7 * 20 + 6] = 43.9;
  daa[8 * 20 + 0] = 0.2;
  daa[8 * 20 + 1] = 41.3;
  daa[8 * 20 + 2] = 179.5;
  daa[8 * 20 + 3] = 0.2;
  daa[8 * 20 + 4] = 12.4;
  daa[8 * 20 + 5] = 313.5;
  daa[8 * 20 + 6] = 15.2;
  daa[8 * 20 + 7] = 0.2;
  daa[9 * 20 + 0] = 25.7;
  daa[9 * 20 + 1] = 1.8;
  daa[9 * 20 + 2] = 21.3;
  daa[9 * 20 + 3] = 6.6;
  daa[9 * 20 + 4] = 63.0;
  daa[9 * 20 + 5] = 10.5;
  daa[9 * 20 + 6] = 6.8;
  daa[9 * 20 + 7] = 2.7;
  daa[9 * 20 + 8] = 0.2;
  daa[10 * 20 + 0] = 3.7;
  daa[10 * 20 + 1] = 1.8;
  daa[10 * 20 + 2] = 12.6;
  daa[10 * 20 + 3] = 1.2;
  daa[10 * 20 + 4] = 78.7;
  daa[10 * 20 + 5] = 16.3;
  daa[10 * 20 + 6] = 1.7;
  daa[10 * 20 + 7] = 1.4;
  daa[10 * 20 + 8] = 5.5;
  daa[10 * 20 + 9] = 514.5;
  daa[11 * 20 + 0] = 0.2;
  daa[11 * 20 + 1] = 208.6;
  daa[11 * 20 + 2] = 467.3;
  daa[11 * 20 + 3] = 1.7;
  daa[11 * 20 + 4] = 0.2;
  daa[11 * 20 + 5] = 349.3;
  daa[11 * 20 + 6] = 106.3;
  daa[11 * 20 + 7] = 0.2;
  daa[11 * 20 + 8] = 0.2;
  daa[11 * 20 + 9] = 3.5;
  daa[11 * 20 + 10] = 3.8;
  daa[12 * 20 + 0] = 120.6;
  daa[12 * 20 + 1] = 5.2;
  daa[12 * 20 + 2] = 78.8;
  daa[12 * 20 + 3] = 0.2;
  daa[12 * 20 + 4] = 312.3;
  daa[12 * 20 + 5] = 67.3;
  daa[12 * 20 + 6] = 0.2;
  daa[12 * 20 + 7] = 55.7;
  daa[12 * 20 + 8] = 0.2;
  daa[12 * 20 + 9] = 514.8;
  daa[12 * 20 + 10] = 885.5;
  daa[12 * 20 + 11] = 105.6;
  daa[13 * 20 + 0] = 13.1;
  daa[13 * 20 + 1] = 4.7;
  daa[13 * 20 + 2] = 19.7;
  daa[13 * 20 + 3] = 0.2;
  daa[13 * 20 + 4] = 184.1;
  daa[13 * 20 + 5] = 0.2;
  daa[13 * 20 + 6] = 0.2;
  daa[13 * 20 + 7] = 0.8;
  daa[13 * 20 + 8] = 13.8;
  daa[13 * 20 + 9] = 117.9;
  daa[13 * 20 + 10] = 262.6;
  daa[13 * 20 + 11] = 10.7;
  daa[13 * 20 + 12] = 321.6;
  daa[14 * 20 + 0] = 49.3;
  daa[14 * 20 + 1] = 0.2;
  daa[14 * 20 + 2] = 16.5;
  daa[14 * 20 + 3] = 0.2;
  daa[14 * 20 + 4] = 0.2;
  daa[14 * 20 + 5] = 39.3;
  daa[14 * 20 + 6] = 7.9;
  daa[14 * 20 + 7] = 0.2;
  daa[14 * 20 + 8] = 0.8;
  daa[14 * 20 + 9] = 0.2;
  daa[14 * 20 + 10] = 12.2;
  daa[14 * 20 + 11] = 16.8;
  daa[14 * 20 + 12] = 5.3;
  daa[14 * 20 + 13] = 14.6;
  daa[15 * 20 + 0] = 673.0;
  daa[15 * 20 + 1] = 2.7;
  daa[15 * 20 + 2] = 398.4;
  daa[15 * 20 + 3] = 44.4;
  daa[15 * 20 + 4] = 664.2;
  daa[15 * 20 + 5] = 52.4;
  daa[15 * 20 + 6] = 31.5;
  daa[15 * 20 + 7] = 226.0;
  daa[15 * 20 + 8] = 10.6;
  daa[15 * 20 + 9] = 7.2;
  daa[15 * 20 + 10] = 8.2;
  daa[15 * 20 + 11] = 144.2;
  daa[15 * 20 + 12] = 111.7;
  daa[15 * 20 + 13] = 36.1;
  daa[15 * 20 + 14] = 86.5;
  daa[16 * 20 + 0] = 243.9;
  daa[16 * 20 + 1] = 0.2;
  daa[16 * 20 + 2] = 165.9;
  daa[16 * 20 + 3] = 0.2;
  daa[16 * 20 + 4] = 182.8;
  daa[16 * 20 + 5] = 43.7;
  daa[16 * 20 + 6] = 43.4;
  daa[16 * 20 + 7] = 0.2;
  daa[16 * 20 + 8] = 18.6;
  daa[16 * 20 + 9] = 203.7;
  daa[16 * 20 + 10] = 47.8;
  daa[16 * 20 + 11] = 69.5;
  daa[16 * 20 + 12] = 288.6;
  daa[16 * 20 + 13] = 13.5;
  daa[16 * 20 + 14] = 46.8;
  daa[16 * 20 + 15] = 660.4;
  daa[17 * 20 + 0] = 0.2;
  daa[17 * 20 + 1] = 0.2;
  daa[17 * 20 + 2] = 7.7;
  daa[17 * 20 + 3] = 0.2;
  daa[17 * 20 + 4] = 21.6;
  daa[17 * 20 + 5] = 6.7;
  daa[17 * 20 + 6] = 11.0;
  daa[17 * 20 + 7] = 1.9;
  daa[17 * 20 + 8] = 0.2;
  daa[17 * 20 + 9] = 0.2;
  daa[17 * 20 + 10] = 21.1;
  daa[17 * 20 + 11] = 16.0;
  daa[17 * 20 + 12] = 70.7;
  daa[17 * 20 + 13] = 53.7;
  daa[17 * 20 + 14] = 0.2;
  daa[17 * 20 + 15] = 2.4;
  daa[17 * 20 + 16] = 0.2;
  daa[18 * 20 + 0] = 1.2;
  daa[18 * 20 + 1] = 3.9;
  daa[18 * 20 + 2] = 251.2;
  daa[18 * 20 + 3] = 0.2;
  daa[18 * 20 + 4] = 72.0;
  daa[18 * 20 + 5] = 86.7;
  daa[18 * 20 + 6] = 7.7;
  daa[18 * 20 + 7] = 8.6;
  daa[18 * 20 + 8] = 191.4;
  daa[18 * 20 + 9] = 12.3;
  daa[18 * 20 + 10] = 19.8;
  daa[18 * 20 + 11] = 117.1;
  daa[18 * 20 + 12] = 70.9;
  daa[18 * 20 + 13] = 791.6;
  daa[18 * 20 + 14] = 18.4;
  daa[18 * 20 + 15] = 30.5;
  daa[18 * 20 + 16] = 46.0;
  daa[18 * 20 + 17] = 37.7;
  daa[19 * 20 + 0] = 339.9;
  daa[19 * 20 + 1] = 0.2;
  daa[19 * 20 + 2] = 22.6;
  daa[19 * 20 + 3] = 0.2;
  daa[19 * 20 + 4] = 350.4;
  daa[19 * 20 + 5] = 0.2;
  daa[19 * 20 + 6] = 13.6;
  daa[19 * 20 + 7] = 2.6;
  daa[19 * 20 + 8] = 0.2;
  daa[19 * 20 + 9] = 1854.5;
  daa[19 * 20 + 10] = 84.7;
  daa[19 * 20 + 11] = 26.1;
  daa[19 * 20 + 12] = 281.3;
  daa[19 * 20 + 13] = 51.9;
  daa[19 * 20 + 14] = 31.7;
  daa[19 * 20 + 15] = 60.6;
  daa[19 * 20 + 16] = 544.1;
  daa[19 * 20 + 17] = 0.2;
  daa[19 * 20 + 18] = 1.6;
  for (i = 0; i < naa; i++)
    for (j = 0; j < i; j++)
      daa[j * naa + i] = daa[i * naa + j];
  pi[0] = 0.054116;
  pi[1] = 0.018227;
  pi[2] = 0.039903;
  pi[3] = 0.020160;
  pi[4] = 0.009709;
  pi[5] = 0.018781;
  pi[6] = 0.024289;
  pi[7] = 0.068183;
  pi[8] = 0.024518;
  pi[9] = 0.092639;
  pi[10] = 0.148658;
  pi[11] = 0.021718;
  pi[12] = 0.061453;
  pi[13] = 0.088668;
  pi[14] = 0.041826;
  pi[15] = 0.091030;
  pi[16] = 0.049194;
  pi[17] = 0.029786;
  pi[18] = 0.039443;
  pi[19] = 0.057701;
  return 1;
}
```