-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathos_integral_operators.f90
216 lines (216 loc) · 17.6 KB
/
os_integral_operators.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
!
! Representation of some "problematic" operators in terms of functions which
! can be handled analytically by the OS integral package in import_gamess.
!
module os_integral_operators
use accuracy
public
!
character(len=clen), save :: rcsid_os_integral_operators = &
"$Id: os_integral_operators.f90,v 1.3 2021/09/29 13:46:03 ps Exp $"
!
! 1/(1+r**2) - needed to implement soft-coulomb integrals and complex
! continuation of the Coulomb interaction
!
real(rk), target :: operator_1_1pr2(2,200) = reshape((/ &
4.5e8_rk, -7.08334995255803351635367391689e-11_rk, &
3.65780144165179293762485725061e8_rk, 7.42758210498556093113898491056e-10_rk, &
2.97322475256665217208455243798e8_rk, -3.9259005803904914704546508383e-9_rk, &
2.41677017473179903539800298518e8_rk, 1.40406104060205928286325882437e-8_rk, &
1.96445898428334000710084537223e8_rk, -3.84663435523469747664800975701e-8_rk, &
1.59680020106164842721688975513e8_rk, 8.66042339444973915457875420693e-8_rk, &
1.29795068388292676461316833238e8_rk, -1.67755964855721442213393257612e-7_rk, &
1.05503241837775557423105782879e8_rk, 2.88801227727631119941857179079e-7_rk, &
8.5757757798478486689600900679e7_rk, -4.52686270993339193309089493962e-7_rk, &
6.97077444684666438000695106546e7_rk, 6.58165237508309818989213007394e-7_rk, &
5.66615751580115742876961642613e7_rk, -9.00601841732684395235292794772e-7_rk, &
4.60570647331635896993123505241e7_rk, 1.17331080678282215687789253685e-6_rk, &
3.7437243950936828535929862384e7_rk, -1.46892807129755848565454931567e-6_rk, &
3.04306677544903546511230452448e7_rk, 1.7804888030629723043270273003e-6_rk, &
2.47354089739558995095698252773e7_rk, -2.10210131732875930809370537152e-6_rk, &
2.01060476899516856833428160186e7_rk, 2.42924909120286704785283485869e-6_rk, &
1.6343095605827772666837668091e7_rk, -2.75881901734130413986030924204e-6_rk, &
1.3284399704011090922974581179e7_rk, 3.0889635056512078690173680452e-6_rk, &
1.07981547530689821701552480336e7_rk, -3.41888464880745479049782775714e-6_rk, &
8.777223560656640702628404532e6_rk, 3.74860082961305206973186152996e-6_rk, &
7.1345202208599874660186515242e6_rk, -4.07873086773607023422168842434e-6_rk, &
5.79925741096789562349849748088e6_rk, 4.41031209881521410955341994626e-6_rk, &
4.71389602629960285340175250876e6_rk, -4.74465680537449359973086706872e-6_rk, &
3.83166570684341007856457812775e6_rk, 5.08324476293782962659510770673e-6_rk, &
3.11454941031545813195273965339e6_rk, -5.42764662762415283909236336635e-6_rk, &
2.53164518292169429536870635151e6_rk, 5.77947204196642701808094800219e-6_rk, &
2.05783453329817565212477420737e6_rk, -6.14033665248686865781394635813e-6_rk, &
1.67270002723975808858164640916e6_rk, 6.5118430652185131181795007511e-6_rk, &
1.35964546024190676742780723321e6_rk, -6.89557173916967316727818912571e-6_rk, &
1.1051806943573693535876574728e6_rk, 7.29307874040403556219851641388e-6_rk, &
898340.341579136702645253620554_rk, -7.70589806893183021196076688388e-6_rk, &
730211.243671584519265906757858_rk, 8.13554690594860882175967154754e-6_rk, &
593548.386624949003760784117423_rk, -8.58353261929578167026435915484e-6_rk, &
482462.698730407608677440134579_rk, 9.05136073189696805673758198809e-6_rk, &
392167.278879844355830580775638_rk, -9.54054332572864641679250248376e-6_rk, &
318771.119567856758523167889654_rk, 0.0000100526075452686664802087807368_rk, &
259111.435713836114013495669209_rk, -0.0000105891039984802315410355156884_rk, &
210617.374022785690073893409754_rk, 0.000011151614945531625333849165767_rk, &
171199.229852768973536802554176_rk, -0.0000117417622274001697791887293095_rk, &
139158.397725585559597301488025_rk, 0.0000123612149270888945083726501284_rk, &
113114.175070800098513895206299_rk, -0.0000130116967818779461617171725812_rk, &
91944.265032272421738550452381_rk, 0.0000136949933805228704332196490672_rk, &
74736.4144859245774095495227059_rk, -0.0000144129591879354126531911408967_rk, &
60749.1032556668579136213349468_rk, 0.0000151675244439874833570616267941_rk, &
49379.5905483626357247718340771_rk, -0.0000159607019842978675876837419056_rk, &
40137.9416657724661170796471444_rk, 0.0000167945940303148945373969740145_rk, &
32625.9157533329097443836888833_rk, -0.0000176713989944554399857181234171_rk, &
26519.8048172780131749545697004_rk, 0.0000185934183440065774717617555729_rk, &
21556.4845095474843946793971193_rk, -0.0000195630635652703191123617321925_rk, &
17522.0755813260735644724436787_rk, 0.0000205828632672428384993435164973_rk, &
14242.7274049124865980970246019_rk, -0.0000216554704620940507669614276781_rk, &
11577.1264077205317983223230671_rk, 0.0000227836700579233124738783364034_rk, &
9410.40659207671199282937577118_rk, -0.0000239703865977481306876001927847_rk, &
7649.19973311727397040161648717_rk, 0.0000252186922774478498448085171931_rk, &
6217.61195806197211454546804815_rk, -0.000026531815274431707407173521991_rk, &
5053.95332974011305134286302946_rk, 0.0000279131484181204098290419854821_rk, &
4108.07950568094772491445966422_rk, -0.0000293662582329080141800349768041_rk, &
3339.23091962221232965378633984_rk, 0.0000308948943840900508266835340997_rk, &
2714.27637151163790344975601363_rk, -0.0000325029995572878712363927517282_rk, &
2206.28533883481480916045450831_rk, 0.0000341947198021522725917482578094_rk, &
1793.36748735226668270067733198_rk, -0.0000359744153715769699433384822795_rk, &
1457.72937347746102431579338862_rk, 0.0000378466720882814944328598902506_rk, &
1184.90768974311596141434713083_rk, -0.0000398163132714223434459422726463_rk, &
963.146012392592210534099516107_rk, 0.0000418884122568511479424877580046_rk, &
782.888193922399869821342096186_rk, -0.0000440683055457513166159618830493_rk, &
636.366569862560595946492768358_rk, 0.000046361606617643590866529905076_rk, &
517.26723481384012350663829542_rk, -0.000048774220445150988656654659265_rk, &
420.457963826955755782799825461_rk, 0.0000513123587494505353287371357138_rk, &
341.767054720048037674690387487_rk, -0.0000539825560370094880357400532256_rk, &
277.803561214239759194462274221_rk, 0.0000567916864600043021121751288804_rk, &
225.811170378988511143882954267_rk, -0.0000597469815447481458810676954249_rk, &
183.549427678520626283192822289_rk, 0.0000628560488345038776934480156198_rk, &
149.197191372634278273962093556_rk, -0.0000661268914952310974651256737658_rk, &
121.274155931826690374809076233_rk, 0.0000695679289351095995139768667246_rk, &
98.5770627561200009262261766856_rk, -0.0000731880184911157203359748237986_rk, &
80.1278493918077787273584014016_rk, 0.0000769964782385724324300549853083_rk, &
65.1315028937360517699087178388_rk, -0.0000810031109826481979091048978395_rk, &
52.9418011514790341423672934358_rk, 0.0000852182294947057342911771457607_rk, &
43.0334659056716805473139677662_rk, -0.0000896526830619680926825415321169_rk, &
34.9795274731198116305985739525_rk, 0.0000943178854312077544530102868018_rk, &
28.4329257821058099369969300969_rk, -0.0000992258415542281362278903972626_rk, &
23.1115548702589030717350617457_rk, 0.00010438963138073963491144241953_rk, &
18.7861063829438959868524360534_rk, -0.000109794206156931781509751332387_rk, &
15.2701882245680293397188980236_rk, 0.000116274671352971588560303529073_rk, &
12.412292556026482258295502879_rk, -0.000111084307955770210283757334721_rk, &
10.0892670234749965925124163457_rk, 0.000214682524824021957071312982616_rk, &
8.20100788081705252339430532175_rk, 0.000331741772872387959122900643559_rk, &
6.66614632210007333234816179503_rk, 0.00190036474573576389135469483433_rk, &
5.41854213938654343893056062127_rk, 0.00482917790615770956014459535296_rk, &
4.40443361091090865457869644651_rk, 0.0113121721988327691455020017934_rk, &
3.58012080258878527137035684085_rk, 0.0205124431251356231732254389116_rk, &
2.91008245177707394467317184311_rk, 0.0330202363569233665925065591885_rk, &
2.36544528609683672790297961273_rk, 0.045846930433632881672735720487_rk, &
1.92273981725187715126248263352_rk, 0.0584433203644373884509119060855_rk, &
1.56288899454782691947061469759_rk, 0.0676552741846815909905823897807_rk, &
1.27038613719972499075468733621_rk, 0.0741108480941836331825414263142_rk, &
1.03262672091191255294878015825_rk, 0.0759670671026262033871144459831_rk, &
0.839365224097723860918919101542_rk, 0.0753697755931279109853221814493_rk, &
0.682273628172674352987579540064_rk, 0.0712142669712150007196044753306_rk, &
0.554582546829112795011632906642_rk, 0.0662588613951307291095428443283_rk, &
0.450789519845878166599624633087_rk, 0.0592404226664553014014696130899_rk, &
0.366421901238616151353197526775_rk, 0.0529220999896855373105895071378_rk, &
0.297844124134089122374222878727_rk, 0.0455174531059927491148924597807_rk, &
0.242101037032263725703670750771_rk, 0.0396985947361165701966000180551_rk, &
0.196790560507112970111155363792_rk, 0.0331582660742622816204039634745_rk, &
0.159960176872529380078254977915_rk, 0.0285991158259111625764314427148_rk, &
0.130022792349158611855958983682_rk, 0.0232866716340156259571752828493_rk, &
0.105688346067187574081718080132_rk, 0.0200940907290853986162389908394_rk, &
0.0859082188023005130776105250773_rk, 0.0159253752891246757398138226431_rk, &
0.0698300459077316624776874601434_rk, 0.0139259153908061619733122721565_rk, &
0.0567609872426470575130866408555_rk, 0.0106582088865549179973836123977_rk, &
0.0461378713257185409429222584401_rk, 0.00960754979945499942471789004031_rk, &
0.0375029271666573773574636201189_rk, 0.00698208424076493281201974825102_rk, &
0.0304840580125247878041995284879_rk, 0.0066563247045933752157644503337_rk, &
0.0247788069656911241335845244577_rk, 0.00445211878300690629875953626256_rk, &
0.0201413235203369961777205429521_rk, 0.00467609831927019943598849857681_rk, &
0.016371769379880851077655728025_rk, 0.00272117833947684896503728561488_rk, &
0.0133077070311375238085299449573_rk, 0.0033693810624505901194758613974_rk, &
0.0108170999919054524255308300173_rk, 0.00153534096146176713500258174581_rk, &
0.00879262309886296924403766436679_rk, 0.00252351434923148253530311859494_rk, &
0.0071470376548715164447590488352_rk, 0.000715592198935830080753631658933_rk, &
0.00580943214167304048663978268608_rk, 0.00199094300273554759726421071406_rk, &
0.00472216650288664677431903895672_rk, 0.000138458178045183850654685266806_rk, &
0.00383838832043968489058497087703_rk, 0.00167123412879440457518180502137_rk, &
0.00312001385158304078788710176129_rk, -0.000280143793152114499253240181288_rk, &
0.00253608692539866882641083815543_rk, 0.00149694294522742421379065501118_rk, &
0.00206144498041722530664554071482_rk, -0.000597093643396545689217603681699_rk, &
0.00167563476027910639720175229192_rk, 0.00142326407338239037504331802462_rk, &
0.00136203094262906040992567336196_rk, -0.000850802041606411742383241815176_rk, &
0.00110711972122731725447429498277_rk, 0.00142074716425879564487701233922_rk, &
0.000899916469419202701640312450768_rk, -0.00106728479480576357532056847557_rk, &
0.000731492390935055838744613353348_rk, 0.0014702891688875443232755952493_rk, &
0.000594589760471014257491844644555_rk, -0.00126431309572921668290377427232_rk, &
0.000483309174009393376987617110184_rk, 0.00155974279825023348475496008078_rk, &
0.000392855331878910299597730649162_rk, -0.00145421878480300736833930501654_rk, &
0.000319330399846060369378947569763_rk, 0.00168164247085117948929532252608_rk, &
0.000259566043760036243705753477214_rk, -0.0016457787975726850571765273815_rk, &
0.00021098689979318065435758855998_rk, 0.00183169068374254780341961270413_rk, &
0.000171499596940697460155348659861_rk, -0.00184547735516986428956911032523_rk, &
0.000139402549540529908309495686867_rk, 0.0020077567328108381651340329546_rk, &
0.00011331263259539684084344175127_rk, -0.00205835116823717427403920031924_rk, &
0.0000921055801921783351274108339623_rk, 0.00220921848365839557272858337112_rk, &
0.0000748675386691388306318057884781_rk, -0.00228855671008257221389413026902_rk, &
0.000060855686861532748861327721684_rk, 0.00243653285614178177771825236811_rk, &
0.0000494662264744054584809583354282_rk, -0.00253975277665188740069045411617_rk, &
0.0000402083632247009717148631755056_rk, 0.00269095817446302860811491466296_rk, &
0.0000326831575488378107368981975731_rk, -0.00281536021518493318621271113172_rk, &
0.0000265663335110825829462943387816_rk, 0.0029743766440561462672023135186_rk, &
0.0000215943051147200466515705917578_rk, -0.00311873931344739509540098700301_rk, &
0.0000175528178622425938934206654905_rk, 0.00328918166162109040136541711198_rk, &
0.0000142677161070138449228793506312_rk, -0.00345331052845501490885528219781_rk, &
0.000011597438343403055280099259722_rk, 0.00363820497482834394425700382958_rk, &
9.42691704265943930091545916608e-6_rk, -0.0038226335133348871347370475062_rk, &
7.66262016643812299980476152323e-6_rk, 0.00402466445673984787647627532882_rk, &
6.22852068702832618982383915706e-6_rk, -0.00423045085135044523447564981576_rk, &
5.06282043297116163124690338804e-6_rk, 0.00445211522896423298968869435841_rk, &
4.11528708412134878419827977145e-6_rk, -0.00468069486011599080523118738626_rk, &
3.34508956202446063542958949756e-6_rk, 0.00492438511124348745117844605226_rk, &
2.71903853831720816483828854409e-6_rk, -0.0051774465666742947803028590104_rk, &
2.21015624119188174596609521362e-6_rk, 0.00544546916373225033153349631966_rk, &
1.79651393006830507487838847679e-6_rk, -0.00572482335848825195991563156047_rk, &
1.46028694296697213850688448901e-6_rk, 0.00601934571550096936437141699489_rk, &
1.18698659671330671554866640461e-6_rk, -0.00632675295469923516791149226388_rk, &
9.64835841039841921955997109525e-7_rk, 0.00664965471930407557684001880852_rk, &
7.8426176229174531282221425268e-7_rk, -0.00698656196544283232585734812277_rk, &
6.37483067720693721234905610468e-7_rk, 0.00733914376534924989651354129407_rk, &
5.18174774252746897877437069936e-7_rk, -0.00770626134438722471302577302632_rk, &
4.21195652508731102128369844768e-7_rk, 0.00808873077304051653077318375951_rk, &
3.4236668110308983158003417906e-7_rk, -0.00848534050330168224367395996074_rk, &
2.78290964380538142294878211366e-7_rk, 0.00889594684202694117664987315706_rk, &
2.26207353491066705428825341496e-7_rk, -0.00931877837667303716650197730881_rk, &
1.83871463047080118041992607384e-7_rk, 0.00975240140050917503357338363912_rk, &
1.49458956136052006146952754239e-7_rk, -0.0101938409996260064063148826593_rk, &
1.2148693004938291382433925851e-7_rk, 0.0106397613572643042210903457835_rk, &
9.87500150836629658594132875969e-8_rk, -0.01108508428788088348450596993_rk, &
8.02684327858130432223945277655e-8_rk, 0.0115236078830291190096088688924_rk, &
6.52457753695726647687729943838e-8_rk, -0.0119469155908784920013184934118_rk, &
5.30346869352248747424829142388e-8_rk, 0.0123446001817378425599875205239_rk, &
4.31089676287149010856940813666e-8_rk, -0.0127033716979733978272328393362_rk, &
3.50408986534297419595820527136e-8_rk, 0.0130070826577604662202720240514_rk, &
2.84828110247310400005225049777e-8_rk, -0.0132361204791565863149064361036_rk, &
2.31521038285681732935473276987e-8_rk, 0.0133675947277084395023737938789_rk, &
1.88190663914241464458873967506e-8_rk, -0.0133754129035536630812651707772_rk, &
1.5296979594909343884130114326e-8_rk, 0.01323133939176737550471229989_rk, &
1.2434069781150546131264250961e-8_rk, -0.0129066098823208164690000886744_rk, &
1.01069685269092135880364294836e-8_rk, 0.0123749751691678909290204642266_rk, &
8.21539645521284876558501434405e-9_rk, -0.011616939032234987846492821192_rk, &
6.67784199947079723683732408319e-9_rk, 0.0106255715037222683657696978511_rk, &
5.42804890950826080536713415169e-9_rk, -0.0094131098978481832771834757995_rk, &
4.4121611392346130932390035751e-9_rk, 0.00801734065856154904541989757332_rk, &
3.58640208353164108438734664737e-9_rk, -0.00650531462665448050800581975452_rk, &
2.91518815810778459491018581983e-9_rk, 0.00497150434096454371674629246845_rk, &
2.36959543275841986175163304411e-9_rk, -0.00352733518622172601212875594832_rk, &
1.9261132422389449651773776696e-9_rk, 0.00228131323310563878615613119893_rk, &
1.56563106538804932628945935111e-9_rk, -0.00131322414655011277817932040945_rk, &
1.27261501512694207138370560793e-9_rk, 0.000651602709000465382790159887218_rk, &
1.03443845266645488573465221257e-9_rk, -0.000266340419154864206570816849064_rk, &
8.40837880769646366946062596795e-10_rk, 0.0000836618523224564821954002770626_rk, &
6.8347066943881130276640662404e-10_rk, -0.000017859320784409720649999870885_rk, &
5.55555555555555555555555555556e-10_rk, 1.93417772746275041940741982965e-6_rk /), (/2,200/))
end module os_integral_operators