ParaMonte Fortran 2.0.0
Parallel Monte Carlo and Machine Learning Library
See the latest version documentation.
pm_matrixDet::setMatDet Interface Reference

Return the determinant of the input square matrix. More...

Detailed Description

Return the determinant of the input square matrix.

This generic interface uses the pivoted LU (LUP) factorization to compute the determinant.
If the matrix is already known to be positive-definite, use the other more appropriate faster generic interfaces of this module (e.g., setMatDetSqrt).

Parameters
[in,out]mat: The input/output contiguous square matrix of shape (1:ndim, 1:ndim) of,
  1. type complex of kind any supported by the processor (e.g., CK, CK32, CK64, or CK128),
  2. type real of kind any supported by the processor (e.g., RK, RK32, RK64, or RK128),
On input, it represents the general matrix whose determinant must be computed.
On output, it is completely overwritten by its LUP factorization.
[out]det: The output scalar of the same type and kind as the input mat representing its determinant.
[out]info: The output scalar of type integer of default kind IK that is non-zero if and only if the Pivoted LU-factorization of the input mat fails.
See the corresponding info argument of setMatLUP for possible meanings of the output non-zero error code.


Possible calling interfaces ⛓

use pm_matrixDet, only: setMatDet
call setMatDet(mat, det, info)
pm_matrixDet::setMatDet
Return the determinant of the input square matrix.
Definition: pm_matrixDet.F90:395
pm_matrixDet
This module contains procedures and generic interfaces relevant to the computation of the determinant...
Definition: pm_matrixDet.F90:95
Warning
The condition size(mat, 1) == size(mat, 2) must hold for the corresponding input arguments.
This condition is verified only if the library is built with the preprocessor macro CHECK_ENABLED=1.
The pure procedure(s) documented herein become impure when the ParaMonte library is compiled with preprocessor macro CHECK_ENABLED=1.
By default, these procedures are pure in release build and impure in debug and testing builds.
See also
getMatDet
setMatDet
getMatDetSqrtLog
setMatDetSqrtLog
getMatDetSqrt
setMatDetSqrt


Example usage ⛓

1program example
2
3 use pm_kind, only: SK, IK, LK
4 use pm_matrixDet, only: getMatDet
5 use pm_matrixDet, only: setMatDet
6 use pm_matrixInv, only: getMatInv
7 use pm_distUnif, only: getUnifRand
8 use pm_io, only: display_type
9
10 implicit none
11
12 integer(IK) :: info, ndim, itry, ntry = 10
13 type(display_type) :: disp
14
15 disp = display_type(file = "main.out.F90")
16
17 call disp%skip()
18 call disp%show("!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
19 call disp%show("! Compute the determinant of the square matrix.")
20 call disp%show("!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%")
21 call disp%skip()
22
23 block
24 use pm_kind, only: TKG => RKH
25 real(TKG), allocatable :: mat(:,:), lup(:,:)
26 real(TKG) :: det
27 do itry = 1, ntry
28 call disp%skip()
29 call disp%show("ndim = getUnifRand(1, 9)")
30 ndim = getUnifRand(1, 9)
31 call disp%show("mat = getUnifRand(1._TKG, 2._TKG, ndim, ndim)")
32 mat = getUnifRand(1._TKG, 2._TKG, ndim, ndim)
33 call disp%show("lup = mat")
34 lup = mat
35 call disp%show("lup ! matrix.")
36 call disp%show( lup )
37 call disp%show("call setMatDet(lup, det, info)")
38 call setMatDet(lup, det, info)
39 call disp%show("lup ! lup factorization.")
40 call disp%show( lup )
41 call disp%show("if (info /= 0) error stop 'determinant computation failed.'")
42 if (info /= 0) error stop 'determinant computation failed.'
43 call disp%show("det * getMatDet(getMatInv(mat)) ! must be one.")
44 call disp%show( det * getMatDet(getMatInv(mat)) )
45 call disp%skip()
46 end do
47 end block
48
49 block
50 use pm_kind, only: TKG => CKH
51 complex(TKG), allocatable :: mat(:,:), lup(:,:)
52 complex(TKG) :: det
53 do itry = 1, ntry
54 call disp%skip()
55 call disp%show("ndim = getUnifRand(1, 9)")
56 ndim = getUnifRand(1, 9)
57 call disp%show("mat = getUnifRand((1._TKG, 1._TKG), (2._TKG, 2._TKG), ndim, ndim)")
58 mat = getUnifRand((1._TKG, 1._TKG), (2._TKG, 2._TKG), ndim, ndim)
59 call disp%show("lup = mat")
60 lup = mat
61 call disp%show("lup ! matrix.")
62 call disp%show( lup )
63 call disp%show("call setMatDet(lup, det, info)")
64 call setMatDet(lup, det, info)
65 call disp%show("if (info /= 0) error stop 'determinant computation failed.'")
66 if (info /= 0) error stop 'determinant computation failed.'
67 call disp%show("lup ! lup factorization.")
68 call disp%show( lup )
69 call disp%show("det * getMatDet(getMatInv(mat)) ! must be one.")
70 call disp%show( det * getMatDet(getMatInv(mat)) )
71 call disp%skip()
72 end do
73 end block
74
75end program example
pm_distUnif::getUnifRand
Generate and return a scalar or a contiguous array of rank 1 of length s1 of randomly uniformly distr...
Definition: pm_distUnif.F90:4150
pm_io::show
This is a generic method of the derived type display_type with pass attribute.
Definition: pm_io.F90:11726
pm_io::skip
This is a generic method of the derived type display_type with pass attribute.
Definition: pm_io.F90:11508
pm_matrixDet::getMatDet
Generate and return the determinant of the input general square matrix.
Definition: pm_matrixDet.F90:216
pm_matrixInv::getMatInv
Generate and return the full inverse of an input matrix of general or triangular form directly or thr...
Definition: pm_matrixInv.F90:362
pm_distUnif
This module contains classes and procedures for computing various statistical quantities related to t...
Definition: pm_distUnif.F90:274
pm_io
This module contains classes and procedures for input/output (IO) or generic display operations on st...
Definition: pm_io.F90:252
pm_io::disp
type(display_type) disp
This is a scalar module variable an object of type display_type for general display.
Definition: pm_io.F90:11393
pm_kind
This module defines the relevant Fortran kind type-parameters frequently used in the ParaMonte librar...
Definition: pm_kind.F90:268
pm_kind::LK
integer, parameter LK
The default logical kind in the ParaMonte library: kind(.true.) in Fortran, kind(....
Definition: pm_kind.F90:541
pm_kind::CKH
integer, parameter CKH
The scalar integer constant of intrinsic default kind, representing the highest-precision complex kin...
Definition: pm_kind.F90:843
pm_kind::IK
integer, parameter IK
The default integer kind in the ParaMonte library: int32 in Fortran, c_int32_t in C-Fortran Interoper...
Definition: pm_kind.F90:540
pm_kind::SK
integer, parameter SK
The default character kind in the ParaMonte library: kind("a") in Fortran, c_char in C-Fortran Intero...
Definition: pm_kind.F90:539
pm_kind::RKH
integer, parameter RKH
The scalar integer constant of intrinsic default kind, representing the highest-precision real kind t...
Definition: pm_kind.F90:858
pm_matrixInv
This module contains abstract and concrete derived types and procedures related to the inversion of s...
Definition: pm_matrixInv.F90:212
pm_io::display_type
Generate and return an object of type display_type.
Definition: pm_io.F90:10282

Example Unix compile command via Intel ifort compiler ⛓
1#!/usr/bin/env sh
2rm main.exe
3ifort -fpp -standard-semantics -O3 -Wl,-rpath,../../../lib -I../../../inc main.F90 ../../../lib/libparamonte* -o main.exe
4./main.exe

Example Windows Batch compile command via Intel ifort compiler ⛓
1del main.exe
2set PATH=..\..\..\lib;%PATH%
3ifort /fpp /standard-semantics /O3 /I:..\..\..\include main.F90 ..\..\..\lib\libparamonte*.lib /exe:main.exe
4main.exe

Example Unix / MinGW compile command via GNU gfortran compiler ⛓
1#!/usr/bin/env sh
2rm main.exe
3gfortran -cpp -ffree-line-length-none -O3 -Wl,-rpath,../../../lib -I../../../inc main.F90 ../../../lib/libparamonte* -o main.exe
4./main.exe

Example output ⛓
1
2!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
3! Compute the determinant of the square matrix.
4!%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5
6
7ndim = getUnifRand(1, 9)
8mat = getUnifRand(1._TKG, 2._TKG, ndim, ndim)
9lup = mat
10lup ! matrix.
11+1.22681961522596335290125423650711481, +1.49389440623510343646024504465687758, +1.32816211615524444108700765863733518, +1.06565510053171045620477029599708936
12+1.08993507772181448640892792882840410, +1.98278276447953944717111619675542958, +1.36237815056562168866369184977689514, +1.10163250880377549467807011116490388
13+1.46361841282694968964788303349776360, +1.22595238046117168424607019764763230, +1.80516686904033091579129089604725386, +1.46971273769540524323863180458792312
14+1.23565139203064559447729281312154751, +1.03569317744681506725615616174661765, +1.51731698095606546633919268398627348, +1.82346206605345182349228077226521288
15call setMatDet(lup, det, info)
16lup ! lup factorization.
17+1.22681961522596335290125423650711481, +1.49389440623510343646024504465687758, +1.32816211615524444108700765863733518, +1.06565510053171045620477029599708936
18+0.888423256520122877698030516231730704, +0.655572231194953493106323558279227663, +0.182408038144321718933100721477672710, +0.154879734062114362405245429901751657
19+1.19301843128532897221541042534299524, -0.848558181845762177374883613214883368, +0.375428817933996381534464200024336886, +0.329791026908357419948286332542830195
20+1.00719892044035793263385825839729971, -0.715337888717884104796703270243897759, +0.825927306253209759872293465272365574, +0.588543326724432737860524689244398956
21if (info /= 0) error stop 'determinant computation failed.'
22det * getMatDet(getMatInv(mat)) ! must be one.
23+1.00000000000000000000000000000000058
24
25
26ndim = getUnifRand(1, 9)
27mat = getUnifRand(1._TKG, 2._TKG, ndim, ndim)
28lup = mat
29lup ! matrix.
30+1.59328162193872098382069302862157234, +1.18733630675410624543669613201305982, +1.32909774562815774929611656690346339, +1.48892223907358986807686465837164042, +1.10960273331059775687423710275002921, +1.48065246788704503130817472568213717
31+1.98841044681667655006219383417318236, +1.64116722142214775558984779237374260, +1.07562780866181277907113045119216984, +1.59048290712167858962288580751578692, +1.16492763031945721898800327343400299, +1.08912884890196322679561281211879282
32+1.68353011220788192065870254079650664, +1.45548318038507860834754545056038709, +1.70920201720020400849303912405633654, +1.11601488202285799410451979614068149, +1.52496308314667171331349725177213026, +1.65276771057217350651543134232731021
33+1.38457217053310067223707515935775642, +1.73263683399827784709398378953650029, +1.65806590606262203698521743407553883, +1.13586072621758299531243764468542558, +1.69384629154809754803861378658557120, +1.26034574855107744092785743175318679
34+1.77535708930620956680109014828533463, +1.74619513207692222844260722976638700, +1.03746581979444175831438737054413905, +1.68796333892132069122798521298046839, +1.51642489226190801391944705896069847, +1.30698582681807607335914905356834450
35+1.95438232704696730894746448922031358, +1.45852987835455523445246818483477392, +1.39280190854724807950022583383964597, +1.81937202474781636976307840918066334, +1.23828665971667953345668559301024822, +1.19666504773617000778973716837477924
36call setMatDet(lup, det, info)
37lup ! lup factorization.
38+1.95438232704696730894746448922031358, +1.45852987835455523445246818483477392, +1.39280190854724807950022583383964597, +1.81937202474781636976307840918066334, +1.23828665971667953345668559301024822, +1.19666504773617000778973716837477924
39+0.708444888889863013107533493547882938, +0.699348796384839547472542408539436792, +0.671342512716277725593687363317525160, -0.153064085704208872914578895294300668, +0.816588436491314906120652442263938682, +0.412574511769243861151319310552755742
40+0.908398047166512550577595933341215138, +0.602373866979268901746022299938787569, -0.632151899472153222388305205274817090, +0.127451149772310887483668272211809554, -0.100323825476896014457789120179791501, -0.285864697294055958297894489761677418E-1
41+0.861412881660500031585901753013926937, +0.284674479906001890962063035165868150, -0.503534721948584144348808302810941315, -0.343466098374895807046575491776751709, +0.175308585256624331532038358440659359, +0.490101308772289454028898342268629032
42+1.01741118884406054313913559366101934, +0.224841459274248643369613275650012188, +0.778879360452158808382811472757641617, +0.947460220084838665980081842678249339, -0.366479761704258942554781181635627070, -0.663221497861821832257364651778886608
43+0.815235381475301056482608326003965489, -0.244349461703181789478091171453322025E-2, -0.308907989318987275394511096138288968, -0.130151009917428227491414509869612005, -0.256299892584185462541109657405793337, +0.391069897547542301643836905067167578
44if (info /= 0) error stop 'determinant computation failed.'
45det * getMatDet(getMatInv(mat)) ! must be one.
46+0.999999999999999999999999999999999904
47
48
49ndim = getUnifRand(1, 9)
50mat = getUnifRand(1._TKG, 2._TKG, ndim, ndim)
51lup = mat
52lup ! matrix.
53+1.17363277684381070677871923815354808, +1.25677899353063919177612772361024309, +1.74118970097446322153261383574485626, +1.63655931224816814599699439040180505, +1.06375935994425497208594154401612328, +1.36613617559510647640085583683877225, +1.82590776527198428033190558974227115
54+1.07475673975466838485137859746362877, +1.89231888297207933158310641125604980, +1.87406046940062171042101052538462896, +1.48900136014797932608069907405266157, +1.95023697561481857435021767758719590, +1.27136265095748367631408656083118140, +1.12842797382092661473389986665318035
55+1.61177441169325830074710715796367873, +1.33101026093421513158257973629327218, +1.04207128754977530103835847954425757, +1.01460548393713286798908368291974037, +1.58832199174890199276086787513457994, +1.92324141899487452293389976702166559, +1.15940715323329504374593846932792284
56+1.99408096801560871926642891476281282, +1.02853878465690343178251190272245741, +1.45297500941849687311821646874086765, +1.65286784337162071157222252296652406, +1.83716474191494444800470888811838657, +1.86657632863177306325709140120948852, +1.86003830613657917800662769743448384
57+1.64136869698913701129374800558232341, +1.67442664962182288000914160678089956, +1.47302235958167312296031882757295585, +1.63968647089504924040635848444992436, +1.10462235501625686993848466090637874, +1.24495378759056363454718425387997981, +1.85642745180617708409825481529415960
58+1.43339453420056360187632992619083814, +1.83449147123770681909288565625097818, +1.11823583187752170441096869372167218, +1.17890082626088246092548477402710321, +1.58193193575528416106288762705123723, +1.87259691888575558834578200439895973, +1.99955144776873937702063001445206429
59+1.13227982543284624289826774945381360, +1.10396368158656528408300914211214680, +1.73651720845681928201381027287740025, +1.95983675068612246353675340696470503, +1.44289946591686500663318069408641211, +1.01214483760732927901719986995939562, +1.38013641844720582302351235386623190
60call setMatDet(lup, det, info)
61lup ! lup factorization.
62+1.99408096801560871926642891476281282, +1.02853878465690343178251190272245741, +1.45297500941849687311821646874086765, +1.65286784337162071157222252296652406, +1.83716474191494444800470888811838657, +1.86657632863177306325709140120948852, +1.86003830613657917800662769743448384
63+0.538973470482596622053795578662288660, +1.33796376467959598637363982147375170, +1.09094548604985093705838726831829076, +0.598149442356891973728262917080783120, +0.960053918816657021747399809504191515, +0.265327529194153165154613390299600000, +0.125916672731924037801743234614630825
64+0.718824640118296174354205884820962303, +0.818521755535925248724075691844833352, -0.819161020804124542240933076333862169, -0.498819638044625069221667585176990891, -0.524492367728703016925161227497086604, +0.313679506215683078213168334162262770, +0.559444545738090610689751746142080261
65+0.567820386230166002916209493509110508, +0.388604239782597697713680559083419972, -0.595172772865556277345197641651193462, +0.491977417017129885333596478630637900, -0.285524727525685943208965505194217080, +0.358508444851483394063646726430312014E-1, +0.608003157722337106495737902746787481
66+0.823120386441745105838302311642566464, +0.618712875174872739454844290614633453, +0.485781183293367619936014009110711360, +0.307761454975489092827529134345333105, -0.658921089141504201377461933758641919, -0.619037909571777740032214960385859068, -0.211401834021823600070863998387474421
67+0.588558235933488968790001898295851954, +0.486877177125790712011584729896440287, -0.433214990602499935049947414474026826, +0.317959458929384150636658549632220162, +0.943026562785607165782209833609399815, +0.846626127272797403472906191952751171, +0.918257900342702316874588082893645171
68+0.808279321424546134944171311707957843, +0.373450793813209863229950949032913065, +0.658909312301311201515606793997823761, -0.439198806356893307141575009198337098, +0.530600583896599929484081890733424638E-1, +0.185849547094092597184754304632660151, -0.652076809068013370540282342214347294
69if (info /= 0) error stop 'determinant computation failed.'
70det * getMatDet(getMatInv(mat)) ! must be one.
71+1.00000000000000000000000000000000058
72
73
74ndim = getUnifRand(1, 9)
75mat = getUnifRand(1._TKG, 2._TKG, ndim, ndim)
76lup = mat
77lup ! matrix.
78+1.77638937700595080313023565853481253, +1.86229296969493781867869565331458926
79+1.61598441566278777443996236541139867, +1.91438502319569646643868747940688137
80call setMatDet(lup, det, info)
81lup ! lup factorization.
82+1.77638937700595080313023565853481253, +1.86229296969493781867869565331458926
83+0.909701688481429327295828797256215967, +0.220253964217116236103654915231444138
84if (info /= 0) error stop 'determinant computation failed.'
85det * getMatDet(getMatInv(mat)) ! must be one.
86+0.999999999999999999999999999999997496
87
88
89ndim = getUnifRand(1, 9)
90mat = getUnifRand(1._TKG, 2._TKG, ndim, ndim)
91lup = mat
92lup ! matrix.
93+1.01959633253147223023283213554488262, +1.27704907596072836597060282544927801, +1.11127635557559039963876934377941745, +1.92495033074274561672982844448933326, +1.69581447963738752174885458853266906, +1.25660284265549331561479208986070078
94+1.04267294523094680302141578511399540, +1.51835286484779846138857431442340225, +1.96830024417052613679742734657786734, +1.01144992336511927560705277635191415, +1.85996052917931746562411739680279574, +1.38889572867042052050992622262886116
95+1.06198232437233768572753025029311722, +1.16700924310250318418567742542558252, +1.51689653367006462547562981001818165, +1.64999469656067081476275854548856236, +1.87398374123169314679083226993803036, +1.95401039519485542157283151688958316
96+1.60315411027367955931023264985121368, +1.92213934128251413206034648138452932, +1.14555742271980621286054378356019321, +1.56110801476691693498966337391260484, +1.28146996168514523029699987351562907, +1.51696575586273012934798007794498755
97+1.39681820375384168436950129778163886, +1.24999042077118459619659776562888118, +1.70655039419792069989314793526270914, +1.52945545471803683765429906419725450, +1.51874955954022517045440639470317956, +1.63851442661900815294093749295304621
98+1.43532342185748402129430363889930608, +1.42890641425731746470775878905873854, +1.18215287860594923376425899337702082, +1.35741106478447038058500507761907342, +1.92410205110747093490891674873015548, +1.23834050555342925806696103950768568
99call setMatDet(lup, det, info)
100lup ! lup factorization.
101+1.60315411027367955931023264985121368, +1.92213934128251413206034648138452932, +1.14555742271980621286054378356019321, +1.56110801476691693498966337391260484, +1.28146996168514523029699987351562907, +1.51696575586273012934798007794498755
102+0.871293779432961929520638055722712049, -0.424757630491641006556372118246069585, +0.708433337798897534424749828379743764, +0.169271752428681638924440854065768315, +0.402212753393762067351169483671874505, +0.316791599922990192609116327084416413
103+0.650388467676977588765388915983692381, -0.631455646126897202863544945682161737, +1.67058713842944097132278743189954429, +0.103010877563501013445938089223015163, +1.28048675849957706985324619214194510, +0.602318539613356679851513704334438178
104+0.635993960903367975826621593182647372, -0.128496956787091088158306689750667410, +0.283577114868022522836423376587246285, +0.924634438630763405465082332285164225, +0.569373697116603594708765955473570380, +0.161724785860325577807114637949561263
105+0.895312192795024163288383656691181213, +0.687470579176941090241368424100067771, -0.197838064800088727833541663771644906, -0.147364024552320189733815515900207205, +0.837511156967257559625031398877114842, -0.194608386815878998429996331292904647
106+0.662433085856632517205578916043707607, +0.250211990463109741831152173065903577, +0.347651951933937266235757429924850686, +0.581526406283990142973890907683935166, +0.176938091267395929830116267644830528, +0.600846218617987299613974648440068747
107if (info /= 0) error stop 'determinant computation failed.'
108det * getMatDet(getMatInv(mat)) ! must be one.
109+0.999999999999999999999999999999999519
110
111
112ndim = getUnifRand(1, 9)
113mat = getUnifRand(1._TKG, 2._TKG, ndim, ndim)
114lup = mat
115lup ! matrix.
116+1.15729081006135263443418229032208201
117call setMatDet(lup, det, info)
118lup ! lup factorization.
119+1.15729081006135263443418229032208201
120if (info /= 0) error stop 'determinant computation failed.'
121det * getMatDet(getMatInv(mat)) ! must be one.
122+1.00000000000000000000000000000000000
123
124
125ndim = getUnifRand(1, 9)
126mat = getUnifRand(1._TKG, 2._TKG, ndim, ndim)
127lup = mat
128lup ! matrix.
129+1.85902527791433253210413113001680354, +1.50387603608432757378205444718840080, +1.72724073160619332075851067836863267, +1.04081425901443749468240621275457374
130+1.01596366615204003549062132524248633, +1.41663126396201809613053260906184783, +1.83303675864318472377648880928495127, +1.10760881753405646855614184757480868
131+1.93206694396950054124477504633850345, +1.89391527923819678144896155624524117, +1.36494901979551560248374167238466036, +1.13620086030605011867408586239194329
132+1.89265398890716707679463316001284014, +1.24152061938194161809613360456364935, +1.26847951394685978455470994375815420, +1.32210206155581277912998932271824390
133call setMatDet(lup, det, info)
134lup ! lup factorization.
135+1.89265398890716707679463316001284014, +1.24152061938194161809613360456364935, +1.26847951394685978455470994375815420, +1.32210206155581277912998932271824390
136+0.536793133930764199113197840364514637, +0.750191519844322186348065338167481373, +1.15212566502467734562341363120696695, +0.397913508535187606015613991505855418
137+1.02082417351155179022559046739081444, +0.835174755361263671051372207710234654, -0.892171802078040587599684537451384575, -0.545760201125398643961417325708849207
138+0.982231981550810552418037793173420407, +0.379122891220677124187737093975339385, -0.498809442174753083776876321794882417E-1, -0.435877822680018224242279137696985598
139if (info /= 0) error stop 'determinant computation failed.'
140det * getMatDet(getMatInv(mat)) ! must be one.
141+0.999999999999999999999999999999999807
142
143
144ndim = getUnifRand(1, 9)
145mat = getUnifRand(1._TKG, 2._TKG, ndim, ndim)
146lup = mat
147lup ! matrix.
148+1.99444185810286561802824121873969921
149call setMatDet(lup, det, info)
150lup ! lup factorization.
151+1.99444185810286561802824121873969921
152if (info /= 0) error stop 'determinant computation failed.'
153det * getMatDet(getMatInv(mat)) ! must be one.
154+0.999999999999999999999999999999999904
155
156
157ndim = getUnifRand(1, 9)
158mat = getUnifRand(1._TKG, 2._TKG, ndim, ndim)
159lup = mat
160lup ! matrix.
161+1.57537326765791815057742253913024572, +1.16073459667968365970423691861761286, +1.28579661393606436753635689288501616, +1.82480558673328573557075263388868656, +1.82219152941920093835065180658164903, +1.27918023476985170729370892640023442
162+1.74636485390579394289765619952449725, +1.96911266594932472792626094350099534, +1.25259665732114185128659578254529626, +1.48839591246553429117844743918531562, +1.81423871753879548851432914199109463, +1.65584915039132523342331928201722439
163+1.45333366921900867742181424969508814, +1.26965861718364082545398620301085412, +1.43520018628394539298543261042711076, +1.14175410067957989413301822127238313, +1.32612171956029192484234150374537017, +1.71573328747799300801970704065897634
164+1.81009543822604882217179967465489059, +1.19257060619131848492761197630444638, +1.36093260543051088055477281520437505, +1.87385370324921374504400943735777264, +1.80502234857280233941921448223808940, +1.53809007011564575452157119123259760
165+1.55899092135997118934725593494458270, +1.90588586987689845287489208020478315, +1.17992898336772597917462484388032011, +1.15378465989213925343503149588338396, +1.73674319698403809523605390930525122, +1.86750114947528356733858563125323568
166+1.75211077357608780083191861096418366, +1.32988164733060733868520190733015430, +1.58773897593931522828073356747225597, +1.12353033734907799579399377687123361, +1.44060358421997246321248276133804574, +1.96536275529050169182576754924350683
167call setMatDet(lup, det, info)
168lup ! lup factorization.
169+1.81009543822604882217179967465489059, +1.19257060619131848492761197630444638, +1.36093260543051088055477281520437505, +1.87385370324921374504400943735777264, +1.80502234857280233941921448223808940, +1.53809007011564575452157119123259760
170+0.861275537431237296199152666455920015, +0.878753980104774176290308010910585682, +0.779102221786870726534824428486227967E-2, -0.460119695441341564117489131999192819, +0.182121603641603619580333934273046847, +0.542781797718781314863813004253984544
171+0.802904442786354933520059839656547585, +0.355205650504135813841902517246410263, +0.339733935935988587084769549931509485, -0.199334247061871174734747398146980145, -0.187829366129843866384963246878026407, +0.287994775236060885696659066601336461
172+0.967965962774433684332817918201932200, +0.199730409541031208591864596482384755, +0.791344060032286998519002488049684880, -0.440654399033171327796463964070171410, -0.194334180540794860833931000644091784, +0.140230934284034720367015612684628705
173+0.964791588899471023685145389065740427, +0.931467275777918541558660036184557998, -0.199205167748718557633314039811128940, -0.157483161441969408157241034463372023, -0.164892917659993536313303296137927891, -0.254216636314337051250932043991476098
174+0.870326080265598661338643060640967401, +0.139753899650709924989653585587969502, +0.295091651449727415652404798896947449, -0.719537199194857546645762181834927015, -0.857389799080707635752733884586240813, -0.337361770810573137139799143567222758
175if (info /= 0) error stop 'determinant computation failed.'
176det * getMatDet(getMatInv(mat)) ! must be one.
177+0.999999999999999999999999999999996341
178
179
180ndim = getUnifRand(1, 9)
181mat = getUnifRand(1._TKG, 2._TKG, ndim, ndim)
182lup = mat
183lup ! matrix.
184+1.91327476778916809899276649693662397, +1.05699651077892179845177542009360730, +1.15893098215783426178403000559239876, +1.39169937869606397442307122611475654, +1.17157781039609795028158266983655841
185+1.04526698097198447633802234745675026, +1.45162822624046632003319123959422719, +1.29783047809118177600938011044288673, +1.29003404382941263188653882240088774, +1.29052926567888388494916377388434287
186+1.11664561834198296890989507508251502, +1.27373951939911466110322227345356092, +1.67363662108513564393953138462898045, +1.50463621626174719757575837737359168, +1.79435931282090638625396300025823377
187+1.70827529713689386662544727746797325, +1.50679833928039720545444946362045188, +1.97262127004174919498740179530134730, +1.87231417354175504805168458559308083, +1.17877090983864258633690519452234777
188+1.93857156775839813324985188567196767, +1.52425554285071403193168298554618103, +1.64063952014458155931587196602574888, +1.07920138434343907528297653728803998, +1.89657001831913679529319734701120038
189call setMatDet(lup, det, info)
190lup ! lup factorization.
191+1.93857156775839813324985188567196767, +1.52425554285071403193168298554618103, +1.64063952014458155931587196602574888, +1.07920138434343907528297653728803998, +1.89657001831913679529319734701120038
192+0.539194424573472778723970899456958214, +0.629758135910149187953056153316725258, +0.413206796094325620911361274111735602, +0.708134674399456764828053541877754786, +0.267909285987996193983794548639898164
193+0.576014647544417706549204384364010390, +0.628409507634846407172832489499970324, +0.468941146856604364402132292410184095, +0.438001849151188568801235918127246873, +0.533550459676979028248202713572487254
194+0.986950804195234774884976096770297083, -0.710381807755577670541014768273969076, -0.355620149724164776131694739095617901, +0.985388977838596325148694405813198431, -0.320184317146162245102262942106244550
195+0.881203111377621413314939989827776996, +0.259813415747767498451352297276659867, +0.894628150994998737387033953467265066, +0.350609648191331565085730333282688450, -0.927168468365129699559186329576971140
196if (info /= 0) error stop 'determinant computation failed.'
197det * getMatDet(getMatInv(mat)) ! must be one.
198+1.00000000000000000000000000000000116
199
200
201ndim = getUnifRand(1, 9)
202mat = getUnifRand((1._TKG, 1._TKG), (2._TKG, 2._TKG), ndim, ndim)
203lup = mat
204lup ! matrix.
205(+1.40310275418108527838196306746443629, +1.62827454247677685020391530843633079), (+1.58161647366706855761460184383796564, +1.13058876299620673278208131788685739), (+1.12312162568685723108939268446761982, +1.45695096841963363123810661956504126)
206(+1.03130043855084853554299943618257614, +1.58091837913118572860984134640684382), (+1.66712421464060512384018440107098967, +1.67660254660155429514603764351699810), (+1.27094504820022294912476419073482093, +1.93545308634703625640888400638991788)
207(+1.41188814471889051809659695483824451, +1.74264186719530909350689004054986436), (+1.83796448618022828261764270558654633, +1.56860737362439611913087313631421109), (+1.59807761443075622412009625432521532, +1.66392812607102838021872147605503325)
208call setMatDet(lup, det, info)
209if (info /= 0) error stop 'determinant computation failed.'
210lup ! lup factorization.
211(+1.41188814471889051809659695483824451, +1.74264186719530909350689004054986436), (+1.83796448618022828261764270558654633, +1.56860737362439611913087313631421109), (+1.59807761443075622412009625432521532, +1.66392812607102838021872147605503325)
212(+0.957913390222172191085207955931778263, -0.290571435313994198388860516851392882E-1), (-0.224573567997602219214235774123869182, -0.318595246319282398926647878546317999), (-0.456047318175848129592214831840194938, -0.905124932942602384525097428049519684E-1)
213(+0.837149874800140473072051839044856317, +0.864558277110919806696660950406759814E-1), (-0.819243123764001815623351505190456950, +0.251439321686170164221288759624457370), (-0.319401172694765101725419066822926358, +0.444849231620938405776086710498487299)
214det * getMatDet(getMatInv(mat)) ! must be one.
215(+1.00000000000000000000000000000000019, +0.674075480535532548569592279904724561E-33)
216
217
218ndim = getUnifRand(1, 9)
219mat = getUnifRand((1._TKG, 1._TKG), (2._TKG, 2._TKG), ndim, ndim)
220lup = mat
221lup ! matrix.
222(+1.17478964862654337844504630829445806, +1.95005960804935313594515217923289927), (+1.23848514717393704814546008451309138, +1.70339628427064880018539415514400358), (+1.96929293866405597591371129173584345, +1.78996798191968859611770830874684921), (+1.81452806436086219325592042198445311, +1.75408524365533189880742869642759736), (+1.96395150405113808777899305467889874, +1.49696156014982975100154964797080116), (+1.89285763570219730697587598017624015, +1.86115093214261426928060858847651195)
223(+1.33346032981356407643126057954632607, +1.69934300935154224609712858546336027), (+1.65034504082428154520083794913119820, +1.73796447690564115145293557366312024), (+1.94667861311619507037828081967156472, +1.71156085840423247302724591132211422), (+1.72102149397892462724095375398767492, +1.45874161805835004201165333939400525), (+1.43114815343755291624791363887857761, +1.60033103003835476384936287003825191), (+1.12090251243699000187232884751407363, +1.62716367742997737019163889072285406)
224(+1.66486298213635729916910861018930853, +1.99997887593218046450539346659633876), (+1.31830246843478302754021862520052707, +1.10061309352676131814173185522399865), (+1.62246638002534667647887205226653024, +1.48496138847530216653591993948567534), (+1.65766260940260215389090425083885766, +1.39900766936676781317614677302415674), (+1.57674555688435476898895159221658293, +1.64590148115008923029595481897037336), (+1.56028999918540256204629768552070355, +1.32400826929251875076405698806705914)
225(+1.64892973172933793888090841797846207, +1.40603324379709346305058560346192253), (+1.42066086588531891883383809796869882, +1.47640949618308683351003792600131735), (+1.07761517996610412352530263991599111, +1.33503068622637675615721962742048890), (+1.70847175775820868118060143388660328, +1.39035947399091001545199132003098037), (+1.38550095081653627748231573047586025, +1.67478726962974853809038468759219836), (+1.72527905067173279634953175789381949, +1.52114896389685394968719042121304445)
226(+1.94544656799996472638550929776484245, +1.54619826325314243670590837948219778), (+1.61316857967918672608541360390753837, +1.75750967689692407859909860893247937), (+1.14522661263793404264342752914069787, +1.79758014296708267613948570609994579), (+1.05519128254675643288497601594449430, +1.50668636018962871555087903178198507), (+1.83281409053689958775271166896902478, +1.34212053877990546969274759106917815), (+1.89573293925899152709525068349984220, +1.36050519152509157041156709093899112)
227(+1.10531028693704016278525077441249077, +1.27425580516552729680450436918342601), (+1.34774234566138363133197734945840616, +1.85129402067119323949108474932451304), (+1.38753125463829557246943219149869404, +1.43932211333487291238840064567227346), (+1.57545265144910723558047029665759860, +1.20748041230342186303029078755174236), (+1.71706563538432915495086953301515443, +1.53115863986605259310110922669139823), (+1.58911194466223820640837964138498998, +1.87139728676565249699604635964367263)
228call setMatDet(lup, det, info)
229if (info /= 0) error stop 'determinant computation failed.'
230lup ! lup factorization.
231(+1.94544656799996472638550929776484245, +1.54619826325314243670590837948219778), (+1.61316857967918672608541360390753837, +1.75750967689692407859909860893247937), (+1.14522661263793404264342752914069787, +1.79758014296708267613948570609994579), (+1.05519128254675643288497601594449430, +1.50668636018962871555087903178198507), (+1.83281409053689958775271166896902478, +1.34212053877990546969274759106917815), (+1.89573293925899152709525068349984220, +1.36050519152509157041156709093899112)
232(+1.02522461947027492395497896075944424, +0.213204698954695248721173496572413720), (+0.391516467645573563458794302053607720E-1, -1.04516421757888782349144837132257477), (+0.831604395104623712086664195945453434, -0.602129724747823681387700241963554402), (+0.897087140128627762560163296126658548, -0.370656020654697848030332512790598518), (-0.161541662146129876382660380318101545E-1, -0.120838113916622952526881954575641839), (-0.931959822982556825916908844090432370E-1, -0.474994318589339621209756024460398489)
233(+0.845551661931926760371040229551491280, +0.201471736427353765052051559679792695), (+0.927724355388349687301416073252351758E-1, +0.609266935622690586102275805248708907), (+0.896484242899331758450059618427462656, -0.489924833568967213328251440484526169), (+0.823304051849557899229083678023056727, -0.540009628837280839963708345031639999), (+0.796845025860500738511195440223065484E-1, +0.117291186072123554595674954558290082), (-0.488686596242428582631496132763116113, +0.195697254943728795272725251731742681)
234(+0.871495154412156274611051069734258688, +0.300850974651334137321175811435525757E-1), (+0.101581951138185219604535802085638255, +0.609378063566382592224677456185669392E-1), (+0.130650207119079409158136593647387562, -0.213620939968707325080653813960911952), (+0.728284166797886042825236043843480188, +0.274956873926048342083419615647353168), (-0.212599108854625446261160636041544551, +0.464952864322709334922166669975940664), (+0.116651802715271676485459670498770749, +0.202410403138599891271537259559597229)
235(+0.858339670524192062833764103041965260, +0.320181654150464525715822292690172832), (+0.322229440799220959449272927048514099, +0.386489601814778329165466832826472938), (+1.02743612281957078672903907148307412, +0.285935898506950192344886321704384566), (+0.478573802486297889010731627868472219E-1, +0.277308802149542173110299013951604730), (+0.869776491582655717884774090379702730, -0.303275856140045973943538554547077943), (+1.15633244502703069097558073806433795, +0.172104086544313320069403318860089681)
236(+0.667246367636440377738457616421964298, +0.124681209113577564887655634474740143), (-0.438642350895150868325382758811460107, +0.485726036032645052545228284296089916), (+1.05807117132245746769261328144616130, -0.586414191424964358407277259686671069E-1), (+0.562379911152764827170893543785593506, -0.861188656456927860465542345253323068E-1), (+0.611985866394395568665602670755225845, +0.170549419736349986320982089165099965), (-0.335279643899635420174665754872184665E-1, -0.778836638913019300975479532832483049E-1)
237det * getMatDet(getMatInv(mat)) ! must be one.
238(+0.999999999999999999999999999999999519, +0.883520361987644447589429881160835407E-32)
239
240
241ndim = getUnifRand(1, 9)
242mat = getUnifRand((1._TKG, 1._TKG), (2._TKG, 2._TKG), ndim, ndim)
243lup = mat
244lup ! matrix.
245(+1.20600942850730226149795722163534822, +1.20335204502830191110186568767425614)
246call setMatDet(lup, det, info)
247if (info /= 0) error stop 'determinant computation failed.'
248lup ! lup factorization.
249(+1.20600942850730226149795722163534822, +1.20335204502830191110186568767425614)
250det * getMatDet(getMatInv(mat)) ! must be one.
251(+1.00000000000000000000000000000000000, -0.481482486096808963263994485646231830E-34)
252
253
254ndim = getUnifRand(1, 9)
255mat = getUnifRand((1._TKG, 1._TKG), (2._TKG, 2._TKG), ndim, ndim)
256lup = mat
257lup ! matrix.
258(+1.73566573003082594870842782602118366, +1.05214673495039012763106129783708815), (+1.97202609449081485340655330937006239, +1.29845572778754108943784999222755612), (+1.58180341211291633299345833317368479, +1.54511208122130332223736396761896289)
259(+1.86431671627461058540931287067777408, +1.64316630411472497661085484633557685), (+1.03820620191282388941706355611796757, +1.83198722862010334765167552229936607), (+1.51694657956171797622763165400552021, +1.78996187171613119892783495210300357)
260(+1.83760650160744667815481732873515121, +1.04775225176042159838540586726148897), (+1.59772835548584740003114451931533980, +1.86597030033664980574303923001048234), (+1.49031036577213572499015676210168173, +1.67952790788269539445205960923830468)
261call setMatDet(lup, det, info)
262if (info /= 0) error stop 'determinant computation failed.'
263lup ! lup factorization.
264(+1.86431671627461058540931287067777408, +1.64316630411472497661085484633557685), (+1.03820620191282388941706355611796757, +1.83198722862010334765167552229936607), (+1.51694657956171797622763165400552021, +1.78996187171613119892783495210300357)
265(+0.803909669702644352632749983627679001, -0.144187167078642429554263932537261998), (+0.873253041008853761451467995940269033, -0.246005089746218841221418192134019721E-1), (+0.104225856919366186103373185490407230, +0.324868654066301808723428126157934883)
266(+0.833514384509196351742889508500624219, -0.172637243473928240362376795400981916), (+0.459411206399800959793923598476685737, +0.606373381528406205324267734678109704), (+0.660086653723000396649548971191734729E-1, +0.237002330488062957608415857165072398)
267det * getMatDet(getMatInv(mat)) ! must be one.
268(+1.00000000000000000000000000000000077, +0.674075480535532548569592279904724561E-33)
269
270
271ndim = getUnifRand(1, 9)
272mat = getUnifRand((1._TKG, 1._TKG), (2._TKG, 2._TKG), ndim, ndim)
273lup = mat
274lup ! matrix.
275(+1.46883447653513210457675044907900561, +1.70165937239416351816672930335190229), (+1.38668073380891635939942611012266707, +1.78244117999954009566573792768867182), (+1.68132521230481089116860836230929427, +1.66276521255916574221088352405771357), (+1.12365573483243614491657629624197718, +1.68789964325788285996337605387996770), (+1.10566004502712268111059883555240234, +1.89415771760047030137552422549882673), (+1.47915377189368945053253492270298108, +1.27983155478393216955538662465267593), (+1.65039960029190820571425410682822451, +1.23559418676607870397523932157051297), (+1.34451335141470014472821895478598973, +1.10254543679316444985915825009278742), (+1.80728516260457047245026207779578558, +1.36818097902985257771372534556040915)
276(+1.76432439983987062436246728572743294, +1.53638356906335374826948408643884123), (+1.47596008761995619240400484475357407, +1.83876559768393145400858541022205470), (+1.32111530741831843516007690106393213, +1.80796177011061479388125154886968966), (+1.81148662644801676099439883995098352, +1.25498656192162497789441473741477379), (+1.21146672052045844421202254689199889, +1.85200932650707896916913196236406143), (+1.76937542109533987361851755254284482, +1.12184907792121049138045768926191542), (+1.36708796274578790564710235681783174, +1.43801607532905599433697841344349348), (+1.25078531009332581618474803188280772, +1.77833468248234235582081314643321592), (+1.33518407630279600738327371019519574, +1.15798470399650078924330773487770038)
277(+1.70517846350897351956476290536847429, +1.01366239394473681100523193577431484), (+1.60563658550557554588271505401159879, +1.01061562836256421724131740624809562), (+1.15707944042599747965526926184182027, +1.73209705467722939505510106160908977), (+1.27042225469145788597741656904329208, +1.04251453928953899178572608956256073), (+1.16147825615961332256848718464836779, +1.43970079983855934134956732011539539), (+1.20829736672795812546069716268584522, +1.61224127193583401012462111673553477), (+1.60376695951287173260436929147365136, +1.11574222124665748483659585616956100), (+1.26476708109055664568002902194512120, +1.90798023085095675161924537368711940), (+1.45689587877070694915258664472269864, +1.99082937272243973663172675413781210)
278(+1.99339239657120635367436644079959510, +1.28543170032937555586652757143629891), (+1.27229956903874843210234211432499131, +1.14154570421405095822180417920782369), (+1.02383444689460390337754230278096480, +1.10919810736331386979314486307677144), (+1.63780452725657757183934138820727060, +1.62416888553566363918110247295378219), (+1.82029993746312624275658588577517465, +1.18612590253606830186114678517613934), (+1.63717416610890783400493590427779526, +1.24711912717908451751789293814764134), (+1.05917800154428371883675960422225598, +1.00996050442944637293038825603913693), (+1.30452663726470788465207779242930690, +1.40706264833780863725267365265891367), (+1.31276234718302580499559113293786736, +1.38209338917050260423498414207828513)
279(+1.38628242138062350181489173220140174, +1.51167502842581989102258412738608658), (+1.40284752037080534778104284813166931, +1.29915430914324752802264445276599012), (+1.13439958310247461167238371660482704, +1.98077936179405028763575465580078693), (+1.39969252971308982028557700299076299, +1.31560699945751871610756979699988882), (+1.79709108810682221857920998424475200, +1.07259059658733185750488053311343170), (+1.14356426267749208740919300191964630, +1.64260669669027528163673096080368804), (+1.24704621378646169477363852588532832, +1.92802090075447506311569628890442281), (+1.94212990681948726490560121214726020, +1.74428380962162597151902601692565907), (+1.93601694714345722616138818253167114, +1.38653866027429229447662334836040305)
280(+1.62900640460247906024061559075749032, +1.45650214369956345766005669523652537), (+1.86555551658488305983275548153925831, +1.03540412523557207671496608324201234), (+1.89763307008417486147829552104468151, +1.22501365537642684267776974167496127), (+1.80728103176540893313821055695356298, +1.14377744655859707415960163792235815), (+1.95335078893459032566235592361964139, +1.24611499386678797296417572779228281), (+1.80031411178202198992684808107044362, +1.81636151017993712137077444922844354), (+1.84830999903005980816510069316570232, +1.66063317968162813023030392381947524), (+1.08562370801486619894484796234257241, +1.48609763798433353562014108474687049), (+1.48850772493893916529138458580254913, +1.87726075933564504437982754605283498)
281(+1.51733397755012355613799886867669559, +1.43593802494784056619158440111428619), (+1.27006443777866097359272803489516430, +1.32991725118024876921119529635975787), (+1.51962965424112978207143614818880466, +1.58713839041108636210739480814330956), (+1.71022773113878383960613681620838777, +1.02958896540167320949636840975620238), (+1.49005338092310393441710582355627049, +1.32033646159871626366242916482268803), (+1.93476616718458560885436336905246950, +1.11817020429908936315384810540429179), (+1.77450909769583054923491201466183561, +1.77764557186984882053817956418935795), (+1.33875839444419178724277423414766327, +1.75892295099090919773939142607634861), (+1.79766906756910387386435222093704893, +1.38081507085553652589719735632732008)
282(+1.33995995756887321240190189738532970, +1.38968981319914920117346161069021222), (+1.02850866871727297177572262163955230, +1.66047600937937923430373459476477472), (+1.64995896939652657117492474207181372, +1.25168931223369015884057754016749623), (+1.59542248484255704262037307231048580, +1.03121570270915257634481733370942197), (+1.63110847296656059485517736400177062, +1.20569841816587176405205100938958652), (+1.52884915329051818699368072568961494, +1.41182985306186275535389390304110650), (+1.70339403063105025494903003669289089, +1.44233301653404618400128050801315433), (+1.64813090719741711011878761646699160, +1.21121293253087477597913343623524401), (+1.06265651094817712742954661986651132, +1.01956217895147051460107108580573708)
283(+1.87170289938901488854669422597606959, +1.39381755291481924997830819171526803), (+1.05183494245265118461067445960605716, +1.39887451971410376384029221235219184), (+1.21555221375491418290063830570147189, +1.83926942261283761927247351737875079), (+1.77711972866686000751463105977859586, +1.31927434223286315374029378762507761), (+1.52913570432461650675357082900122091, +1.58394635410781286531532572901276648), (+1.06234697633616551516478866294375268, +1.25045640950442243429127780711286996), (+1.98353965497071222943363717050426732, +1.85014687728757126261291686672212254), (+1.84852039854096408510072167943893194, +1.61900765755425882936932386987166889), (+1.06242322118691768370533799733691534, +1.66986426522824326885888760677863443)
284call setMatDet(lup, det, info)
285if (info /= 0) error stop 'determinant computation failed.'
286lup ! lup factorization.
287(+1.99339239657120635367436644079959510, +1.28543170032937555586652757143629891), (+1.27229956903874843210234211432499131, +1.14154570421405095822180417920782369), (+1.02383444689460390337754230278096480, +1.10919810736331386979314486307677144), (+1.63780452725657757183934138820727060, +1.62416888553566363918110247295378219), (+1.82029993746312624275658588577517465, +1.18612590253606830186114678517613934), (+1.63717416610890783400493590427779526, +1.24711912717908451751789293814764134), (+1.05917800154428371883675960422225598, +1.00996050442944637293038825603913693), (+1.30452663726470788465207779242930690, +1.40706264833780863725267365265891367), (+1.31276234718302580499559113293786736, +1.38209338917050260423498414207828513)
288(+0.909976966788610223236766342820995702, +0.143869769200089384828249195676113719), (+0.872026130901175793610439779980501117, -0.186421617486675968514204291218723876), (+1.12554738130886470046461725034787001, +0.683701004964634251055656633673406885E-1), (+0.550585438563797099108902018386469338, -0.569809388784848915734223468356780824), (+0.467567433036509693287349850261480315, -0.951183890301432919712079045412710046E-1), (+0.489946070993837539084374091910068819, +0.445971960187017389500715640471600637), (+1.02978519856904089731301849269062591, +0.589208688660637763220866183588217090), (+0.100968294008208798889709065158587507, +0.180210909497032520308305517906783835E-1), (+0.492765683047962416192463046668269646, +0.430720793335891317955874015387398168)
289(+0.792296233017995193320305674582469760, +0.186238855864426443462689668489452977), (+0.133900831355628900660626253401589527, +0.623885021725243058502422014354451056), (+0.937297931876385177138250033793549514, -0.529168891622207958538674839977121654), (+0.171060091592835883540841717057144962, -0.827834085799074979897296665733003215), (+0.287843814346163937018015579413469912, -0.352047130731399160197032819139714367), (+0.676615307905155560337158058795838469, -0.246549408450264288858749587684705971), (+1.28201456004481305572721254884135776, -0.276478080021934659404258654788366909), (+0.874332454587369186939748592534807366, -0.211956697097314673401247053271902195), (+0.482697856808336068035326075258564911, -0.685055565712322558378701645123747679)
290(+0.835785122738864611292281230291328357, -0.304417221575923311748829127401325435E-1), (+0.534228930226639132959566672034094817, +0.223446885606012451979151328478458748), (-0.507992764405231207965669625226656871, +0.298062912489012816092173767017460988), (-0.729183054113069125115747232003080300, -0.555222397817562438037103419097095139), (-0.625760641168583687155221653389890098, +0.185473029449951387941740345543235355), (-0.898567724631560170684963444361690646E-1, -0.548912315054955018375787326205319264E-1), (+0.838138866476563150562480588916411145, -0.763569021145773367428900700390370274), (+0.462694418358920411386702901987825794, +0.371223254676150616514184219326744431), (+0.191646234003475835025031306561574096, +0.435709609392137210686496057378823776E-1)
291(+0.836580772099375655325085310735083854, +0.218876917978488029910188068766769629), (+0.629779233747996776223951377254602302, +0.209955379816641164898009141508925305), (-0.391567129113497444020932707287106445, +0.365064161219727542572824127361540327), (+0.635042765828833994512224035879880301, +0.263821107187712911970385883256465096), (+0.649950719477203643866858302770618193, -0.552011864066717794347499895592469690), (+0.494931652405561662532578512069304432E-1, -0.427763765307655245368813930527358968), (-0.275451619416676467189797020983887377, -0.484996320896025107722507781029737615E-1), (+1.16804527361200188920278032026932165, -0.510910577940337069718238944507546311), (+0.749102867194655520909309689349639101, -0.954435962827687971435150754261199755)
292(+0.976174518860202979629388205625305743, +0.141255629323664196462173440597494171), (+0.305717981643394445098046387335961845, +0.689991632262181489806405927702723444), (+0.245874800164955135679509722445909417, -0.926519554738046361980100972346023328E-1), (+0.435863165417539853229904250195026655, +0.419839176059843641621330798425003046), (-0.595762138460394435423754171821579728, +0.168807611494224364340193356015985949), (+0.335178684692581279961864955943312609, -0.879457041382457361384501326337575260), (-0.580246434236508912636755623936866968, -0.402879036945010210266957623350230730), (+0.526137798574648123328553185981191900, -0.579144877191117771863390933438284204), (+0.560194640047488512167723923570367331, -1.42965985554903375601569668648865051)
293(+0.909238845800588401808159576950716414, +0.267330675961148432582916158103868873), (+0.491927696218251051114358771653391001, +0.568887464367567729107927247460780462), (+0.564534306791826738119065886418937339, +0.568432119186597342249557686934453852E-2), (+0.408141802474316676225387135858698184, -0.594320622108629231492172442541447979), (-0.372145072831206711721096867753875741, -0.532442956092173008374484445957369762), (+0.914218986896119720502967117013968132, -0.428338858589127999398720317858076383E-1), (+0.643417840829905466598702972053208269, +0.295784038739352061534866381791771570), (-0.158630234639959343370347474741446753, +0.630095100428302039925493804324517279), (+0.941320595395678890813030176388759979, +1.12316776776221808132452713528069105)
294(+0.981647126950568573892933479182513540, +0.662073446363670427346279834310526539E-1), (-0.178771723275320253186750480897203424, +0.184302927782235633503080684107689375), (+0.180071846549437069041846541238405023, +0.621693796402918930817585279948118877), (+0.598525483656353230683357442192389603, +0.290494974225551143987918738783665464), (-0.748950885031405256730438250182371370E-1, +0.168879867923677780743978688373399019), (+0.220835524541710083868709373765304653, -0.764408424091464459800043944493866966), (+0.666555292432675343320615450123353423, -0.525944236967133407481479884101181156), (+0.326206921453453110515306090235799154, -0.934852272291078981133526892975368570), (-0.936242431037584680992498865028771742, +0.246297600435220083280853417460342179)
295(+0.865710517574535269644737436320011406, +0.162098683081557283706095903092038331), (+0.356090921360698980426723658056336910, +0.231431240996784555918189575624298061), (+0.265920333416595543548331391989195613, +0.338009575441053099026012935483996091), (+0.351868620027240013481509914443042607, +0.286337555330247000575556097917576592), (-0.298891856108769695416118286186052245E-1, +0.274291973725264318906651678239437255E-1), (+0.774631483038622357257714076909942921, +0.149063332514062343641213234682341570), (+0.496565221528694226878704004898916465, +0.326636464617223829722860172037602137), (+0.866520261215929311025046040278779099E-1, -0.183581146619164533208317439629558533), (-0.320918951404790177912069310979137957, -0.430618836703030496449802021672178550)
296det * getMatDet(getMatInv(mat)) ! must be one.
297(+1.00000000000000000000000000000000231, +0.231111593326468302366717353110191278E-32)
298
299
300ndim = getUnifRand(1, 9)
301mat = getUnifRand((1._TKG, 1._TKG), (2._TKG, 2._TKG), ndim, ndim)
302lup = mat
303lup ! matrix.
304(+1.80566578869581058141578270226458632, +1.03884695155080054125097825775429894), (+1.15657351243865293915204415641490956, +1.73040512680505975112531093552235292), (+1.52125295715086220963862659009134223, +1.98256408637839426199592045099042021), (+1.88189089917940502086837246312062031, +1.71221235676273972134626316857025709), (+1.58506591665148923868282312652522056, +1.81765146807590263634796193848866992), (+1.96339238616479188142783521384843985, +1.25457830899117333249647236752836071), (+1.80604846341477692682724488082687876, +1.18798431223466436290978566657102126)
305(+1.96279506656562760045369566908758798, +1.03009709703377503621996420724608149), (+1.17602818309053092459412865778213555, +1.76393184718810064968131888594498993), (+1.74891716195562106724367384334963313, +1.04873807323980366272404394247901157), (+1.08360991865485685263893722943676216, +1.46270837014184979848663471620905123), (+1.40278821589418830184901453248334768, +1.34892778236191229933853966292822617), (+1.95458108601306140886176602015074110, +1.35603775686494750103440913897717204), (+1.12574167726702576227250703146151029, +1.80686277383083757030674259742620705)
306(+1.20101049461011404425413577836751763, +1.16953360233378067226301981209698788), (+1.71378134983611226832312173971860827, +1.67556438852168953797023361925636860), (+1.46046483351688112770813608690734771, +1.66893603956722358683112155438282420), (+1.47643593425568032276118170222245936, +1.88370756616323893753988737672160472), (+1.14969978194765613479635469357786553, +1.38172188726345288088901577893234780), (+1.60392596636197288770149461945702941, +1.64788733728805426086037654273932975), (+1.68086591183865215998746813671475935, +1.86800098870868882740579197009539597)
307(+1.94431192160766706672463912886831269, +1.50993979300309517198620060462704552), (+1.54983154881776410638525065117080328, +1.57233407382165437717785444356349390), (+1.66601649774064439809006296307031906, +1.55315940491376316475936632849975328), (+1.78344692267055442906720721512528777, +1.40426427981056988203292650885480439), (+1.68024397344163921341662619044679474, +1.70487389688250024864603499474440258), (+1.86829894243484818093040802963092940, +1.72437089611661151413717253324086784), (+1.99264075822641513341209776818822137, +1.05290532140994638676322540044473005)
308(+1.60528031350022199822650737253202601, +1.53315626369714754726738329036714684), (+1.20340259980903530665744750951457162, +1.60399057371669060243829246439707914), (+1.39567858909953679363806449093487973, +1.06590979205172762581444208252767970), (+1.19730435417464887901420678088515136, +1.62826681073500588673724157773432395), (+1.23738162950198268919519893613789596, +1.45327171331049724831189038801168408), (+1.09300982137093189203624270978995728, +1.13319576842184693412028264481125882), (+1.78908432409693505756790983387022124, +1.50241498357085140926419522675794276)
309(+1.84222074653759030160488351424038110, +1.24866178347867250808759680223571179), (+1.19290773156709553752285868596922404, +1.01177538046590640260800947742045336), (+1.58790169132059790529390830469689802, +1.09029999610667932363037207680404613), (+1.75247114464861751031822081182277771, +1.64655127239592659318247855464408841), (+1.31654904736862769109719667699053337, +1.47909372260255592781157498053216553), (+1.66522957528547098580425071560551653, +1.42263048792725272031739477399496535), (+1.72532931878715909961908559314546496, +1.07040689210511759187000176938088381)
310(+1.62128217191025575310036694135233298, +1.03758872390851991395402243187264441), (+1.41709670158139194986675484561169371, +1.43807558449396617303641048616450542), (+1.24936655327067482724950197314922582, +1.34899612274884518263300926263302268), (+1.71245099913194160753503096906018901, +1.18647331340513145767949947364931078), (+1.16573510782465996663721788926681524, +1.39847639495772501739910553610905761), (+1.89662072348140873860639756008485509, +1.12688583539011761921397821346235402), (+1.98348212400908593582193903012540289, +1.05472127121861968911960494879104527)
311call setMatDet(lup, det, info)
312if (info /= 0) error stop 'determinant computation failed.'
313lup ! lup factorization.
314(+1.94431192160766706672463912886831269, +1.50993979300309517198620060462704552), (+1.54983154881776410638525065117080328, +1.57233407382165437717785444356349390), (+1.66601649774064439809006296307031906, +1.55315940491376316475936632849975328), (+1.78344692267055442906720721512528777, +1.40426427981056988203292650885480439), (+1.68024397344163921341662619044679474, +1.70487389688250024864603499474440258), (+1.86829894243484818093040802963092940, +1.72437089611661151413717253324086784), (+1.99264075822641513341209776818822137, +1.05290532140994638676322540044473005)
315(+0.676713473151130625196008630740288633, +0.759842076364987699841135362576604679E-1), (+0.784462018375653636144621785366668382, +0.493782014265006122549074614248295987), (+0.451064609719275137205125520242349455, +0.491301200820530090370635233938700407), (+0.376255281648238301340891446004026228, +0.797909206869682898857221317607814716), (+0.142199539113475069598546350615487970, +0.100338744261428368007797084848381194), (+0.470647856355162571536534679291220923, +0.339021104407234051521126184911687058), (+0.412423240160340251232479896132372314, +1.00407654264003023978211681979051591)
316(+0.886375207747854051772061851227379572, -0.158552800736414286108141891372445679), (-0.541133822252667056446920957243158435E-1, +0.819293016626621834782908811209711851), (+0.452871943546789205002359043261513065, -0.406760450288996253283286288810035846), (-0.457613655120721820127023633119992613E-1, +0.235688073115135111182184562317182941), (-0.266949185902572573336116560579345728, -0.689620480452649731494011799355843451E-2), (+0.328389358363968626714537651144807418, -0.243430747875694685338146592552806082), (+0.375237395843828834849042053557620144E-1, +0.905970871083018938293619162496513683)
317(+0.838143574822210072510699972358812696, -0.116596201362420116683741496779799446), (+0.435491400789268588174429306859860018E-1, +0.728861437775456714940595563456222540), (-0.231228652884414207381178235761376628, +0.951330052360220394051297052221354587), (+1.00218990497495920092685705259192789, +0.532226182325090059353070445847041996), (-0.233480877644696562396700731240666364E-1, +0.728981186483933765749116995916898130), (+0.267387678492412084063356671026670578, -0.699351498725430522285504590640942684), (+1.59759052913497223220466516080855293, +0.367295722416746601923210125056901050)
318(+0.778674394562435023646634597267782806, -0.710599614693782891493599498349626698E-1), (+0.276107642374067376918899029668355108, +0.239059737062740777895254770210689939), (-0.218030263199621760309542426189943159, -0.163825029877295151928215990827836515), (+0.184809961532172928553956304738186580, -0.144669585605621046239228044156926445), (-0.437271227272210186533341676523074192, -0.547029439565059790921241773637722643E-1), (+0.433625124970789985411143887524194083, -0.120542533288604511089928215451672578), (-0.542111617581094513824371984626313873E-2, +0.367540081455898680701600627942150247)
319(+0.897013164319254095439453360360452921, +0.919196092240721817577189558096877078E-1), (-0.922760570030132227858035210919317088E-2, +0.709820550717226595044449991779310708E-1), (+0.659047693175481802392041779405149693, -0.529605889836567662588020556145632252), (-0.237287831812250134068709687433319715, +0.131781977013427321951288336108654217), (-0.316056229206755178231829031976637752, +0.497870927323411096784012789423041156), (-0.435148731242978768364166682381463898, -0.736422761383201936390663732794544534), (+0.277747522909319892218735598242891915, -0.226948840648451448503677584318903644)
320(+0.902147685942674194871170368762137429, -0.583892457944504142694024338501969123E-1), (-0.452959578424529824767083563244443598, -0.117973050745961543888728095607972377), (+0.103580984415211232879714739522475236, +0.230272557377804112944054333577284013), (+0.515194627660076987612445041300173511, +0.600250446805873079191720198530297025), (-0.460787361405812924206638856138364656, +0.510425900208761184858241819861102917), (+0.213776295012744305615291603125065695, -0.494356307384510292659634255885957349), (-0.225393596952524273616577778729296251, -0.152377330270501464043782570641555465)
321det * getMatDet(getMatInv(mat)) ! must be one.
322(+1.00000000000000000000000000000000520, +0.529630734706489859590393934210855013E-33)
323
324
325ndim = getUnifRand(1, 9)
326mat = getUnifRand((1._TKG, 1._TKG), (2._TKG, 2._TKG), ndim, ndim)
327lup = mat
328lup ! matrix.
329(+1.23895217035497171331755247040739747, +1.70651416073090537912457595928248102)
330call setMatDet(lup, det, info)
331if (info /= 0) error stop 'determinant computation failed.'
332lup ! lup factorization.
333(+1.23895217035497171331755247040739747, +1.70651416073090537912457595928248102)
334det * getMatDet(getMatInv(mat)) ! must be one.
335(+1.00000000000000000000000000000000000, +0.481482486096808963263994485646231830E-34)
336
337
338ndim = getUnifRand(1, 9)
339mat = getUnifRand((1._TKG, 1._TKG), (2._TKG, 2._TKG), ndim, ndim)
340lup = mat
341lup ! matrix.
342(+1.92934930998234157898354463947338753, +1.36326516675869695672233111640206305), (+1.08324453881835616796965640658563947, +1.16293941545466969352063417407487948), (+1.01756044882365019909498995752847169, +1.42147630635850051831171578853017909), (+1.47900088369805907344242632042668830, +1.37349761533980489221795956211129948), (+1.51398207198393268662618676353616467, +1.48645728320548156225710111461579566), (+1.74931185181255573807729380769670899, +1.53714153592910541045884960763539463), (+1.39266563398217839145516275850359583, +1.46394590925368282728923894651114691)
343(+1.17782405764340492964855003392419364, +1.57959278608829774444009220157801097), (+1.00798820694036915573697979005987338, +1.54169837755157786651449657073153491), (+1.33162589426591819917402493936803438, +1.91577166715892370977676772346000240), (+1.83000978798487419772636711873208159, +1.81183317687388185657943975594331721), (+1.19358404895203804109507102201072477, +1.48522088917644932113013887814179791), (+1.37715018322048181281340219881883492, +1.34455857044801450958779361876951538), (+1.06034289284177777211972969871696880, +1.80327258194337284300468849888978690)
344(+1.82092989991803034557149824172082062, +1.25813783342369942956574809411789774), (+1.26486709937043291934509042391188969, +1.97178076386621390926056453489555694), (+1.37157005607648607097151480010767944, +1.88272039672014858329608847717427617), (+1.79445749856272362905606487972633023, +1.05982494307413525590258472090723971), (+1.88664221598738560498551621968498663, +1.30319580488912616046166251360301553), (+1.22436662187543663811067417632545333, +1.35970634358642523209942678507143590), (+1.58720379303655404722928792040607457, +1.44418142544911710524038010194972847)
345(+1.14115444941959225786263234250743904, +1.95876754776668524442132224124640938), (+1.62444574847795378006481567104138889, +1.16977764551585906035349869448342191), (+1.37321065212601675840377499592345127, +1.94504401114376863478580387574399285), (+1.58756237717166056285074317240558680, +1.72389656166561595685384562500951545), (+1.82256580487317738862987349509860482, +1.19432596032373731176673274049339903), (+1.10770316960360758320492252829237627, +1.32270739531361879989364597789241955), (+1.54845666137128546852530012712619597, +1.44296740553442816434236852262312293)
346(+1.90171293755310610366954653783314711, +1.53684081110383070255238072138921194), (+1.27040189729706022077439656192072393, +1.68826891842437233718040992279843602), (+1.05414649023897276776294251472567495, +1.03081327686856007950725852020995097), (+1.05557699203242324067413042025972769, +1.60756708749255407629649499020197657), (+1.05537004059059899494163310658291870, +1.11646524853998012887927053488923299), (+1.36842092468632132686029467062326678, +1.63508241838850550067288597216724480), (+1.88422508513142107297836480035639990, +1.45143347623459223279727117411096548)
347(+1.21672737370378845024195629165524356, +1.05960813569554324452220122033883912), (+1.21237168694807629938645208946716392, +1.79556306670892261050936444850680067), (+1.44445622609165821747095988057131317, +1.57255556318711863821479317031019933), (+1.16262292003804036916356690335113185, +1.22429658917521676709175333949379300), (+1.69810694415079955795335276452536221, +1.35678990827419943962864777506835561), (+1.76389035353295511682011632923176408, +1.79903089890691693028787619159259938), (+1.31319819287866878996926390095313364, +1.68122173687159041975641384286116274)
348(+1.67572907779898223415384388887694821, +1.07701866635232842913514588902507455), (+1.94588957283429827010727839778591625, +1.15125217579008702098834197731247380), (+1.36478562077034029026530263856320157, +1.08789492064917483262028114110615382), (+1.40532853616405521631368586218420416, +1.98533217020772557461041652195357869), (+1.97635724314087087974718858281501750, +1.30563220732497044303669193442988594), (+1.93651172162451453251452970085222751, +1.38774590052343590995239808539644809), (+1.20409946420354760916227131502836523, +1.71152672402541955987695303079175869)
349call setMatDet(lup, det, info)
350if (info /= 0) error stop 'determinant computation failed.'
351lup ! lup factorization.
352(+1.90171293755310610366954653783314711, +1.53684081110383070255238072138921194), (+1.27040189729706022077439656192072393, +1.68826891842437233718040992279843602), (+1.05414649023897276776294251472567495, +1.03081327686856007950725852020995097), (+1.05557699203242324067413042025972769, +1.60756708749255407629649499020197657), (+1.05537004059059899494163310658291870, +1.11646524853998012887927053488923299), (+1.36842092468632132686029467062326678, +1.63508241838850550067288597216724480), (+1.88422508513142107297836480035639990, +1.45143347623459223279727117411096548)
353(+0.866531047410470388898350622991170296, +0.329727614318349456856715358966405207), (+1.08027194458077177948897314602925791, -0.712046375498218760616443660159515568), (+0.799647592404719039388817748786408865, +0.704231095285678317929713943423635856), (+1.20293140125902830323775968728175892, -0.171631137540965070363276443110586567E-1), (+1.27618432106500572916793596529006440, -0.121110486598000652067685406705438902), (+0.461055777465884391702298710991910042, -0.545348452175148328194688266722313450), (+0.394294822255893290099685537676809907, -0.436025807032925503481478883560411023)
354(+0.659427378262076667911552244870128032, +0.242798100130668925367064429901851198E-1), (-0.887613600836275479093217657130426562E-2, +0.597170912460881882566117300112786580), (+1.20199522679453705676084481926353243, +0.395939158614600895142477008898730116), (+0.506006045521010819314467957948840466, -0.579914354180821979894061017591115928), (+0.968278535123745270378829158487647372, -0.168237176200003278781057813706804508), (+0.579641782275932130987566995094957194, +0.407418100079159396557310753606325654), (-0.150953000588209290584047083451806579, +0.439026514676021992708227725941749977)
355(+0.780722656261981332404379138305074828, +0.199686471229516243820643731220379605), (+0.240764495225381035773061222842048701, +0.130876302067610536689126185281439651), (+0.615729311291568591583448789772927574, +0.318215477677402785305034275662191190), (+0.538936420221740500712872063073677763, +0.388733050963776766161868034435458578), (-0.380267762206886967666046786458833738, +0.604322199959384279830001999341072600E-1), (+0.225662061955559192831618783934406610, -0.569593447240693633569960780136210254), (-0.402291114344943721079181091325014901E-1, +0.124941097480986470462452385524309090)
356(+0.809910439070373424138861694293980087, -0.881756370474086105134284056034141147E-1), (+0.539952117806113563942103939882340726, +0.259559642582210295584501393632426739), (+0.686635201007005287501711331041084345E-1, -0.223807349381252393385159466078774464), (+0.368145868406369718543407665920314906, +0.896978087386021713361756204112826172), (+0.468010860634985999472337160722124618, +0.775704777872700681516910899229247505), (-0.431432403882227355009015126641534276, +0.467960426320269344728468037081271510), (-0.737022988497620123829872362039655953, +0.761386967759763416224472238329942715)
357(+0.964170684175594733727650087933817592, -0.623183905478110325346671553824423221E-1), (+0.474928188210935566043231645789108864E-2, -0.353882070226763522664539037108336463), (-0.461303730977431626199859420932874629E-1, +0.658230776263042264708525887240896034), (-0.356591759535533645008594561824697619E-1, -0.572808464567758554626174477720038186E-1), (+0.465771450514293555468401363771438402, -0.209573535724585947321077254855328045), (+0.957288880648306525788265436106057471, -0.466838073162579083264557269330731518), (+0.950728283543329400290915498012337607E-1, -0.637846656664403530448028058432059591E-1)
358(+0.902658059235546469250734152396979156, -0.678882222084384319724330019733705120E-1), (-0.224913320208913508213358601738255140, +0.346160827568647729502823997622889437), (+0.804495119489034106916378419115437030, +0.488243017642674253771981386703648172), (-0.835511507191622717648600565292944010E-1, -0.905304744679322302735349934884733369), (-0.647039790720048855724689136404107682, -0.591127635751043845421007382153800328), (-0.125271248911727131583698735160411081, -0.804734160022329264052990721725081171), (-0.918794161811214180190554141354617789, -0.152590039297130129481007819534889782)
359det * getMatDet(getMatInv(mat)) ! must be one.
360(+1.00000000000000000000000000000000116, +0.962964972193617926527988971292463659E-34)
361
362
363ndim = getUnifRand(1, 9)
364mat = getUnifRand((1._TKG, 1._TKG), (2._TKG, 2._TKG), ndim, ndim)
365lup = mat
366lup ! matrix.
367(+1.90938610344036060377676743557887033, +1.13137662442532403380472535359983212), (+1.86178687591959388218191718160183300, +1.39607995741807628312540330779467862), (+1.09277392407255302699364494462850695, +1.80060937552341719417606782368573466)
368(+1.34981952887116114272470144294175217, +1.90813831025098440861747856424291723), (+1.43478391556092880541753372858244385, +1.95282753314376188052113354054995130), (+1.41953851893597906639061592322890020, +1.89156396144040979437345621078384191)
369(+1.40472376511568749679062826323079236, +1.75132125098320519971825815817737648), (+1.56736688347451892656294620161530161, +1.96382330922212165358585296319063072), (+1.88296093666534926922028336875547245, +1.39991146728814083622642822161773415)
370call setMatDet(lup, det, info)
371if (info /= 0) error stop 'determinant computation failed.'
372lup ! lup factorization.
373(+1.90938610344036060377676743557887033, +1.13137662442532403380472535359983212), (+1.86178687591959388218191718160183300, +1.39607995741807628312540330779467862), (+1.09277392407255302699364494462850695, +1.80060937552341719417606782368573466)
374(+0.946768833260505908401055124123958245, +0.356223983757026728883371207259148771), (+0.302002259255210116731645980871678801, -0.211548210475941486101413745062663115E-1), (+1.48977688849287737698469214735550692, -0.694121650913016495326039287437178765)
375(+0.961504761583031765635771377211176886, +0.429622011621559048172938020586769545), (+0.849201149140632387881067754936746208, -0.567578263774561598674209699149767135), (+0.271260725465689647276798569827460756, +1.12580362499694226033084719426042091)
376det * getMatDet(getMatInv(mat)) ! must be one.
377(+0.999999999999999999999999999999999711, +0.529630734706489859590393934210855013E-33)
378
379
380ndim = getUnifRand(1, 9)
381mat = getUnifRand((1._TKG, 1._TKG), (2._TKG, 2._TKG), ndim, ndim)
382lup = mat
383lup ! matrix.
384(+1.89636919885334843190535393834516817, +1.72188264262066646939183361797655735), (+1.66579455095478036594951281609024692, +1.33906447315671263634970897562742815), (+1.20107997471888487212571623860678956, +1.79258862846355231788132829932476449), (+1.35185765920584613101134515977352294, +1.85617051141696379367572516655924470), (+1.37136447478593868863293722885284234, +1.20255855286921330800381520121215510), (+1.01219007004011811696604517949990408, +1.72313377029384361261312584283011769)
385(+1.17090802162271759968428048619499583, +1.89856696159298958085208566420746101), (+1.72468329701074332869965576650703835, +1.59504919546204325855373308894072305), (+1.21327357340589785713216830987058118, +1.57340863677898972102880435792515380), (+1.99798973655461680722625325225480393, +1.05067464977793745275729456294496525), (+1.93496153761767597437615831531054126, +1.86987578240282823074773555173211835), (+1.42694412572649524638344646818208191, +1.91155784412762938251099632366300103)
386(+1.81023633038097719565286214604780721, +1.90225772172114810572564601567209980), (+1.33983162786275688032695391718923270, +1.26506356466982484132802809147825110), (+1.87657645485795706061033560577250889, +1.21480137963559319588348228757985612), (+1.26530002468156162424522489120242319, +1.71170999647224810337135052203217397), (+1.14178751123862340088432478631178346, +1.94095360734894705831257823317923748), (+1.17598746338897432316909001873075594, +1.06243129797140873914552632497274928)
387(+1.04091700394777165450260017810285584, +1.45633622825733838073954615859436244), (+1.41601777723776200420110770079668780, +1.59245672010766763740596145610332976), (+1.87495253055114358755376843413827132, +1.79975341391881536507903180750808312), (+1.40636096976675621012725292724484940, +1.45498005492024940002921948296308744), (+1.32018294453949403747936752049384935, +1.91355092541638407651648769394244002), (+1.81629383631197786889357518429095701, +1.30705204009813594314253072742081970)
388(+1.49308377019865758035002238548149415, +1.08855592577732162257361031681780267), (+1.92334746745224789254355292398924375, +1.62718878962416429788937239630214232), (+1.27169270751960642631985815243879862, +1.73839812542062384562268875713477401), (+1.47005930221121948562567518288116037, +1.63213170208011303349611075599365013), (+1.16215287215969554437530627874896444, +1.27805408956728930011222830258776414), (+1.70316827296833200610651292568859698, +1.24522189723360903061276153610405978)
389(+1.08832604349533347545081277308517356, +1.52955031370543680095732299627554506), (+1.33840738931746822154890912994708047, +1.21428130454597817917818587848402445), (+1.50047456431930028031035396624941983, +1.75436946014560378153414936216186956), (+1.54287352879864680519646795195682372, +1.19819925431872938246911179665994166), (+1.43975831135284324356280615083216660, +1.10389409991921589181240675380286625), (+1.95575789388388930989561930931133824, +1.55893715560079423630593187695521934)
390call setMatDet(lup, det, info)
391if (info /= 0) error stop 'determinant computation failed.'
392lup ! lup factorization.
393(+1.81023633038097719565286214604780721, +1.90225772172114810572564601567209980), (+1.33983162786275688032695391718923270, +1.26506356466982484132802809147825110), (+1.87657645485795706061033560577250889, +1.21480137963559319588348228757985612), (+1.26530002468156162424522489120242319, +1.71170999647224810337135052203217397), (+1.14178751123862340088432478631178346, +1.94095360734894705831257823317923748), (+1.17598746338897432316909001873075594, +1.06243129797140873914552632497274928)
394(+0.692266072325896960186200167299655540, -0.126123065788949585933653619065716792), (+0.836273793659628979699398102838652565, +0.920411877014636213006789515133610847), (-0.180611978628454030581831664326344112, +1.13411192135825477682085710649783023), (+0.378248911314382326223341005820493646, +0.606756464116967028985684726545991351), (+0.126933096810802516515100019406782812, +0.784035006380024201369469916738848499E-1), (+0.755074958093265118595524940625583400, +0.658055899682824755190739313743120335)
395(+0.675020753298340200459858397411448608, +0.951659101534486045761124609273729587E-1), (+0.705392511498477715798455807832642216, -0.457330869797952770719362917217832379E-1), (+0.799368055386817744653301579757084708, -0.710283462522655100784974530896305886E-2), (+0.420590838018206427314346246168204214, -0.231576120522998386871796344709182051), (+0.641042005139693650299792297529043495, +0.445207511689452666646003850095142697), (+0.560865985595928866206622661745587764, +0.483191526823717245611111237835870846E-1)
396(+0.831143593317422864013839020745442859, +0.175401210337567606338138632400552305), (+0.634094392491778047482607458777923291, -0.328883175236768362880338708771267784), (-0.484075751987232664301747365702704042, -0.684781041303178560677657219643609184), (+1.16935898954480647600849428548680173, -0.678366224362215368760844009592045206), (+1.22558784475377257880332170581415650, +0.702912227934399901956692516832349421), (+0.179082688717648347470459400333225294, +1.06077767070103741433422950815608253)
397(+0.707664117672794250092899333917877554, +0.101307646078887857820545822378669709), (+0.389428491985982973729288027950540965, -0.209415368156317464398657738441438278), (+0.158102647054882203746691857441922338, +0.283012470294352950787506813973771941), (+0.406543018382997098544824150603698202, -0.898037128509589275975632209331631306E-1), (+0.225809968540521657484195626379423606, -0.816784830417808802019190662738917877), (+0.556265392009975058796159809434725087, +0.827379072241163688740652434836867254E-2)
398(+0.972852156762300859657579237510510229, -0.711138554193089448200749454756909922E-1), (+0.268473480755586100700511797492250851, -0.519936558866753663688031150412329601E-1), (-0.907197132877118994838685784832082331, +0.530298586137001174586340193578238835), (+0.189540776746308378153636713973519574, -0.142699738297499063078568787096440788), (+0.762924731896336627686088285726628769, +0.486241637105294809963830570434388010), (-0.515599929020135009429779050333593998, -0.701259449710603665548021691700659015E-1)
399det * getMatDet(getMatInv(mat)) ! must be one.
400(+1.00000000000000000000000000000000135, +0.174537401210093249183198001046759038E-33)
401
402
Test:
test_pm_matrixDet


Final Remarks ⛓

If you believe this algorithm or its documentation can be improved, we appreciate your contribution and help to edit this page's documentation and source file on GitHub.
For details on the naming abbreviations, see this page.
For details on the naming conventions, see this page.
This software is distributed under the MIT license with additional terms outlined below.

  1. If you use any parts or concepts from this library to any extent, please acknowledge the usage by citing the relevant publications of the ParaMonte library.
  2. If you regenerate any parts/ideas from this library in a programming environment other than those currently supported by this ParaMonte library (i.e., other than C, C++, Fortran, MATLAB, Python, R), please also ask the end users to cite this original ParaMonte library.

This software is available to the public under a highly permissive license.
Help us justify its continued development and maintenance by acknowledging its benefit to society, distributing it, and contributing to it.

Copyright
Computational Data Science Lab
Author:
Amir Shahmoradi, Apr 21, 2017, 1:43 PM, Institute for Computational Engineering and Sciences (ICES), The University of Texas Austin

Definition at line 395 of file pm_matrixDet.F90.

The documentation for this interface was generated from the following file: