19 #ifndef _libint2_src_lib_libint_osvrrxsxsderiv_h_ 20 #define _libint2_src_lib_libint_osvrrxsxsderiv_h_ 25 #include <util_types.h> 26 #include <libint2/cgshell_ordering.h> 30 template <
int part,
int La,
int Lc,
45 static void compute(
const Libint_t* inteval,
46 LIBINT2_REALTYPE* target,
47 const LIBINT2_REALTYPE* src0,
48 const LIBINT2_REALTYPE* src1,
49 const LIBINT2_REALTYPE* src2,
50 const LIBINT2_REALTYPE* src3,
51 const LIBINT2_REALTYPE* src4,
52 const LIBINT2_REALTYPE* src5,
53 const LIBINT2_REALTYPE* src6,
54 const LIBINT2_REALTYPE* src7,
55 const LIBINT2_REALTYPE* src8,
56 const LIBINT2_REALTYPE* src9,
57 const LIBINT2_REALTYPE* src10,
58 const LIBINT2_REALTYPE* src11,
59 const LIBINT2_REALTYPE* src12,
60 const LIBINT2_REALTYPE* src13,
61 const LIBINT2_REALTYPE* src14,
62 const LIBINT2_REALTYPE* src15,
63 const LIBINT2_REALTYPE* src16,
64 const LIBINT2_REALTYPE* src17,
65 const LIBINT2_REALTYPE* src18,
66 const LIBINT2_REALTYPE* src19,
67 const LIBINT2_REALTYPE* src20,
68 const LIBINT2_REALTYPE* src21,
69 const LIBINT2_REALTYPE* src22
102 template <
int La,
int Lc,
125 LIBINT2_REALTYPE* target,
126 const LIBINT2_REALTYPE* src0,
127 const LIBINT2_REALTYPE* src1,
128 const LIBINT2_REALTYPE* src2,
129 const LIBINT2_REALTYPE* src3,
130 const LIBINT2_REALTYPE* src4,
131 const LIBINT2_REALTYPE* src5,
132 const LIBINT2_REALTYPE* src6,
133 const LIBINT2_REALTYPE* src7,
134 const LIBINT2_REALTYPE* src8,
135 const LIBINT2_REALTYPE* src9,
136 const LIBINT2_REALTYPE* src10,
137 const LIBINT2_REALTYPE* src11,
138 const LIBINT2_REALTYPE* src12,
139 const LIBINT2_REALTYPE* src13,
140 const LIBINT2_REALTYPE* src14,
141 const LIBINT2_REALTYPE* src15,
142 const LIBINT2_REALTYPE* src16,
143 const LIBINT2_REALTYPE* src17,
144 const LIBINT2_REALTYPE* src18,
145 const LIBINT2_REALTYPE* src19,
146 const LIBINT2_REALTYPE* src20,
147 const LIBINT2_REALTYPE* src21,
148 const LIBINT2_REALTYPE* src22
152 assert(not (La < 2 || Lc < 1));
154 const unsigned int veclen = vectorize ? inteval->veclen : 1;
156 const unsigned int Nc = INT_NCART(Lc);
157 const unsigned int NcV = Nc * veclen;
160 FOR_CART(ax, ay, az, La)
162 int a[3]; a[0] = ax; a[1] = ay; a[2] = az;
164 enum XYZ {x=0, y=1, z=2};
167 if (ay != 0) xyz = y;
168 if (ax != 0) xyz = x;
172 const LIBINT2_REALTYPE *PA, *WP;
175 #if LIBINT2_DEFINED(eri,PA_x) 176 if (not unit_b) PA = inteval->PA_x;
181 #if LIBINT2_DEFINED(eri,PA_y) 182 if (not unit_b) PA = inteval->PA_y;
187 #if LIBINT2_DEFINED(eri,PA_z) 188 if (not unit_b) PA = inteval->PA_z;
194 const unsigned int iam1 = INT_CARTINDEX(La-1,a[0],a[1]);
195 const unsigned int am10c0_offset = iam1 * NcV;
196 const LIBINT2_REALTYPE* src0_ptr = unit_b ? 0 : src0 + am10c0_offset;
197 const LIBINT2_REALTYPE* src1_ptr = src1 + am10c0_offset;
202 const unsigned int iam2 = INT_CARTINDEX(La-2,a[0],a[1]);
203 const unsigned int am20c0_offset = iam2 * NcV;
205 const LIBINT2_REALTYPE* src2_ptr = src2 + am20c0_offset;
206 const LIBINT2_REALTYPE* src3_ptr = src3 + am20c0_offset;
207 const LIBINT2_REALTYPE axyz = (LIBINT2_REALTYPE)a[xyz];
210 for(
unsigned int c = 0; c < Nc; ++c) {
211 for(
unsigned int v=0; v<veclen; ++v, ++cv) {
212 LIBINT2_REALTYPE value = WP[v] * src1_ptr[cv] + axyz * inteval->oo2z[v] * (src2_ptr[cv] - inteval->roz[v] * src3_ptr[cv]);
213 if (not unit_b) value += PA[v] * src0_ptr[cv];
217 #if LIBINT2_FLOP_COUNT 218 inteval->nflops[0] += (unit_b ? 6 : 8) * NcV;
224 for(
unsigned int c = 0; c < Nc; ++c) {
225 for(
unsigned int v=0; v<veclen; ++v, ++cv) {
226 LIBINT2_REALTYPE value = WP[v] * src1_ptr[cv];
228 value += PA[v] * src0_ptr[cv];
232 #if LIBINT2_FLOP_COUNT 233 inteval->nflops[0] += (unit_b ? 1 : 3) * NcV;
238 const unsigned int Ncm1 = INT_NCART(Lc-1);
239 const unsigned int Ncm1V = Ncm1 * veclen;
240 const unsigned int am10cm10_offset = iam1 * Ncm1V;
241 const LIBINT2_REALTYPE* src4_ptr = src4 + am10cm10_offset;
245 FOR_CART(cx, cy, cz, Lc-1)
247 int c[3]; c[0] = cx; c[1] = cy; c[2] = cz;
250 const unsigned int cc = INT_CARTINDEX(Lc,c[0],c[1]);
251 const unsigned int cc_offset = cc * veclen;
252 LIBINT2_REALTYPE* tptr = target + cc_offset;
253 const LIBINT2_REALTYPE cxyz = (LIBINT2_REALTYPE)c[xyz];
254 for(
unsigned int v=0; v<veclen; ++v) {
255 tptr[v] += cxyz * inteval->oo2ze[v] * src4_ptr[v];
257 #if LIBINT2_FLOP_COUNT 258 inteval->nflops[0] += 3 * veclen;
266 #define dcontrA(target,srcA,srcB,id,ecoef1,ecoef2) { \ 267 const LIBINT2_REALTYPE* srcA_ptr = srcA + am10c0_offset; \ 268 const LIBINT2_REALTYPE* srcB_ptr = srcB + am10c0_offset; \ 269 const LIBINT2_REALTYPE di = (LIBINT2_REALTYPE)id; \ 270 const LIBINT2_REALTYPE* c1 = inteval->ecoef1; \ 271 const LIBINT2_REALTYPE* c2 = inteval->ecoef2; \ 272 unsigned int cv = 0; \ 273 for(unsigned int c = 0; c < Nc; ++c) { \ 274 for(unsigned int v=0; v<veclen; ++v, ++cv) { \ 275 target[cv] -= di * (c1[v] * srcA_ptr[cv] + c2[v] * srcB_ptr[cv]); \ 280 #define dcontrB(target,srcA,srcB,id,ecoef1,ecoef2) { \ 281 const LIBINT2_REALTYPE* srcA_ptr = srcA + am10c0_offset; \ 282 const LIBINT2_REALTYPE* srcB_ptr = srcB + am10c0_offset; \ 283 const LIBINT2_REALTYPE di = (LIBINT2_REALTYPE)id; \ 284 const LIBINT2_REALTYPE* c1 = inteval->ecoef1; \ 285 const LIBINT2_REALTYPE* c2 = inteval->ecoef2; \ 286 unsigned int cv = 0; \ 287 for(unsigned int c = 0; c < Nc; ++c) { \ 288 for(unsigned int v=0; v<veclen; ++v, ++cv) { \ 289 target[cv] += di * (c1[v] * srcA_ptr[cv] - c2[v] * srcB_ptr[cv]); \ 294 #define dcontrCD(target,srcA,id,ecoef1) { \ 295 const LIBINT2_REALTYPE* srcA_ptr = srcA + am10c0_offset; \ 296 const LIBINT2_REALTYPE di = (LIBINT2_REALTYPE)id; \ 297 const LIBINT2_REALTYPE* c1 = inteval->ecoef1; \ 298 unsigned int cv = 0; \ 299 for(unsigned int c = 0; c < Nc; ++c) { \ 300 for(unsigned int v=0; v<veclen; ++v, ++cv) { \ 301 target[cv] += di * c1[v] * srcA_ptr[cv]; \ 307 #if LIBINT2_DEFINED(any,rho12_over_alpha1) && LIBINT2_DEFINED(any,alpha1_rho_over_zeta2) 308 if (Da_x > 0 && xyz == x){
309 dcontrA(target,src5,src6,Da_x,rho12_over_alpha1,alpha1_rho_over_zeta2);
310 #if LIBINT2_FLOP_COUNT 311 inteval->nflops[0] += 5 * NcV;
314 if (Da_y > 0 && xyz == y){
315 dcontrA(target,src11,src12,Da_y,rho12_over_alpha1,alpha1_rho_over_zeta2);
316 #if LIBINT2_FLOP_COUNT 317 inteval->nflops[0] += 5 * NcV;
320 if (Da_z > 0 && xyz == z){
321 dcontrA(target,src17,src18,Da_z,rho12_over_alpha1,alpha1_rho_over_zeta2);
322 #if LIBINT2_FLOP_COUNT 323 inteval->nflops[0] += 5 * NcV;
328 #if LIBINT2_DEFINED(any,rho12_over_alpha1) && LIBINT2_DEFINED(any,alpha2_rho_over_zeta2) 329 if (Db_x > 0 && xyz == x){
330 dcontrB(target,src7,src8,Db_x,rho12_over_alpha1,alpha2_rho_over_zeta2);
331 #if LIBINT2_FLOP_COUNT 332 inteval->nflops[0] += 5 * NcV;
335 if (Db_y > 0 && xyz == y){
336 dcontrB(target,src13,src14,Db_y,rho12_over_alpha1,alpha2_rho_over_zeta2);
337 #if LIBINT2_FLOP_COUNT 338 inteval->nflops[0] += 5 * NcV;
341 if (Db_z > 0 && xyz == z){
342 dcontrB(target,src19,src20,Db_z,rho12_over_alpha1,alpha2_rho_over_zeta2);
343 #if LIBINT2_FLOP_COUNT 344 inteval->nflops[0] += 5 * NcV;
349 #if LIBINT2_DEFINED(any,alpha3_over_zetapluseta) 350 if (Dc_x > 0 && xyz == x){
351 dcontrCD(target,src9,Dc_x,alpha3_over_zetapluseta);
352 #if LIBINT2_FLOP_COUNT 353 inteval->nflops[0] += 3 * NcV;
356 if (Dc_y > 0 && xyz == y){
357 dcontrCD(target,src15,Dc_y,alpha3_over_zetapluseta);
358 #if LIBINT2_FLOP_COUNT 359 inteval->nflops[0] += 3 * NcV;
362 if (Dc_z > 0 && xyz == z){
363 dcontrCD(target,src21,Dc_z,alpha3_over_zetapluseta);
364 #if LIBINT2_FLOP_COUNT 365 inteval->nflops[0] += 3 * NcV;
370 #if LIBINT2_DEFINED(any,alpha4_over_zetapluseta) 371 if (Dd_x > 0 && xyz == x){
372 dcontrCD(target,src10,Dd_x,alpha4_over_zetapluseta);
373 #if LIBINT2_FLOP_COUNT 374 inteval->nflops[0] += 3 * NcV;
377 if (Dd_y > 0 && xyz == y){
378 dcontrCD(target,src16,Dd_y,alpha4_over_zetapluseta);
379 #if LIBINT2_FLOP_COUNT 380 inteval->nflops[0] += 3 * NcV;
383 if (Dd_z > 0 && xyz == z){
384 dcontrCD(target,src22,Dd_z,alpha4_over_zetapluseta);
385 #if LIBINT2_FLOP_COUNT 386 inteval->nflops[0] += 3 * NcV;
403 template <
int part,
int La,
int Lc,
418 static void compute(
const Libint_t* inteval,
419 LIBINT2_REALTYPE* target,
420 const LIBINT2_REALTYPE* src1,
421 const LIBINT2_REALTYPE* src4,
422 const LIBINT2_REALTYPE* src5,
423 const LIBINT2_REALTYPE* src6,
424 const LIBINT2_REALTYPE* src7,
425 const LIBINT2_REALTYPE* src8,
426 const LIBINT2_REALTYPE* src9,
427 const LIBINT2_REALTYPE* src10,
428 const LIBINT2_REALTYPE* src11,
429 const LIBINT2_REALTYPE* src12,
430 const LIBINT2_REALTYPE* src13,
431 const LIBINT2_REALTYPE* src14,
432 const LIBINT2_REALTYPE* src15,
433 const LIBINT2_REALTYPE* src16,
434 const LIBINT2_REALTYPE* src17,
435 const LIBINT2_REALTYPE* src18,
436 const LIBINT2_REALTYPE* src19,
437 const LIBINT2_REALTYPE* src20,
438 const LIBINT2_REALTYPE* src21,
439 const LIBINT2_REALTYPE* src22
469 template <
int La,
int Lc,
490 LIBINT2_REALTYPE* target,
491 const LIBINT2_REALTYPE* src1,
492 const LIBINT2_REALTYPE* src4,
493 const LIBINT2_REALTYPE* src5,
494 const LIBINT2_REALTYPE* src6,
495 const LIBINT2_REALTYPE* src7,
496 const LIBINT2_REALTYPE* src8,
497 const LIBINT2_REALTYPE* src9,
498 const LIBINT2_REALTYPE* src10,
499 const LIBINT2_REALTYPE* src11,
500 const LIBINT2_REALTYPE* src12,
501 const LIBINT2_REALTYPE* src13,
502 const LIBINT2_REALTYPE* src14,
503 const LIBINT2_REALTYPE* src15,
504 const LIBINT2_REALTYPE* src16,
505 const LIBINT2_REALTYPE* src17,
506 const LIBINT2_REALTYPE* src18,
507 const LIBINT2_REALTYPE* src19,
508 const LIBINT2_REALTYPE* src20,
509 const LIBINT2_REALTYPE* src21,
510 const LIBINT2_REALTYPE* src22
514 assert(not (La < 1 || Lc < 1));
516 const unsigned int veclen = vectorize ? inteval->veclen : 1;
518 const unsigned int Nc = INT_NCART(Lc);
519 const unsigned int NcV = Nc * veclen;
522 FOR_CART(ax, ay, az, La)
524 int a[3]; a[0] = ax; a[1] = ay; a[2] = az;
526 enum XYZ {x=0, y=1, z=2};
529 if (ay != 0) xyz = y;
530 if (ax != 0) xyz = x;
534 const LIBINT2_REALTYPE *WP;
547 const unsigned int iam1 = INT_CARTINDEX(La-1,a[0],a[1]);
548 const unsigned int am10c0_offset = iam1 * NcV;
549 const LIBINT2_REALTYPE* src1_ptr = src1 + am10c0_offset;
553 for(
unsigned int c = 0; c < Nc; ++c) {
554 for(
unsigned int v=0; v<veclen; ++v, ++cv) {
555 target[cv] = WP[v] * src1_ptr[cv];
558 #if LIBINT2_FLOP_COUNT 559 inteval->nflops[0] += NcV;
564 const unsigned int Ncm1 = INT_NCART(Lc-1);
565 const unsigned int Ncm1V = Ncm1 * veclen;
566 const unsigned int am10cm10_offset = iam1 * Ncm1V;
567 const LIBINT2_REALTYPE* src4_ptr = src4 + am10cm10_offset;
571 FOR_CART(cx, cy, cz, Lc-1)
573 int c[3]; c[0] = cx; c[1] = cy; c[2] = cz;
576 const unsigned int cc = INT_CARTINDEX(Lc,c[0],c[1]);
577 const unsigned int cc_offset = cc * veclen;
578 LIBINT2_REALTYPE* tptr = target + cc_offset;
579 const LIBINT2_REALTYPE cxyz = (LIBINT2_REALTYPE)c[xyz];
580 for(
unsigned int v=0; v<veclen; ++v) {
581 tptr[v] += cxyz * inteval->oo2ze[v] * src4_ptr[v];
583 #if LIBINT2_FLOP_COUNT 584 inteval->nflops[0] += 3 * veclen;
592 #define dcontrA(target,srcA,srcB,id,ecoef1,ecoef2) { \ 593 const LIBINT2_REALTYPE* srcA_ptr = srcA + am10c0_offset; \ 594 const LIBINT2_REALTYPE* srcB_ptr = srcB + am10c0_offset; \ 595 const LIBINT2_REALTYPE di = (LIBINT2_REALTYPE)id; \ 596 const LIBINT2_REALTYPE* c1 = inteval->ecoef1; \ 597 const LIBINT2_REALTYPE* c2 = inteval->ecoef2; \ 598 unsigned int cv = 0; \ 599 for(unsigned int c = 0; c < Nc; ++c) { \ 600 for(unsigned int v=0; v<veclen; ++v, ++cv) { \ 601 target[cv] -= di * (c1[v] * srcA_ptr[cv] + c2[v] * srcB_ptr[cv]); \ 606 #define dcontrB(target,srcA,srcB,id,ecoef1,ecoef2) { \ 607 const LIBINT2_REALTYPE* srcA_ptr = srcA + am10c0_offset; \ 608 const LIBINT2_REALTYPE* srcB_ptr = srcB + am10c0_offset; \ 609 const LIBINT2_REALTYPE di = (LIBINT2_REALTYPE)id; \ 610 const LIBINT2_REALTYPE* c1 = inteval->ecoef1; \ 611 const LIBINT2_REALTYPE* c2 = inteval->ecoef2; \ 612 unsigned int cv = 0; \ 613 for(unsigned int c = 0; c < Nc; ++c) { \ 614 for(unsigned int v=0; v<veclen; ++v, ++cv) { \ 615 target[cv] += di * (c1[v] * srcA_ptr[cv] - c2[v] * srcB_ptr[cv]); \ 620 #define dcontrCD(target,srcA,id,ecoef1) { \ 621 const LIBINT2_REALTYPE* srcA_ptr = srcA + am10c0_offset; \ 622 const LIBINT2_REALTYPE di = (LIBINT2_REALTYPE)id; \ 623 const LIBINT2_REALTYPE* c1 = inteval->ecoef1; \ 624 unsigned int cv = 0; \ 625 for(unsigned int c = 0; c < Nc; ++c) { \ 626 for(unsigned int v=0; v<veclen; ++v, ++cv) { \ 627 target[cv] += di * c1[v] * srcA_ptr[cv]; \ 633 #if LIBINT2_DEFINED(any,rho12_over_alpha1) && LIBINT2_DEFINED(any,alpha1_rho_over_zeta2) 634 if (Da_x > 0 && xyz == x){
635 dcontrA(target,src5,src6,Da_x,rho12_over_alpha1,alpha1_rho_over_zeta2);
636 #if LIBINT2_FLOP_COUNT 637 inteval->nflops[0] += 5 * NcV;
640 if (Da_y > 0 && xyz == y){
641 dcontrA(target,src11,src12,Da_y,rho12_over_alpha1,alpha1_rho_over_zeta2);
642 #if LIBINT2_FLOP_COUNT 643 inteval->nflops[0] += 5 * NcV;
646 if (Da_z > 0 && xyz == z){
647 dcontrA(target,src17,src18,Da_z,rho12_over_alpha1,alpha1_rho_over_zeta2);
648 #if LIBINT2_FLOP_COUNT 649 inteval->nflops[0] += 5 * NcV;
654 #if LIBINT2_DEFINED(any,rho12_over_alpha1) && LIBINT2_DEFINED(any,alpha2_rho_over_zeta2) 655 if (Db_x > 0 && xyz == x){
656 dcontrB(target,src7,src8,Db_x,rho12_over_alpha1,alpha2_rho_over_zeta2);
657 #if LIBINT2_FLOP_COUNT 658 inteval->nflops[0] += 5 * NcV;
661 if (Db_y > 0 && xyz == y){
662 dcontrB(target,src13,src14,Db_y,rho12_over_alpha1,alpha2_rho_over_zeta2);
663 #if LIBINT2_FLOP_COUNT 664 inteval->nflops[0] += 5 * NcV;
667 if (Db_z > 0 && xyz == z){
668 dcontrB(target,src19,src20,Db_z,rho12_over_alpha1,alpha2_rho_over_zeta2);
669 #if LIBINT2_FLOP_COUNT 670 inteval->nflops[0] += 5 * NcV;
675 #if LIBINT2_DEFINED(any,alpha3_over_zetapluseta) 676 if (Dc_x > 0 && xyz == x){
677 dcontrCD(target,src9,Dc_x,alpha3_over_zetapluseta);
678 #if LIBINT2_FLOP_COUNT 679 inteval->nflops[0] += 3 * NcV;
682 if (Dc_y > 0 && xyz == y){
683 dcontrCD(target,src15,Dc_y,alpha3_over_zetapluseta);
684 #if LIBINT2_FLOP_COUNT 685 inteval->nflops[0] += 3 * NcV;
688 if (Dc_z > 0 && xyz == z){
689 dcontrCD(target,src21,Dc_z,alpha3_over_zetapluseta);
690 #if LIBINT2_FLOP_COUNT 691 inteval->nflops[0] += 3 * NcV;
696 #if LIBINT2_DEFINED(any,alpha4_over_zetapluseta) 697 if (Dd_x > 0 && xyz == x){
698 dcontrCD(target,src10,Dd_x,alpha4_over_zetapluseta);
699 #if LIBINT2_FLOP_COUNT 700 inteval->nflops[0] += 3 * NcV;
703 if (Dd_y > 0 && xyz == y){
704 dcontrCD(target,src16,Dd_y,alpha4_over_zetapluseta);
705 #if LIBINT2_FLOP_COUNT 706 inteval->nflops[0] += 3 * NcV;
709 if (Dd_z > 0 && xyz == z){
710 dcontrCD(target,src22,Dd_z,alpha4_over_zetapluseta);
711 #if LIBINT2_FLOP_COUNT 712 inteval->nflops[0] += 3 * NcV;
730 #endif // header guard libint2::OSVRR_xs_xs_deriv< 0, La, Lc, Da_x, Da_y, Da_z, Db_x, Db_y, Db_z, Dc_x, Dc_y, Dc_z, Dd_x, Dd_y, Dd_z, unit_b, vectorize >::compute static void compute(const Libint_t *inteval, LIBINT2_REALTYPE *target, const LIBINT2_REALTYPE *src0, const LIBINT2_REALTYPE *src1, const LIBINT2_REALTYPE *src2, const LIBINT2_REALTYPE *src3, const LIBINT2_REALTYPE *src4, const LIBINT2_REALTYPE *src5, const LIBINT2_REALTYPE *src6, const LIBINT2_REALTYPE *src7, const LIBINT2_REALTYPE *src8, const LIBINT2_REALTYPE *src9, const LIBINT2_REALTYPE *src10, const LIBINT2_REALTYPE *src11, const LIBINT2_REALTYPE *src12, const LIBINT2_REALTYPE *src13, const LIBINT2_REALTYPE *src14, const LIBINT2_REALTYPE *src15, const LIBINT2_REALTYPE *src16, const LIBINT2_REALTYPE *src17, const LIBINT2_REALTYPE *src18, const LIBINT2_REALTYPE *src19, const LIBINT2_REALTYPE *src20, const LIBINT2_REALTYPE *src21, const LIBINT2_REALTYPE *src22)
Definition: OSVRR_xs_xs_deriv.h:124
Definition: OSVRR_xs_xs_deriv.h:44
Defaults definitions for various parameters assumed by Libint.
Definition: algebra.cc:23
libint2::OSAVRR_xs_xs_deriv< 0, La, Lc, Da_x, Da_y, Da_z, Db_x, Db_y, Db_z, Dc_x, Dc_y, Dc_z, Dd_x, Dd_y, Dd_z, vectorize >::compute static void compute(const Libint_t *inteval, LIBINT2_REALTYPE *target, const LIBINT2_REALTYPE *src1, const LIBINT2_REALTYPE *src4, const LIBINT2_REALTYPE *src5, const LIBINT2_REALTYPE *src6, const LIBINT2_REALTYPE *src7, const LIBINT2_REALTYPE *src8, const LIBINT2_REALTYPE *src9, const LIBINT2_REALTYPE *src10, const LIBINT2_REALTYPE *src11, const LIBINT2_REALTYPE *src12, const LIBINT2_REALTYPE *src13, const LIBINT2_REALTYPE *src14, const LIBINT2_REALTYPE *src15, const LIBINT2_REALTYPE *src16, const LIBINT2_REALTYPE *src17, const LIBINT2_REALTYPE *src18, const LIBINT2_REALTYPE *src19, const LIBINT2_REALTYPE *src20, const LIBINT2_REALTYPE *src21, const LIBINT2_REALTYPE *src22)
Definition: OSVRR_xs_xs_deriv.h:489
Definition: test_eri_rys.cc:46
Definition: OSVRR_xs_xs_deriv.h:417