Skip to content

Instantly share code, notes, and snippets.

@tailriver
Last active April 16, 2019 10:14
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save tailriver/b262665516558bdb99b8d189116c5c3d to your computer and use it in GitHub Desktop.
Save tailriver/b262665516558bdb99b8d189116c5c3d to your computer and use it in GitHub Desktop.
MUMPS 5.1.2 factorization issue
Entering DMUMPS 5.1.2 from C interface with JOB, N, NNZ = 1 48 1440
executing #MPI = 1 and #OMP = 1
***********CONTROL PARAMETERS (ICNTL)**************
ICNTL(1) Output stream for error messages = 6
ICNTL(2) Output stream for diagnostic messages = 6
ICNTL(3) Output stream for global information = 6
ICNTL(4) Level of printing = 3
ICNTL(5) Matrix format ( keep(55) ) = 0
ICNTL(6) Maximum transversal ( keep(23) ) = 7
ICNTL(7) Ordering = 3
ICNTL(12) LDLT ordering strat ( keep(95) ) = 0
ICNTL(13) Parallel root (0=on, 1=off) = 0
ICNTL(18) Distributed matrix ( keep(54) ) = 0
ICNTL(19) Schur option ( keep(60) 0=off,else=on ) = 0
ICNTL(22) Out-off-core option (0=Off, >0=ON) = 0
ICNTL(8) Scaling strategy = 7
ICNTL(14) Percent of memory increase = 20
=================================================
MUMPS compiled with option -Dmetis
MUMPS compiled with option -Dptscotch
MUMPS compiled with option -DBLR_MT
=================================================
L U Solver for unsymmetric matrices
Type of parallelism: Working host
****** ANALYSIS STEP ********
Resetting candidate strategy to 0 because NSLAVES=1
Entering analysis phase with ...
N NNZ LIW INFO(1)
48 1440 2929 0
Matrix entries: IRN() ICN()
1 1 1 2 1 3
1 4 1 5 1 6
1 7 1 8 1 9
1 10
... Structural symmetry (in percent)= 100
Average density of rows/columns = 29
... No column permutation
Ordering based on SCOTCH
ELAPSED TIME SPENT IN SCOTCH reordering = 0.0006
NFSIZ(.) = 0 0 0 0 0 0 0 0 0 0
FILS (.) = -25 1 2 3 4 5 6 7 8 9
FRERE(.) = 49 49 49 49 49 49 49 49 49 49
Estimated INTEGER space for factors : 195
INFO(3), est. real space to store factors : 1440
Estimated number of entries in factors : 1440
Current value of space relaxation parameter : 20
Estimated size of IS (In Core factorization): 337
Estimated size of S (In Core factorization): 5432
Estimated size of S (OOC factorization) : 8080
Leaving analysis phase with ...
INFOG(1) = 0
INFOG(2) = 0
-- (20) Number of entries in factors (estim.) = 1440
-- (3) Storage of factors (REAL, estimated) = 1440
-- (4) Storage of factors (INT , estimated) = 195
-- (5) Maximum frontal size (estimated) = 24
-- (6) Number of nodes in the tree = 3
-- (32) Type of analysis effectively used = 1
-- (7) Ordering option effectively used = 3
ICNTL(6) Maximum transversal option = 0
ICNTL(7) Pivot order option = 3
Percentage of memory relaxation (effective) = 20
Number of level 2 nodes = 0
Number of split nodes = 0
RINFOG(1) Operations during elimination (estim)= 2.462D+04
Estimated space in MBYTES for IC factorization : 1
** Rank of proc needing largest memory in IC facto : 0
** Estimated corresponding MBYTES for IC facto : 1
** Estimated avg. MBYTES per work. proc at facto (IC) : 1
** TOTAL space in MBYTES for IC factorization : 1
** Rank of proc needing largest memory for OOC facto : 0
** Estimated corresponding MBYTES for OOC facto : 1
** Estimated avg. MBYTES per work. proc at facto (OOC) : 1
** TOTAL space in MBYTES for OOC factorization : 1
ELAPSED TIME IN ANALYSIS DRIVER= 0.0012
***********CONTROL PARAMETERS (ICNTL)**************
ICNTL(1) Output stream for error messages = 6
ICNTL(2) Output stream for diagnostic messages = 6
ICNTL(3) Output stream for global information = 6
ICNTL(4) Level of printing = 3
ICNTL(5) Matrix format ( keep(55) ) = 0
ICNTL(6) Maximum transversal ( keep(23) ) = 7
ICNTL(7) Ordering = 3
ICNTL(12) LDLT ordering strat ( keep(95) ) = 0
ICNTL(13) Parallel root (0=on, 1=off) = 0
ICNTL(18) Distributed matrix ( keep(54) ) = 0
ICNTL(19) Schur option ( keep(60) 0=off,else=on ) = 0
ICNTL(22) Out-off-core option (0=Off, >0=ON) = 0
ICNTL(8) Scaling strategy = 7
ICNTL(14) Percent of memory increase = 20
Entering DMUMPS 5.1.2 from C interface with JOB, N, NNZ = 2 48 1440
executing #MPI = 1 and #OMP = 1
***********CONTROL PARAMETERS (ICNTL)**************
ICNTL(1) Output stream for error messages = 6
ICNTL(2) Output stream for diagnostic messages = 6
ICNTL(3) Output stream for global information = 6
ICNTL(4) Level of printing = 3
ICNTL(8) Scaling strategy = 7
ICNTL(14) Percent of memory increase = 20
****** FACTORIZATION STEP ********
****** FACTORIZATION STEP ********
GLOBAL STATISTICS PRIOR NUMERICAL FACTORIZATION ...
NUMBER OF WORKING PROCESSES = 1
OUT-OF-CORE OPTION (ICNTL(22)) = 0
REAL SPACE FOR FACTORS = 1440
INTEGER SPACE FOR FACTORS = 195
MAXIMUM FRONTAL SIZE (ESTIMATED) = 24
NUMBER OF NODES IN THE TREE = 3
MEMORY ALLOWED (MB -- 0: N/A ) = 0
RELATIVE THRESHOLD FOR PIVOTING, CNTL(1) = 0.1000D-01
Convergence error after scaling for ONE-NORM (option 7/8) = 0.38D-01
Maximum effective relaxed size of S = 5432
Average effective relaxed size of S = 5432
ELAPSED TIME FOR MATRIX DISTRIBUTION = 0.0001
Allocated buffers
------------------
Size of reception buffer in bytes ...... = 140000
Size of async. emission buffer (bytes).. = 140012
Small emission buffer (bytes) .......... = 20
STATISTICS PRIOR NUMERICAL FACTORIZATION ...
Size of internal working array S = 5432
Size of internal working array IS = 337
MINIMUM (ICNTL(14)=0) size of S = 4512
MINIMUM (ICNTL(14)=0) size of IS = 217
REAL SPACE FOR ORIGINAL MATRIX = 1440
INTEGER SPACE FOR ORIGINAL MATRIX = 1536
REAL SPACE FOR FACTORS = 1440
INTEGER SPACE FOR FACTORS = 195
MAXIMUM FRONTAL SIZE (ESTIMATED) = 24
** Space in MBYTES used during factorization : 1
** Memory relaxation parameter ( ICNTL(14) ) : 20
** Rank of processor needing largest memory in facto : 0
** Space in MBYTES used by this processor for facto : 1
** Avg. Space in MBYTES per working proc during facto : 1
LEAVING FACTORIZATION PHASE WITH ...
INFO (1) = 0
--- (2) = 0
NUMBER OF NODES IN THE TREE = 3
INFO (9) REAL SPACE FOR FACTORS = 1440
--- (10) INTEGER SPACE FOR FACTORS = 195
--- (11) MAXIMUM SIZE OF FRONTAL MATRICES = 24
--- (12) NUMBER OF OFF DIAGONAL PIVOTS = 0
--- (13) NUMBER OF DELAYED PIVOTS = 0
--- (14) NUMBER OF MEMORY COMPRESSES = 0
--- (25) NUMBER OF ENTRIES IN FACTORS = 0
RINFO(2) OPERATIONS DURING NODE ASSEMBLY = 2.880D+02
-----(3) OPERATIONS DURING NODE ELIMINATION = 2.462D+04
ELAPSED TIME FOR FACTORIZATION = 0.0064
** Effective minimum Space in MBYTES for facto : 1
Maximum effective space used in S (KEEP8(67)) 1440
Average effective space used in S (KEEP8(67)) 1440
** EFF Min: Rank of processor needing largest memory : 0
** EFF Min: Space in MBYTES used by this processor : 1
** EFF Min: Avg. Space in MBYTES per working proc : 1
GLOBAL STATISTICS
RINFOG(2) OPERATIONS IN NODE ASSEMBLY = 2.880D+02
------(3) OPERATIONS IN NODE ELIMINATION= 2.462D+04
INFOG (9) REAL SPACE FOR FACTORS = 1440
INFOG(10) INTEGER SPACE FOR FACTORS = 195
INFOG(11) MAXIMUM FRONT SIZE = 24
INFOG(29) NUMBER OF ENTRIES IN FACTORS = 1440
INFOG(12) NUMBER OF OFF DIAGONAL PIVOTS = 0
INFOG(13) NUMBER OF DELAYED PIVOTS = 0
INFOG(14) NUMBER OF MEMORY COMPRESS = 0
ELAPSED TIME IN FACTORIZATION DRIVER= 0.0068
***********CONTROL PARAMETERS (ICNTL)**************
ICNTL(1) Output stream for error messages = 6
ICNTL(2) Output stream for diagnostic messages = 6
ICNTL(3) Output stream for global information = 6
ICNTL(4) Level of printing = 3
ICNTL(8) Scaling strategy = 7
ICNTL(14) Percent of memory increase = 20
Entering DMUMPS 5.1.2 from C interface with JOB, N, NNZ = 1 48 1440
executing #MPI = 1 and #OMP = 1
***********CONTROL PARAMETERS (ICNTL)**************
ICNTL(1) Output stream for error messages = 6
ICNTL(2) Output stream for diagnostic messages = 6
ICNTL(3) Output stream for global information = 6
ICNTL(4) Level of printing = 3
ICNTL(5) Matrix format ( keep(55) ) = 0
ICNTL(6) Maximum transversal ( keep(23) ) = 7
ICNTL(7) Ordering = 3
ICNTL(12) LDLT ordering strat ( keep(95) ) = 0
ICNTL(13) Parallel root (0=on, 1=off) = 0
ICNTL(18) Distributed matrix ( keep(54) ) = 0
ICNTL(19) Schur option ( keep(60) 0=off,else=on ) = 0
ICNTL(22) Out-off-core option (0=Off, >0=ON) = 0
ICNTL(8) Scaling strategy = 7
ICNTL(14) Percent of memory increase = 20
=================================================
MUMPS compiled with option -Dmetis
MUMPS compiled with option -Dptscotch
MUMPS compiled with option -DBLR_MT
=================================================
L U Solver for unsymmetric matrices
Type of parallelism: Working host
****** ANALYSIS STEP ********
Resetting candidate strategy to 0 because NSLAVES=1
Entering analysis phase with ...
N NNZ LIW INFO(1)
48 1440 2929 0
Matrix entries: IRN() ICN()
1 1 1 2 1 3
1 4 1 5 1 6
1 7 1 8 1 9
1 10
... Structural symmetry (in percent)= 100
Average density of rows/columns = 29
... No column permutation
Ordering based on SCOTCH
ELAPSED TIME SPENT IN SCOTCH reordering = 0.0006
NFSIZ(.) = 0 0 0 0 0 0 0 0 0 0
FILS (.) = -25 1 2 3 4 5 6 7 8 9
FRERE(.) = 49 49 49 49 49 49 49 49 49 49
Estimated INTEGER space for factors : 195
INFO(3), est. real space to store factors : 1440
Estimated number of entries in factors : 1440
Current value of space relaxation parameter : 20
Estimated size of IS (In Core factorization): 337
Estimated size of S (In Core factorization): 5432
Estimated size of S (OOC factorization) : 8080
Leaving analysis phase with ...
INFOG(1) = 0
INFOG(2) = 0
-- (20) Number of entries in factors (estim.) = 1440
-- (3) Storage of factors (REAL, estimated) = 1440
-- (4) Storage of factors (INT , estimated) = 195
-- (5) Maximum frontal size (estimated) = 24
-- (6) Number of nodes in the tree = 3
-- (32) Type of analysis effectively used = 1
-- (7) Ordering option effectively used = 3
ICNTL(6) Maximum transversal option = 0
ICNTL(7) Pivot order option = 3
Percentage of memory relaxation (effective) = 20
Number of level 2 nodes = 0
Number of split nodes = 0
RINFOG(1) Operations during elimination (estim)= 2.462D+04
Estimated space in MBYTES for IC factorization : 1
** Rank of proc needing largest memory in IC facto : 0
** Estimated corresponding MBYTES for IC facto : 1
** Estimated avg. MBYTES per work. proc at facto (IC) : 1
** TOTAL space in MBYTES for IC factorization : 1
** Rank of proc needing largest memory for OOC facto : 0
** Estimated corresponding MBYTES for OOC facto : 1
** Estimated avg. MBYTES per work. proc at facto (OOC) : 1
** TOTAL space in MBYTES for OOC factorization : 1
ELAPSED TIME IN ANALYSIS DRIVER= 0.0013
***********CONTROL PARAMETERS (ICNTL)**************
ICNTL(1) Output stream for error messages = 6
ICNTL(2) Output stream for diagnostic messages = 6
ICNTL(3) Output stream for global information = 6
ICNTL(4) Level of printing = 3
ICNTL(5) Matrix format ( keep(55) ) = 0
ICNTL(6) Maximum transversal ( keep(23) ) = 7
ICNTL(7) Ordering = 3
ICNTL(12) LDLT ordering strat ( keep(95) ) = 0
ICNTL(13) Parallel root (0=on, 1=off) = 0
ICNTL(18) Distributed matrix ( keep(54) ) = 0
ICNTL(19) Schur option ( keep(60) 0=off,else=on ) = 0
ICNTL(22) Out-off-core option (0=Off, >0=ON) = 0
ICNTL(8) Scaling strategy = 7
ICNTL(14) Percent of memory increase = 20
Entering DMUMPS 5.1.2 from C interface with JOB, N, NNZ = 2 48 1440
executing #MPI = 1 and #OMP = 1
***********CONTROL PARAMETERS (ICNTL)**************
ICNTL(1) Output stream for error messages = 6
ICNTL(2) Output stream for diagnostic messages = 6
ICNTL(3) Output stream for global information = 6
ICNTL(4) Level of printing = 3
ICNTL(8) Scaling strategy = 7
ICNTL(14) Percent of memory increase = 20
****** FACTORIZATION STEP ********
****** FACTORIZATION STEP ********
GLOBAL STATISTICS PRIOR NUMERICAL FACTORIZATION ...
NUMBER OF WORKING PROCESSES = 1
OUT-OF-CORE OPTION (ICNTL(22)) = 0
REAL SPACE FOR FACTORS = 1440
INTEGER SPACE FOR FACTORS = 195
MAXIMUM FRONTAL SIZE (ESTIMATED) = 24
NUMBER OF NODES IN THE TREE = 3
MEMORY ALLOWED (MB -- 0: N/A ) = 0
RELATIVE THRESHOLD FOR PIVOTING, CNTL(1) = 0.1000D-01
Convergence error after scaling for ONE-NORM (option 7/8) = 0.38D-01
Maximum effective relaxed size of S = 5432
Average effective relaxed size of S = 5432
ELAPSED TIME FOR MATRIX DISTRIBUTION = 0.0001
Allocated buffers
------------------
Size of reception buffer in bytes ...... = 140000
Size of async. emission buffer (bytes).. = 140012
Small emission buffer (bytes) .......... = 20
STATISTICS PRIOR NUMERICAL FACTORIZATION ...
Size of internal working array S = 5432
Size of internal working array IS = 337
MINIMUM (ICNTL(14)=0) size of S = 4512
MINIMUM (ICNTL(14)=0) size of IS = 217
REAL SPACE FOR ORIGINAL MATRIX = 1440
INTEGER SPACE FOR ORIGINAL MATRIX = 1536
REAL SPACE FOR FACTORS = 1440
INTEGER SPACE FOR FACTORS = 195
MAXIMUM FRONTAL SIZE (ESTIMATED) = 24
** Space in MBYTES used during factorization : 1
** Memory relaxation parameter ( ICNTL(14) ) : 20
** Rank of processor needing largest memory in facto : 0
** Space in MBYTES used by this processor for facto : 1
** Avg. Space in MBYTES per working proc during facto : 1
LEAVING FACTORIZATION PHASE WITH ...
INFO (1) = 0
--- (2) = 0
NUMBER OF NODES IN THE TREE = 3
INFO (9) REAL SPACE FOR FACTORS = 1464
--- (10) INTEGER SPACE FOR FACTORS = 197
--- (11) MAXIMUM SIZE OF FRONTAL MATRICES = 25
--- (12) NUMBER OF OFF DIAGONAL PIVOTS = 0
--- (13) NUMBER OF DELAYED PIVOTS = 1
--- (14) NUMBER OF MEMORY COMPRESSES = 0
--- (25) NUMBER OF ENTRIES IN FACTORS = 0
RINFO(2) OPERATIONS DURING NODE ASSEMBLY = 3.130D+02
-----(3) OPERATIONS DURING NODE ELIMINATION = 2.549D+04
ELAPSED TIME FOR FACTORIZATION = 0.0051
** Effective minimum Space in MBYTES for facto : 1
Maximum effective space used in S (KEEP8(67)) 1464
Average effective space used in S (KEEP8(67)) 1464
** EFF Min: Rank of processor needing largest memory : 0
** EFF Min: Space in MBYTES used by this processor : 1
** EFF Min: Avg. Space in MBYTES per working proc : 1
GLOBAL STATISTICS
RINFOG(2) OPERATIONS IN NODE ASSEMBLY = 3.130D+02
------(3) OPERATIONS IN NODE ELIMINATION= 2.549D+04
INFOG (9) REAL SPACE FOR FACTORS = 1464
INFOG(10) INTEGER SPACE FOR FACTORS = 197
INFOG(11) MAXIMUM FRONT SIZE = 25
INFOG(29) NUMBER OF ENTRIES IN FACTORS = 1464
INFOG(12) NUMBER OF OFF DIAGONAL PIVOTS = 0
INFOG(13) NUMBER OF DELAYED PIVOTS = 1
INFOG(14) NUMBER OF MEMORY COMPRESS = 0
ELAPSED TIME IN FACTORIZATION DRIVER= 0.0055
***********CONTROL PARAMETERS (ICNTL)**************
ICNTL(1) Output stream for error messages = 6
ICNTL(2) Output stream for diagnostic messages = 6
ICNTL(3) Output stream for global information = 6
ICNTL(4) Level of printing = 3
ICNTL(8) Scaling strategy = 7
ICNTL(14) Percent of memory increase = 20

Correct

OS

Ubuntu 16.04.6 LTS

uname

Linux 9d776e84ab3c 4.4.0-145-generic #171-Ubuntu SMP Tue Mar 26 12:43:40 UTC 2019 x86_64 x86_64 x86_64 GNU/Linux

cpuinfo

processor	: 0
vendor_id	: GenuineIntel
cpu family	: 6
model		: 79
model name	: Intel(R) Xeon(R) CPU E5-1620 v4 @ 3.50GHz
stepping	: 1
microcode	: 0xb00002e
cpu MHz		: 1437.050
cache size	: 10240 KB
physical id	: 0
siblings	: 8
core id		: 0
cpu cores	: 4
apicid		: 0
initial apicid	: 0
fpu		: yes
fpu_exception	: yes
cpuid level	: 20
wp		: yes
flags		: fpu vme de pse tsc msr pae mce cx8 apic sep mtrr pge mca cmov pat pse36 clflush dts acpi mmx fxsr sse sse2 ss ht tm pbe syscall nx pdpe1gb rdtscp lm constant_tsc arch_perfmon pebs bts rep_good nopl xtopology nonstop_tsc aperfmperf pni pclmulqdq dtes64 monitor ds_cpl vmx smx est tm2 ssse3 sdbg fma cx16 xtpr pdcm pcid dca sse4_1 sse4_2 x2apic movbe popcnt tsc_deadline_timer aes xsave avx f16c rdrand lahf_lm abm 3dnowprefetch epb invpcid_single intel_pt ssbd ibrs ibpb stibp kaiser tpr_shadow vnmi flexpriority ept vpid fsgsbase tsc_adjust bmi1 hle avx2 smep bmi2 erms invpcid rtm cqm rdseed adx smap xsaveopt cqm_llc cqm_occup_llc cqm_mbm_total cqm_mbm_local dtherm ida arat pln pts flush_l1d
bugs		: cpu_meltdown spectre_v1 spectre_v2 spec_store_bypass l1tf
bogomips	: 6997.94
clflush size	: 64
cache_alignment	: 64
address sizes	: 46 bits physical, 48 bits virtual
power management:

Incorrect

OS

CentOS Linux release 7.6.1810 (Core)

uname

Linux d504e39aae71 3.10.0-957.5.1.el7.x86_64 #1 SMP Fri Feb 1 14:54:57 UTC 2019 x86_64 x86_64 x86_64 GNU/Linux

cpuinfo

processor	: 0
vendor_id	: GenuineIntel
cpu family	: 6
model		: 58
model name	: Intel(R) Xeon(R) CPU E3-1225 V2 @ 3.20GHz
stepping	: 9
microcode	: 0x20
cpu MHz		: 1674.023
cache size	: 8192 KB
physical id	: 0
siblings	: 4
core id		: 0
cpu cores	: 4
apicid		: 0
initial apicid	: 0
fpu		: yes
fpu_exception	: yes
cpuid level	: 13
wp		: yes
flags		: fpu vme de pse tsc msr pae mce cx8 apic sep mtrr pge mca cmov pat pse36 clflush dts acpi mmx fxsr sse sse2 ss ht tm pbe syscall nx rdtscp lm constant_tsc arch_perfmon pebs bts rep_good nopl xtopology nonstop_tsc aperfmperf eagerfpu pni pclmulqdq dtes64 monitor ds_cpl vmx smx est tm2 ssse3 cx16 xtpr pdcm pcid sse4_1 sse4_2 x2apic popcnt tsc_deadline_timer aes xsave avx f16c rdrand lahf_lm epb ssbd ibrs ibpb stibp tpr_shadow vnmi flexpriority ept vpid fsgsbase smep erms xsaveopt dtherm ida arat pln pts spec_ctrl intel_stibp flush_l1d
bogomips	: 6385.75
clflush size	: 64
cache_alignment	: 64
address sizes	: 36 bits physical, 48 bits virtual
power management:

Docker

Both results get from the same application, libraries (including libdmumps.so) and Docker image (CentOS 6 based).

  • Version: 18.09.5 (both machine)
  • The app even works without container. Virtualization is not root cause of the issue.
1 1 1
1 2 0
1 3 0
1 4 0
1 5 0
1 6 0
1 7 0
1 8 0
1 9 0
1 10 0
1 11 0
1 12 0
1 13 0
1 14 0
1 15 0
1 16 0
1 17 0
1 18 0
1 19 0
1 20 0
1 21 0
1 22 0
1 23 0
1 24 0
2 1 0
2 2 1
2 3 0
2 4 0
2 5 0
2 6 0
2 7 0
2 8 0
2 9 0
2 10 0
2 11 0
2 12 0
2 13 0
2 14 0
2 15 0
2 16 0
2 17 0
2 18 0
2 19 0
2 20 0
2 21 0
2 22 0
2 23 0
2 24 0
3 1 0
3 2 0
3 3 1
3 4 0
3 5 0
3 6 0
3 7 0
3 8 0
3 9 0
3 10 0
3 11 0
3 12 0
3 13 0
3 14 0
3 15 0
3 16 0
3 17 0
3 18 0
3 19 0
3 20 0
3 21 0
3 22 0
3 23 0
3 24 0
4 1 0
4 2 0
4 3 0
4 4 1.23611e+07
4 5 0
4 6 0
4 7 1.80556e+06
4 8 -2.91667e+06
4 9 0
4 10 -4.30556e+06
4 11 2.08333e+06
4 12 0
4 13 -4.30556e+06
4 14 833333
4 15 2.08333e+06
4 16 1.80556e+06
4 17 -1.66667e+06
4 18 -2.91667e+06
4 19 -4.02778e+06
4 20 1.66667e+06
4 21 1.66667e+06
4 22 -972222
4 23 -833333
4 24 -833333
5 1 0
5 2 0
5 3 0
5 4 0
5 5 1
5 6 0
5 7 0
5 8 0
5 9 0
5 10 0
5 11 0
5 12 0
5 13 0
5 14 0
5 15 0
5 16 0
5 17 0
5 18 0
5 19 0
5 20 0
5 21 0
5 22 0
5 23 0
5 24 0
6 1 0
6 2 0
6 3 0
6 4 0
6 5 0
6 6 1
6 7 0
6 8 0
6 9 0
6 10 0
6 11 0
6 12 0
6 13 0
6 14 0
6 15 0
6 16 0
6 17 0
6 18 0
6 19 0
6 20 0
6 21 0
6 22 0
6 23 0
6 24 0
7 1 0
7 2 0
7 3 0
7 4 1.80556e+06
7 5 0
7 6 0
7 7 1.23611e+07
7 8 2.08333e+06
7 9 0
7 10 -2.36111e+06
7 11 -2.91667e+06
7 12 0
7 13 -4.02778e+06
7 14 -1.66667e+06
7 15 1.66667e+06
7 16 -972222
7 17 833333
7 18 -833333
7 19 -4.30556e+06
7 20 -833333
7 21 2.08333e+06
7 22 1.80556e+06
7 23 1.66667e+06
7 24 -2.91667e+06
8 1 0
8 2 0
8 3 0
8 4 -2.91667e+06
8 5 0
8 6 0
8 7 2.08333e+06
8 8 1.23611e+07
8 9 0
8 10 2.91667e+06
8 11 1.80556e+06
8 12 0
8 13 -1.66667e+06
8 14 -4.02778e+06
8 15 1.66667e+06
8 16 -833333
8 17 -4.30556e+06
8 18 2.08333e+06
8 19 833333
8 20 -972222
8 21 -833333
8 22 1.66667e+06
8 23 1.80556e+06
8 24 -2.91667e+06
9 1 0
9 2 0
9 3 0
9 4 0
9 5 0
9 6 0
9 7 0
9 8 0
9 9 1
9 10 0
9 11 0
9 12 0
9 13 0
9 14 0
9 15 0
9 16 0
9 17 0
9 18 0
9 19 0
9 20 0
9 21 0
9 22 0
9 23 0
9 24 0
10 1 0
10 2 0
10 3 0
10 4 -4.30556e+06
10 5 0
10 6 0
10 7 -2.36111e+06
10 8 2.91667e+06
10 9 0
10 10 1.23611e+07
10 11 -2.08333e+06
10 12 0
10 13 -972222
10 14 -833333
10 15 833333
10 16 -4.02778e+06
10 17 1.66667e+06
10 18 -1.66667e+06
10 19 1.80556e+06
10 20 -1.66667e+06
10 21 2.91667e+06
10 22 -4.30556e+06
10 23 833333
10 24 -2.08333e+06
11 1 0
11 2 0
11 3 0
11 4 2.08333e+06
11 5 0
11 6 0
11 7 -2.91667e+06
11 8 1.80556e+06
11 9 0
11 10 -2.08333e+06
11 11 1.23611e+07
11 12 0
11 13 833333
11 14 -4.30556e+06
11 15 2.08333e+06
11 16 1.66667e+06
11 17 -4.02778e+06
11 18 1.66667e+06
11 19 -1.66667e+06
11 20 1.80556e+06
11 21 -2.91667e+06
11 22 -833333
11 23 -972222
11 24 -833333
12 1 0
12 2 0
12 3 0
12 4 0
12 5 0
12 6 0
12 7 0
12 8 0
12 9 0
12 10 0
12 11 0
12 12 1
12 13 0
12 14 0
12 15 0
12 16 0
12 17 0
12 18 0
12 19 0
12 20 0
12 21 0
12 22 0
12 23 0
12 24 0
13 1 0
13 2 0
13 3 0
13 4 -4.30556e+06
13 5 0
13 6 0
13 7 -4.02778e+06
13 8 -1.66667e+06
13 9 0
13 10 -972222
13 11 833333
13 12 0
13 13 1.23611e+07
13 14 2.08333e+06
13 15 -2.08333e+06
13 16 -2.36111e+06
13 17 -2.91667e+06
13 18 2.91667e+06
13 19 1.80556e+06
13 20 2.91667e+06
13 21 -1.66667e+06
13 22 -4.30556e+06
13 23 -2.08333e+06
13 24 833333
13 25 0
13 26 0
13 27 0
13 28 0
13 29 0
13 30 0
13 31 0
13 32 0
13 33 0
13 34 0
13 35 0
13 36 0
14 1 0
14 2 0
14 3 0
14 4 833333
14 5 0
14 6 0
14 7 -1.66667e+06
14 8 -4.02778e+06
14 9 0
14 10 -833333
14 11 -4.30556e+06
14 12 0
14 13 2.08333e+06
14 14 1.23611e+07
14 15 -2.08333e+06
14 16 2.91667e+06
14 17 1.80556e+06
14 18 -1.66667e+06
14 19 -2.91667e+06
14 20 -2.36111e+06
14 21 2.91667e+06
14 22 -2.08333e+06
14 23 -4.30556e+06
14 24 833333
14 25 0
14 26 0
14 27 0
14 28 0
14 29 0
14 30 0
14 31 0
14 32 0
14 33 0
14 34 0
14 35 0
14 36 0
15 1 0
15 2 0
15 3 0
15 4 2.08333e+06
15 5 0
15 6 0
15 7 1.66667e+06
15 8 1.66667e+06
15 9 0
15 10 833333
15 11 2.08333e+06
15 12 0
15 13 -2.08333e+06
15 14 -2.08333e+06
15 15 4.84092e+09
15 16 -2.91667e+06
15 17 -1.66667e+06
15 18 2.89894e+09
15 19 -1.66667e+06
15 20 -2.91667e+06
15 21 2.89894e+09
15 22 -833333
15 23 -833333
15 24 1.73731e+09
15 25 0
15 26 0
15 27 -6.95313e+09
15 28 0
15 29 0
15 30 -2.31771e+09
15 31 0
15 32 0
15 33 -2.31771e+09
15 34 0
15 35 0
15 36 -7.72569e+08
16 1 0
16 2 0
16 3 0
16 4 1.80556e+06
16 5 0
16 6 0
16 7 -972222
16 8 -833333
16 9 0
16 10 -4.02778e+06
16 11 1.66667e+06
16 12 0
16 13 -2.36111e+06
16 14 2.91667e+06
16 15 -2.91667e+06
16 16 1.23611e+07
16 17 -2.08333e+06
16 18 2.08333e+06
16 19 -4.30556e+06
16 20 2.08333e+06
16 21 -833333
16 22 1.80556e+06
16 23 -2.91667e+06
16 24 1.66667e+06
16 25 0
16 26 0
16 27 0
16 28 0
16 29 0
16 30 0
16 31 0
16 32 0
16 33 0
16 34 0
16 35 0
16 36 0
17 1 0
17 2 0
17 3 0
17 4 -1.66667e+06
17 5 0
17 6 0
17 7 833333
17 8 -4.30556e+06
17 9 0
17 10 1.66667e+06
17 11 -4.02778e+06
17 12 0
17 13 -2.91667e+06
17 14 1.80556e+06
17 15 -1.66667e+06
17 16 -2.08333e+06
17 17 1.23611e+07
17 18 -2.08333e+06
17 19 2.08333e+06
17 20 -4.30556e+06
17 21 833333
17 22 2.91667e+06
17 23 -2.36111e+06
17 24 2.91667e+06
17 25 0
17 26 0
17 27 0
17 28 0
17 29 0
17 30 0
17 31 0
17 32 0
17 33 0
17 34 0
17 35 0
17 36 0
18 1 0
18 2 0
18 3 0
18 4 -2.91667e+06
18 5 0
18 6 0
18 7 -833333
18 8 2.08333e+06
18 9 0
18 10 -1.66667e+06
18 11 1.66667e+06
18 12 0
18 13 2.91667e+06
18 14 -1.66667e+06
18 15 2.89894e+09
18 16 2.08333e+06
18 17 -2.08333e+06
18 18 4.84092e+09
18 19 833333
18 20 -833333
18 21 1.73731e+09
18 22 1.66667e+06
18 23 -2.91667e+06
18 24 2.89894e+09
18 25 0
18 26 0
18 27 -2.31771e+09
18 28 0
18 29 0
18 30 -6.95313e+09
18 31 0
18 32 0
18 33 -7.72569e+08
18 34 0
18 35 0
18 36 -2.31771e+09
19 1 0
19 2 0
19 3 0
19 4 -4.02778e+06
19 5 0
19 6 0
19 7 -4.30556e+06
19 8 833333
19 9 0
19 10 1.80556e+06
19 11 -1.66667e+06
19 12 0
19 13 1.80556e+06
19 14 -2.91667e+06
19 15 -1.66667e+06
19 16 -4.30556e+06
19 17 2.08333e+06
19 18 833333
19 19 1.23611e+07
19 20 -2.08333e+06
19 21 -2.08333e+06
19 22 -2.36111e+06
19 23 2.91667e+06
19 24 2.91667e+06
19 25 0
19 26 0
19 27 0
19 28 0
19 29 0
19 30 0
19 31 0
19 32 0
19 33 0
19 34 0
19 35 0
19 36 0
20 1 0
20 2 0
20 3 0
20 4 1.66667e+06
20 5 0
20 6 0
20 7 -833333
20 8 -972222
20 9 0
20 10 -1.66667e+06
20 11 1.80556e+06
20 12 0
20 13 2.91667e+06
20 14 -2.36111e+06
20 15 -2.91667e+06
20 16 2.08333e+06
20 17 -4.30556e+06
20 18 -833333
20 19 -2.08333e+06
20 20 1.23611e+07
20 21 2.08333e+06
20 22 -2.91667e+06
20 23 1.80556e+06
20 24 1.66667e+06
20 25 0
20 26 0
20 27 0
20 28 0
20 29 0
20 30 0
20 31 0
20 32 0
20 33 0
20 34 0
20 35 0
20 36 0
21 1 0
21 2 0
21 3 0
21 4 1.66667e+06
21 5 0
21 6 0
21 7 2.08333e+06
21 8 -833333
21 9 0
21 10 2.91667e+06
21 11 -2.91667e+06
21 12 0
21 13 -1.66667e+06
21 14 2.91667e+06
21 15 2.89894e+09
21 16 -833333
21 17 833333
21 18 1.73731e+09
21 19 -2.08333e+06
21 20 2.08333e+06
21 21 4.84092e+09
21 22 -2.91667e+06
21 23 1.66667e+06
21 24 2.89894e+09
21 25 0
21 26 0
21 27 -2.31771e+09
21 28 0
21 29 0
21 30 -7.72569e+08
21 31 0
21 32 0
21 33 -6.95313e+09
21 34 0
21 35 0
21 36 -2.31771e+09
22 1 0
22 2 0
22 3 0
22 4 -972222
22 5 0
22 6 0
22 7 1.80556e+06
22 8 1.66667e+06
22 9 0
22 10 -4.30556e+06
22 11 -833333
22 12 0
22 13 -4.30556e+06
22 14 -2.08333e+06
22 15 -833333
22 16 1.80556e+06
22 17 2.91667e+06
22 18 1.66667e+06
22 19 -2.36111e+06
22 20 -2.91667e+06
22 21 -2.91667e+06
22 22 1.23611e+07
22 23 2.08333e+06
22 24 2.08333e+06
22 25 0
22 26 0
22 27 0
22 28 0
22 29 0
22 30 0
22 31 0
22 32 0
22 33 0
22 34 0
22 35 0
22 36 0
23 1 0
23 2 0
23 3 0
23 4 -833333
23 5 0
23 6 0
23 7 1.66667e+06
23 8 1.80556e+06
23 9 0
23 10 833333
23 11 -972222
23 12 0
23 13 -2.08333e+06
23 14 -4.30556e+06
23 15 -833333
23 16 -2.91667e+06
23 17 -2.36111e+06
23 18 -2.91667e+06
23 19 2.91667e+06
23 20 1.80556e+06
23 21 1.66667e+06
23 22 2.08333e+06
23 23 1.23611e+07
23 24 2.08333e+06
23 25 0
23 26 0
23 27 0
23 28 0
23 29 0
23 30 0
23 31 0
23 32 0
23 33 0
23 34 0
23 35 0
23 36 0
24 1 0
24 2 0
24 3 0
24 4 -833333
24 5 0
24 6 0
24 7 -2.91667e+06
24 8 -2.91667e+06
24 9 0
24 10 -2.08333e+06
24 11 -833333
24 12 0
24 13 833333
24 14 833333
24 15 1.73731e+09
24 16 1.66667e+06
24 17 2.91667e+06
24 18 2.89894e+09
24 19 2.91667e+06
24 20 1.66667e+06
24 21 2.89894e+09
24 22 2.08333e+06
24 23 2.08333e+06
24 24 4.84092e+09
24 25 0
24 26 0
24 27 -7.72569e+08
24 28 0
24 29 0
24 30 -2.31771e+09
24 31 0
24 32 0
24 33 -2.31771e+09
24 34 0
24 35 0
24 36 -6.95313e+09
25 13 0
25 14 0
25 15 0
25 16 0
25 17 0
25 18 0
25 19 0
25 20 0
25 21 0
25 22 0
25 23 0
25 24 0
25 25 6.18056e+06
25 26 1.04167e+06
25 27 1.04167e+06
25 28 -1.18056e+06
25 29 -1.45833e+06
25 30 -1.45833e+06
25 31 902778
25 32 1.45833e+06
25 33 833333
25 34 -2.15278e+06
25 35 -1.04167e+06
25 36 -416667
25 37 0
25 38 833333
25 39 0
25 40 -2.15278e+06
25 41 -416667
25 42 0
25 43 0
25 44 416667
25 45 0
25 46 -2.01389e+06
25 47 -833333
25 48 0
26 13 0
26 14 0
26 15 0
26 16 0
26 17 0
26 18 0
26 19 0
26 20 0
26 21 0
26 22 0
26 23 0
26 24 0
26 25 1.04167e+06
26 26 6.18056e+06
26 27 1.04167e+06
26 28 1.45833e+06
26 29 902778
26 30 833333
26 31 -1.45833e+06
26 32 -1.18056e+06
26 33 -1.45833e+06
26 34 -1.04167e+06
26 35 -2.15278e+06
26 36 -416667
26 37 0
26 38 902778
26 39 0
26 40 416667
26 41 -486111
26 42 0
26 43 0
26 44 -2.15278e+06
26 45 0
26 46 -833333
26 47 -2.01389e+06
26 48 0
27 13 0
27 14 0
27 15 -6.95313e+09
27 16 0
27 17 0
27 18 -2.31771e+09
27 19 0
27 20 0
27 21 -2.31771e+09
27 22 0
27 23 0
27 24 -7.72569e+08
27 25 1.04167e+06
27 26 1.04167e+06
27 27 1.23673e+10
27 28 1.45833e+06
27 29 833333
27 30 902778
27 31 833333
27 32 1.45833e+06
27 33 902778
27 34 416667
27 35 416667
27 36 -486111
27 37 0
27 38 -1.45833e+06
27 39 0
27 40 -1.04167e+06
27 41 -416667
27 42 0
27 43 0
27 44 -1.04167e+06
27 45 0
27 46 -833333
27 47 -833333
27 48 0
28 13 0
28 14 0
28 15 0
28 16 0
28 17 0
28 18 0
28 19 0
28 20 0
28 21 0
28 22 0
28 23 0
28 24 0
28 25 -1.18056e+06
28 26 1.45833e+06
28 27 1.45833e+06
28 28 6.18056e+06
28 29 -1.04167e+06
28 30 -1.04167e+06
28 31 -2.15278e+06
28 32 1.04167e+06
28 33 416667
28 34 902778
28 35 -1.45833e+06
28 36 -833333
28 37 0
28 38 416667
28 39 0
28 40 902778
28 41 -833333
28 42 0
28 43 0
28 44 833333
28 45 0
28 46 -486111
28 47 -416667
28 48 0
29 13 0
29 14 0
29 15 0
29 16 0
29 17 0
29 18 0
29 19 0
29 20 0
29 21 0
29 22 0
29 23 0
29 24 0
29 25 -1.45833e+06
29 26 902778
29 27 833333
29 28 -1.04167e+06
29 29 6.18056e+06
29 30 1.04167e+06
29 31 1.04167e+06
29 32 -2.15278e+06
29 33 -416667
29 34 1.45833e+06
29 35 -1.18056e+06
29 36 -1.45833e+06
29 37 0
29 38 -486111
29 39 0
29 40 -833333
29 41 902778
29 42 0
29 43 0
29 44 -2.01389e+06
29 45 0
29 46 416667
29 47 -2.15278e+06
29 48 0
30 13 0
30 14 0
30 15 -2.31771e+09
30 16 0
30 17 0
30 18 -6.95313e+09
30 19 0
30 20 0
30 21 -7.72569e+08
30 22 0
30 23 0
30 24 -2.31771e+09
30 25 -1.45833e+06
30 26 833333
30 27 902778
30 28 -1.04167e+06
30 29 1.04167e+06
30 30 1.23673e+10
30 31 -416667
30 32 416667
30 33 -486111
30 34 -833333
30 35 1.45833e+06
30 36 902778
30 37 0
30 38 -416667
30 39 0
30 40 1.45833e+06
30 41 -1.45833e+06
30 42 0
30 43 0
30 44 -833333
30 45 0
30 46 416667
30 47 -1.04167e+06
30 48 0
31 13 0
31 14 0
31 15 0
31 16 0
31 17 0
31 18 0
31 19 0
31 20 0
31 21 0
31 22 0
31 23 0
31 24 0
31 25 902778
31 26 -1.45833e+06
31 27 833333
31 28 -2.15278e+06
31 29 1.04167e+06
31 30 -416667
31 31 6.18056e+06
31 32 -1.04167e+06
31 33 1.04167e+06
31 34 -1.18056e+06
31 35 1.45833e+06
31 36 -1.45833e+06
31 37 0
31 38 -416667
31 39 0
31 40 -2.01389e+06
31 41 833333
31 42 0
31 43 0
31 44 -833333
31 45 0
31 46 -2.15278e+06
31 47 416667
31 48 0
32 13 0
32 14 0
32 15 0
32 16 0
32 17 0
32 18 0
32 19 0
32 20 0
32 21 0
32 22 0
32 23 0
32 24 0
32 25 1.45833e+06
32 26 -1.18056e+06
32 27 1.45833e+06
32 28 1.04167e+06
32 29 -2.15278e+06
32 30 416667
32 31 -1.04167e+06
32 32 6.18056e+06
32 33 -1.04167e+06
32 34 -1.45833e+06
32 35 902778
32 36 -833333
32 37 0
32 38 -2.15278e+06
32 39 0
32 40 833333
32 41 -2.01389e+06
32 42 0
32 43 0
32 44 902778
32 45 0
32 46 -416667
32 47 -486111
32 48 0
33 13 0
33 14 0
33 15 -2.31771e+09
33 16 0
33 17 0
33 18 -7.72569e+08
33 19 0
33 20 0
33 21 -6.95313e+09
33 22 0
33 23 0
33 24 -2.31771e+09
33 25 833333
33 26 -1.45833e+06
33 27 902778
33 28 416667
33 29 -416667
33 30 -486111
33 31 1.04167e+06
33 32 -1.04167e+06
33 33 1.23673e+10
33 34 1.45833e+06
33 35 -833333
33 36 902778
33 37 0
33 38 1.04167e+06
33 39 0
33 40 -833333
33 41 833333
33 42 0
33 43 0
33 44 1.45833e+06
33 45 0
33 46 -1.04167e+06
33 47 416667
33 48 0
34 13 0
34 14 0
34 15 0
34 16 0
34 17 0
34 18 0
34 19 0
34 20 0
34 21 0
34 22 0
34 23 0
34 24 0
34 25 -2.15278e+06
34 26 -1.04167e+06
34 27 416667
34 28 902778
34 29 1.45833e+06
34 30 -833333
34 31 -1.18056e+06
34 32 -1.45833e+06
34 33 1.45833e+06
34 34 6.18056e+06
34 35 1.04167e+06
34 36 -1.04167e+06
34 37 0
34 38 -833333
34 39 0
34 40 -486111
34 41 416667
34 42 0
34 43 0
34 44 -416667
34 45 0
34 46 902778
34 47 833333
34 48 0
35 13 0
35 14 0
35 15 0
35 16 0
35 17 0
35 18 0
35 19 0
35 20 0
35 21 0
35 22 0
35 23 0
35 24 0
35 25 -1.04167e+06
35 26 -2.15278e+06
35 27 416667
35 28 -1.45833e+06
35 29 -1.18056e+06
35 30 1.45833e+06
35 31 1.45833e+06
35 32 902778
35 33 -833333
35 34 1.04167e+06
35 35 6.18056e+06
35 36 -1.04167e+06
35 37 0
35 38 -2.01389e+06
35 39 0
35 40 -416667
35 41 -2.15278e+06
35 42 0
35 43 0
35 44 -486111
35 45 0
35 46 833333
35 47 902778
35 48 0
36 13 0
36 14 0
36 15 -7.72569e+08
36 16 0
36 17 0
36 18 -2.31771e+09
36 19 0
36 20 0
36 21 -2.31771e+09
36 22 0
36 23 0
36 24 -6.95313e+09
36 25 -416667
36 26 -416667
36 27 -486111
36 28 -833333
36 29 -1.45833e+06
36 30 902778
36 31 -1.45833e+06
36 32 -833333
36 33 902778
36 34 -1.04167e+06
36 35 -1.04167e+06
36 36 1.23673e+10
36 37 0
36 38 833333
36 39 0
36 40 416667
36 41 1.04167e+06
36 42 0
36 43 0
36 44 416667
36 45 0
36 46 1.45833e+06
36 47 1.45833e+06
36 48 0
37 25 0
37 26 0
37 27 0
37 28 0
37 29 0
37 30 0
37 31 0
37 32 0
37 33 0
37 34 0
37 35 0
37 36 0
37 37 1
37 38 0
37 39 0
37 40 0
37 41 0
37 42 0
37 43 0
37 44 0
37 45 0
37 46 0
37 47 0
37 48 0
38 25 833333
38 26 902778
38 27 -1.45833e+06
38 28 416667
38 29 -486111
38 30 -416667
38 31 -416667
38 32 -2.15278e+06
38 33 1.04167e+06
38 34 -833333
38 35 -2.01389e+06
38 36 833333
38 37 0
38 38 6.18056e+06
38 39 0
38 40 1.45833e+06
38 41 902778
38 42 0
38 43 0
38 44 -1.18056e+06
38 45 0
38 46 -1.04167e+06
38 47 -2.15278e+06
38 48 0
39 25 0
39 26 0
39 27 0
39 28 0
39 29 0
39 30 0
39 31 0
39 32 0
39 33 0
39 34 0
39 35 0
39 36 0
39 37 0
39 38 0
39 39 1
39 40 0
39 41 0
39 42 0
39 43 0
39 44 0
39 45 0
39 46 0
39 47 0
39 48 0
40 25 -2.15278e+06
40 26 416667
40 27 -1.04167e+06
40 28 902778
40 29 -833333
40 30 1.45833e+06
40 31 -2.01389e+06
40 32 833333
40 33 -833333
40 34 -486111
40 35 -416667
40 36 416667
40 37 0
40 38 1.45833e+06
40 39 0
40 40 6.18056e+06
40 41 -1.04167e+06
40 42 0
40 43 0
40 44 1.04167e+06
40 45 0
40 46 902778
40 47 -1.45833e+06
40 48 0
41 25 -416667
41 26 -486111
41 27 -416667
41 28 -833333
41 29 902778
41 30 -1.45833e+06
41 31 833333
41 32 -2.01389e+06
41 33 833333
41 34 416667
41 35 -2.15278e+06
41 36 1.04167e+06
41 37 0
41 38 902778
41 39 0
41 40 -1.04167e+06
41 41 6.18056e+06
41 42 0
41 43 0
41 44 -2.15278e+06
41 45 0
41 46 1.45833e+06
41 47 -1.18056e+06
41 48 0
42 25 0
42 26 0
42 27 0
42 28 0
42 29 0
42 30 0
42 31 0
42 32 0
42 33 0
42 34 0
42 35 0
42 36 0
42 37 0
42 38 0
42 39 0
42 40 0
42 41 0
42 42 1
42 43 0
42 44 0
42 45 0
42 46 0
42 47 0
42 48 0
43 25 0
43 26 0
43 27 0
43 28 0
43 29 0
43 30 0
43 31 0
43 32 0
43 33 0
43 34 0
43 35 0
43 36 0
43 37 0
43 38 0
43 39 0
43 40 0
43 41 0
43 42 0
43 43 1
43 44 0
43 45 0
43 46 0
43 47 0
43 48 0
44 25 416667
44 26 -2.15278e+06
44 27 -1.04167e+06
44 28 833333
44 29 -2.01389e+06
44 30 -833333
44 31 -833333
44 32 902778
44 33 1.45833e+06
44 34 -416667
44 35 -486111
44 36 416667
44 37 0
44 38 -1.18056e+06
44 39 0
44 40 1.04167e+06
44 41 -2.15278e+06
44 42 0
44 43 0
44 44 6.18056e+06
44 45 0
44 46 -1.45833e+06
44 47 902778
44 48 0
45 25 0
45 26 0
45 27 0
45 28 0
45 29 0
45 30 0
45 31 0
45 32 0
45 33 0
45 34 0
45 35 0
45 36 0
45 37 0
45 38 0
45 39 0
45 40 0
45 41 0
45 42 0
45 43 0
45 44 0
45 45 1
45 46 0
45 47 0
45 48 0
46 25 -2.01389e+06
46 26 -833333
46 27 -833333
46 28 -486111
46 29 416667
46 30 416667
46 31 -2.15278e+06
46 32 -416667
46 33 -1.04167e+06
46 34 902778
46 35 833333
46 36 1.45833e+06
46 37 0
46 38 -1.04167e+06
46 39 0
46 40 902778
46 41 1.45833e+06
46 42 0
46 43 0
46 44 -1.45833e+06
46 45 0
46 46 6.18056e+06
46 47 1.04167e+06
46 48 0
47 25 -833333
47 26 -2.01389e+06
47 27 -833333
47 28 -416667
47 29 -2.15278e+06
47 30 -1.04167e+06
47 31 416667
47 32 -486111
47 33 416667
47 34 833333
47 35 902778
47 36 1.45833e+06
47 37 0
47 38 -2.15278e+06
47 39 0
47 40 -1.45833e+06
47 41 -1.18056e+06
47 42 0
47 43 0
47 44 902778
47 45 0
47 46 1.04167e+06
47 47 6.18056e+06
47 48 0
48 25 0
48 26 0
48 27 0
48 28 0
48 29 0
48 30 0
48 31 0
48 32 0
48 33 0
48 34 0
48 35 0
48 36 0
48 37 0
48 38 0
48 39 0
48 40 0
48 41 0
48 42 0
48 43 0
48 44 0
48 45 0
48 46 0
48 47 0
48 48 1
#
# This file is part of MUMPS 5.1.2, released
# on Mon Oct 2 07:37:01 UTC 2017
#
#Begin orderings
# NOTE that PORD is distributed within MUMPS by default. It is recommended to
# install other orderings. For that, you need to obtain the corresponding package
# and modify the variables below accordingly.
# For example, to have Metis available within MUMPS:
# 1/ download Metis and compile it
# 2/ uncomment (suppress # in first column) lines
# starting with LMETISDIR, LMETIS
# 3/ add -Dmetis in line ORDERINGSF
# ORDERINGSF = -Dpord -Dmetis
# 4/ Compile and install MUMPS
# make clean; make (to clean up previous installation)
#
# Metis/ParMetis and SCOTCH/PT-SCOTCH (ver 6.0 and later) orderings are recommended.
#
SCOTCHDIR = /usr/local
ISCOTCH = -I$(SCOTCHDIR)/include
#
# You have to choose one among the following two lines depending on
# the type of analysis you want to perform. If you want to perform only
# sequential analysis choose the first (remember to add -Dscotch in the ORDERINGSF
# variable below); for both parallel and sequential analysis choose the second
# line (remember to add -Dptscotch in the ORDERINGSF variable below)
LSCOTCH = -L$(SCOTCHDIR)/lib -lptesmumps -lptscotch -lptscotcherr -lscotch
#LSCOTCH = -L$(SCOTCHDIR)/lib -lptesmumps -lptscotch -lptscotcherr
LPORDDIR = $(topdir)/PORD/lib/
IPORD = -I$(topdir)/PORD/include/
LPORD = -L$(LPORDDIR) -lpord
LMETISDIR = /usr/local/lib
IMETIS = -I/usr/local/include
# You have to choose one among the following two lines depending on
# the type of analysis you want to perform. If you want to perform only
# sequential analysis choose the first (remember to add -Dmetis in the ORDERINGSF
# variable below); for both parallel and sequential analysis choose the second
# line (remember to add -Dparmetis in the ORDERINGSF variable below)
LMETIS = -L$(LMETISDIR) -lmetis
#LMETIS = -L$(LMETISDIR) -lparmetis -lmetis
# The following variables will be used in the compilation process.
# Please note that -Dptscotch and -Dparmetis imply -Dscotch and -Dmetis respectively.
# If you want to use Metis 4.X or an older version, you should use -Dmetis4 instead of -Dmetis
# or in addition with -Dparmetis (if you are using parmetis 3.X or older).
#ORDERINGSF = -Dscotch -Dmetis -Dpord -Dptscotch -Dparmetis
ORDERINGSF = -Dpord -Dmetis -Dptscotch
ORDERINGSC = $(ORDERINGSF)
LORDERINGS = $(LMETIS) $(LPORD) $(LSCOTCH)
IORDERINGSF = $(ISCOTCH)
IORDERINGSC = $(IMETIS) $(IPORD) $(ISCOTCH)
#End orderings
########################################################################
################################################################################
PLAT =
LIBEXT = .so
OUTC = -o
OUTF = -o
RM = /bin/rm -f
CC = mpiicc
FC = mpiifort
FL = mpiifort
AR = mpiicc -shared -o
#RANLIB = ranlib
RANLIB = echo
# Make this variable point to the path where the Intel MKL library is
# installed. It is set to the default install directory for Intel MKL.
MKLROOT=/opt/intel/mkl/lib/intel64
LAPACK = -L$(MKLROOT) -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core
SCALAP = -L$(MKLROOT) -lmkl_scalapack_lp64 -lmkl_blacs_intelmpi_lp64
INCPAR = -I/usr/local/include
LIBPAR = $(SCALAP) $(LAPACK)
INCSEQ = -I$(topdir)/libseq
LIBSEQ = $(LAPACK) -L$(topdir)/libseq -lmpiseq
LIBBLAS = -L$(MKLROOT) -lmkl_intel_lp64 -lmkl_intel_thread -lmkl_core
LIBOTHERS = -lpthread
#Preprocessor defs for calling Fortran from C (-DAdd_ or -DAdd__ or -DUPPER)
CDEFS = -DAdd_
#Begin Optimized options
OPTF = -O3 -g -fPIC -nofor_main -DBLR_MT -qopenmp # or -openmp for old compilers
OPTL = -O3 -g -nofor_main -qopenmp
OPTC = -O3 -g -fPIC -qopenmp
#End Optimized options
INCS = $(INCPAR)
LIBS = $(LIBPAR)
LIBSEQNEEDED =
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment