/***************************************************************************
JSPICE3 adaptation of Spice3f2 - Copyright (c) Stephen R. Whiteley 1992
Copyright 1990 Regents of the University of California.  All rights reserved.
Authors: 1988 Jaijeet S Roychowdhury
         1993 Stephen R. Whiteley
****************************************************************************/

#include "spice.h"
#include <math.h>
#define JOB char
#include "distodef.h"
#include "distoext.h"

/*
 * all subFns are local to this file so they need not be renamed to
 * the awful 7 letter standard; however, for reasons of uniformity,
 * they are being renamed, losing all readability in the process.
 * the renaming convention is as follows: 
 *   example: 3v3F1m2
 * 3v => 3 variable term xyz
 * 2F1m2 => Two F1 minus F2
 * therefore the old name would be : S3v3F1minusF2
 * for the imaginary sub functions, the v is replaced by an i
 *
 */


double
S2v2F1(cxy, r1h1x, i1h1x, r1h1y, i1h1y)

double cxy, r1h1x, i1h1x, r1h1y, i1h1y;
{
    return (cxy*(r1h1x*r1h1y - i1h1x*i1h1y));
}


double
S2i2F1(cxy, r1h1x, i1h1x, r1h1y, i1h1y)

double cxy, r1h1x, i1h1x, r1h1y, i1h1y;
{
    return (cxy*(r1h1x*i1h1y + i1h1x*r1h1y));
}


double
S2v3F1(cxy, r1h1x, i1h1x, r1h1y, i1h1y, r2h11x, i2h11x, r2h11y, i2h11y)

double cxy, r1h1x, i1h1x, r1h1y, i1h1y;
double r2h11x, i2h11x, r2h11y, i2h11y;
{
    return (cxy*(r1h1x*r2h11y - i1h1x*i2h11y + r1h1y*r2h11x - i1h1y*i2h11x));
}


double
S2i3F1(cxy, r1h1x, i1h1x, r1h1y, i1h1y, r2h11x, i2h11x, r2h11y, i2h11y)

double cxy, r1h1x, i1h1x, r1h1y, i1h1y;
double r2h11x, i2h11x, r2h11y, i2h11y;
{
    return (cxy*(r1h1x*i2h11y + i1h1x*r2h11y + r1h1y*i2h11x + i1h1y*r2h11x));
}


double
S2vF12(cxy, r1h1x, i1h1x, r1h1y, i1h1y, r1h2x, i1h2x, r1h2y, i1h2y)

double cxy, r1h1x, i1h1x, r1h1y, i1h1y;
double r1h2x, i1h2x, r1h2y, i1h2y;
{
    return (cxy*(r1h1x*r1h2y - i1h1x*i1h2y + r1h1y*r1h2x - i1h1y*i1h2x));
}


double
S2iF12(cxy, r1h1x, i1h1x, r1h1y, i1h1y, r1h2x, i1h2x, r1h2y, i1h2y)

double cxy, r1h1x, i1h1x, r1h1y, i1h1y;
double r1h2x, i1h2x, r1h2y, i1h2y;
{
    return (cxy*(r1h1x*i1h2y + i1h1x*r1h2y + r1h1y*i1h2x + i1h1y*r1h2x));
}


double
S2v2F12(cxy, r1h1x, i1h1x, r1h1y, i1h1y, r1h2x, i1h2x, r1h2y, i1h2y,
    r2h11x, i2h11x, r2h11y, i2h11y, h2f1f2x, ih2f1f2x, h2f1f2y, ih2f1f2y)

double cxy, r1h1x, i1h1x, r1h1y, i1h1y;
double r1h2x, i1h2x, r1h2y, i1h2y;
double r2h11x, i2h11x, r2h11y, i2h11y;
double h2f1f2x, ih2f1f2x, h2f1f2y, ih2f1f2y;
{
    return ( cxy * (
        2*(r1h1x*h2f1f2y - i1h1x*ih2f1f2y
          +r1h1y*h2f1f2x - i1h1y*ih2f1f2x)
        + r1h2x*r2h11y - i1h2x*i2h11y
        + r1h2y*r2h11x - i1h2y*i2h11x
        ));
}


double
S2i2F12(cxy, r1h1x, i1h1x, r1h1y, i1h1y, r1h2x, i1h2x, r1h2y, i1h2y,
    r2h11x, i2h11x, r2h11y, i2h11y, h2f1f2x, ih2f1f2x, h2f1f2y, ih2f1f2y)

double cxy, r1h1x, i1h1x, r1h1y, i1h1y;
double r1h2x, i1h2x, r1h2y, i1h2y;
double r2h11x, i2h11x, r2h11y, i2h11y;
double h2f1f2x, ih2f1f2x, h2f1f2y, ih2f1f2y;
{
    return ( cxy * (
        2*(r1h1x*ih2f1f2y + i1h1x*h2f1f2y
          +r1h1y*ih2f1f2x + i1h1y*h2f1f2x)
        + r1h2x*i2h11y + i1h2x*r2h11y
        + r1h2y*i2h11x + i1h2y*r2h11x
        ));
}


double
S3v3F1(cxyz, r1h1x, i1h1x, r1h1y, i1h1y, r1h1z, i1h1z)

double cxyz, r1h1x, i1h1x, r1h1y;
double i1h1y, r1h1z,i1h1z;
{
    return ( cxyz * (
        (r1h1x*r1h1y - i1h1x*i1h1y)*r1h1z - (i1h1x*r1h1y + r1h1x*i1h1y)*i1h1z
           ));
}


double
S3i3F1(cxyz, r1h1x, i1h1x, r1h1y, i1h1y, r1h1z, i1h1z)

double cxyz, r1h1x, i1h1x, r1h1y;
double i1h1y, r1h1z,i1h1z;
{
    return( cxyz * (
       (r1h1x*r1h1y - i1h1x*i1h1y)*i1h1z + (i1h1x*r1h1y + r1h1x*i1h1y)*r1h1z
           ));
}


double
S3v2F12(cxyz, r1h1x, i1h1x, r1h1y, i1h1y, r1h1z, i1h1z, r1h2x, i1h2x, r1h2y,
    i1h2y, r1h2z, i1h2z)

double cxyz, r1h1x, i1h1x, r1h1y, i1h1y;
double r1h1z, i1h1z, r1h2x, i1h2x, r1h2y;
double i1h2y, r1h2z, i1h2z;
{
    return ( cxyz * (
         (r1h1x*r1h1y - i1h1x*i1h1y)*r1h2z - (i1h1x*r1h1y + r1h1x*i1h1y)*i1h2z
         +
         (r1h1x*r1h1z - i1h1x*i1h1z)*r1h2y - (i1h1x*r1h1z + r1h1x*i1h1z)*i1h2y
         +
         (r1h1z*r1h1y - i1h1z*i1h1y)*r1h2x - (i1h1z*r1h1y + r1h1z*i1h1y)*i1h2x
         ));
}


double
S3i2F12(cxyz, r1h1x, i1h1x, r1h1y, i1h1y, r1h1z, i1h1z, r1h2x, i1h2x, r1h2y,
    i1h2y, r1h2z, i1h2z)

double cxyz, r1h1x, i1h1x, r1h1y, i1h1y;
double r1h1z, i1h1z, r1h2x, i1h2x, r1h2y;
double i1h2y, r1h2z, i1h2z;
{
    return ( cxyz * (
         (r1h1x*r1h1y - i1h1x*i1h1y)*i1h2z + (i1h1x*r1h1y + r1h1x*i1h1y)*r1h2z
         +
         (r1h1x*r1h1z - i1h1x*i1h1z)*i1h2y + (i1h1x*r1h1z + r1h1x*i1h1z)*r1h2y
         +
         (r1h1z*r1h1y - i1h1z*i1h1y)*i1h2x + (i1h1z*r1h1y + r1h1z*i1h1y)*r1h2x
         ));
}


/* the load functions */
/* also renamed... */

double
DFn2F1(cxx,cyy,czz,cxy,cyz,cxz, r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z)

double cxx,cyy,czz,cxy,cyz,cxz;
double r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z;
{
    double temp;
    
    temp = S2v2F1(cxx,r1h1x,i1h1x,r1h1x,i1h1x) +
           S2v2F1(cyy,r1h1y,i1h1y,r1h1y,i1h1y) +
           S2v2F1(czz,r1h1z,i1h1z,r1h1z,i1h1z) +
           S2v2F1(cxy,r1h1x,i1h1x,r1h1y,i1h1y) +
           S2v2F1(cyz,r1h1y,i1h1y,r1h1z,i1h1z) +
           S2v2F1(cxz,r1h1x,i1h1x,r1h1z,i1h1z);

     return (temp);
}


double
DFi2F1(cxx,cyy,czz,cxy,cyz,cxz, r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z)

double cxx,cyy,czz,cxy,cyz,cxz;
double r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z;
{
    double temp;
    
    temp = S2i2F1(cxx,r1h1x,i1h1x,r1h1x,i1h1x) +
           S2i2F1(cyy,r1h1y,i1h1y,r1h1y,i1h1y) +
           S2i2F1(czz,r1h1z,i1h1z,r1h1z,i1h1z) +
           S2i2F1(cxy,r1h1x,i1h1x,r1h1y,i1h1y) +
           S2i2F1(cyz,r1h1y,i1h1y,r1h1z,i1h1z) +
           S2i2F1(cxz,r1h1x,i1h1x,r1h1z,i1h1z);

    return (temp);
}


double
DFn3F1(cxx,cyy,czz,cxy,cyz,cxz,cxxx,cyyy,czzz,cxxy,cxxz,cxyy,cyyz,cxzz,cyzz,
    cxyz, r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z, r2h11x,i2h11x,r2h11y,i2h11y,
    r2h11z,i2h11z)

double cxx,cyy,czz,cxy,cyz,cxz,cxxx,cyyy,czzz,cxxy,cxxz,cxyy,cyyz,cxzz,cyzz;
double cxyz;
double r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z;
double r2h11x,i2h11x,r2h11y,i2h11y,r2h11z,i2h11z;
{
    double temp;

    temp = S2v3F1(cxx,r1h1x,i1h1x,r1h1x,i1h1x,r2h11x,i2h11x,r2h11x,i2h11x)
          +S2v3F1(cyy,r1h1y,i1h1y,r1h1y,i1h1y,r2h11y,i2h11y,r2h11y,i2h11y)
          +S2v3F1(czz,r1h1z,i1h1z,r1h1z,i1h1z,r2h11z,i2h11z,r2h11z,i2h11z);
    temp +=
          S2v3F1(cxy,r1h1x,i1h1x,r1h1y,i1h1y,r2h11x,i2h11x,r2h11y,i2h11y)
          +S2v3F1(cyz,r1h1y,i1h1y,r1h1z,i1h1z,r2h11y,i2h11y,r2h11z,i2h11z)
          +S2v3F1(cxz,r1h1x,i1h1x,r1h1z,i1h1z,r2h11x,i2h11x,r2h11z,i2h11z)
          +S3v3F1(cxxx,r1h1x,i1h1x,r1h1x,i1h1x,r1h1x,i1h1x);
    temp +=
          S3v3F1(cyyy,r1h1y,i1h1y,r1h1y,i1h1y,r1h1y,i1h1y)
          +S3v3F1(czzz,r1h1z,i1h1z,r1h1z,i1h1z,r1h1z,i1h1z)
          +S3v3F1(cxxy,r1h1x,i1h1x,r1h1x,i1h1x,r1h1y,i1h1y)
          +S3v3F1(cxxz,r1h1x,i1h1x,r1h1x,i1h1x,r1h1z,i1h1z)
          +S3v3F1(cxyy,r1h1x,i1h1x,r1h1y,i1h1y,r1h1y,i1h1y);
    temp +=
          S3v3F1(cyyz,r1h1y,i1h1y,r1h1y,i1h1y,r1h1z,i1h1z)
          +S3v3F1(cxzz,r1h1x,i1h1x,r1h1z,i1h1z,r1h1z,i1h1z)
          +S3v3F1(cyzz,r1h1y,i1h1y,r1h1z,i1h1z,r1h1z,i1h1z)
          +S3v3F1(cxyz,r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z);

    return (temp);
}


double
DFi3F1(cxx,cyy,czz,cxy,cyz,cxz,cxxx,cyyy,czzz,cxxy,cxxz,cxyy,cyyz,cxzz,cyzz,
    cxyz, r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z, r2h11x,i2h11x,r2h11y,i2h11y,
    r2h11z,i2h11z)

double cxx,cyy,czz,cxy,cyz,cxz,cxxx,cyyy,czzz,cxxy,cxxz,cxyy,cyyz,cxzz,cyzz;
double cxyz;
double r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z;
double r2h11x,i2h11x,r2h11y,i2h11y,r2h11z,i2h11z;
{
    double temp;

    temp = S2i3F1(cxx,r1h1x,i1h1x,r1h1x,i1h1x,r2h11x,i2h11x,r2h11x,i2h11x)
          +S2i3F1(cyy,r1h1y,i1h1y,r1h1y,i1h1y,r2h11y,i2h11y,r2h11y,i2h11y)
          +S2i3F1(czz,r1h1z,i1h1z,r1h1z,i1h1z,r2h11z,i2h11z,r2h11z,i2h11z)
          +S2i3F1(cxy,r1h1x,i1h1x,r1h1y,i1h1y,r2h11x,i2h11x,r2h11y,i2h11y);
    temp +=
          S2i3F1(cyz,r1h1y,i1h1y,r1h1z,i1h1z,r2h11y,i2h11y,r2h11z,i2h11z)
          +S2i3F1(cxz,r1h1x,i1h1x,r1h1z,i1h1z,r2h11x,i2h11x,r2h11z,i2h11z)
          +S3i3F1(cxxx,r1h1x,i1h1x,r1h1x,i1h1x,r1h1x,i1h1x)
          +S3i3F1(cyyy,r1h1y,i1h1y,r1h1y,i1h1y,r1h1y,i1h1y);
    temp +=
          S3i3F1(czzz,r1h1z,i1h1z,r1h1z,i1h1z,r1h1z,i1h1z)
          +S3i3F1(cxxy,r1h1x,i1h1x,r1h1x,i1h1x,r1h1y,i1h1y)
          +S3i3F1(cxxz,r1h1x,i1h1x,r1h1x,i1h1x,r1h1z,i1h1z)
          +S3i3F1(cxyy,r1h1x,i1h1x,r1h1y,i1h1y,r1h1y,i1h1y);
    temp +=
          S3i3F1(cyyz,r1h1y,i1h1y,r1h1y,i1h1y,r1h1z,i1h1z)
          +S3i3F1(cxzz,r1h1x,i1h1x,r1h1z,i1h1z,r1h1z,i1h1z)
          +S3i3F1(cyzz,r1h1y,i1h1y,r1h1z,i1h1z,r1h1z,i1h1z)
          +S3i3F1(cxyz,r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z);

    return (temp);
}


double
DFnF12(cxx,cyy,czz,cxy,cyz,cxz,r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z,
    r1h2x,i1h2x,r1h2y,i1h2y,r1h2z,i1h2z)

double cxx,cyy,czz,cxy,cyz,cxz;
double r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z;
double r1h2x,i1h2x,r1h2y,i1h2y,r1h2z,i1h2z;
{
    double temp;

    temp = S2vF12(cxx,r1h1x,i1h1x,r1h1x,i1h1x,r1h2x,i1h2x,r1h2x,i1h2x)
          +S2vF12(cyy,r1h1y,i1h1y,r1h1y,i1h1y,r1h2y,i1h2y,r1h2y,i1h2y)
          +S2vF12(czz,r1h1z,i1h1z,r1h1z,i1h1z,r1h2z,i1h2z,r1h2z,i1h2z);
    temp +=
          S2vF12(cxy,r1h1x,i1h1x,r1h1y,i1h1y,r1h2x,i1h2x,r1h2y,i1h2y)
          +S2vF12(cyz,r1h1y,i1h1y,r1h1z,i1h1z,r1h2y,i1h2y,r1h2z,i1h2z)
          +S2vF12(cxz,r1h1x,i1h1x,r1h1z,i1h1z,r1h2x,i1h2x,r1h2z,i1h2z);

    return (0.5*temp);
}


double
DFiF12(cxx,cyy,czz,cxy,cyz,cxz, r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z,
    r1h2x,i1h2x,r1h2y,i1h2y,r1h2z,i1h2z)

double cxx,cyy,czz,cxy,cyz,cxz;
double r1h1x,i1h1x,r1h1y,i1h1y,r1h1z,i1h1z;
double r1h2x,i1h2x,r1h2y,i1h2y,r1h2z,i1h2z;
{
    double temp;

    temp = S2iF12(cxx,r1h1x,i1h1x,r1h1x,i1h1x,r1h2x,i1h2x,r1h2x,i1h2x)
          +S2iF12(cyy,r1h1y,i1h1y,r1h1y,i1h1y,r1h2y,i1h2y,r1h2y,i1h2y)
          +S2iF12(czz,r1h1z,i1h1z,r1h1z,i1h1z,r1h2z,i1h2z,r1h2z,i1h2z);
    temp +=
          S2iF12(cxy,r1h1x,i1h1x,r1h1y,i1h1y,r1h2x,i1h2x,r1h2y,i1h2y)
          +S2iF12(cyz,r1h1y,i1h1y,r1h1z,i1h1z,r1h2y,i1h2y,r1h2z,i1h2z)
          +S2iF12(cxz,r1h1x,i1h1x,r1h1z,i1h1z,r1h2x,i1h2x,r1h2z,i1h2z);

    return (temp*0.5); /* divided by two to scale down */
}


double
DFn2F12(p)

/*
 * a structure because a standard C compiler can handle only
 * 32 variables.
 *
 */
DpassStr* p;
{
double temp;

    temp = S2v2F12(p->cxx,p->r1h1x, p->i1h1x,
        p->r1h1x, p->i1h1x, p->r1h2x, p->i1h2x,
        p->r1h2x, p->i1h2x, p->r2h11x, p->i2h11x,
        p->r2h11x, p->i2h11x, p->h2f1f2x, p->ih2f1f2x,
        p->h2f1f2x, p->ih2f1f2x);
    temp += S2v2F12(p->cyy, p->r1h1y, p->i1h1y,
        p->r1h1y, p->i1h1y, p->r1h2y, p->i1h2y,
        p->r1h2y, p->i1h2y, p->r2h11y, p->i2h11y,
        p->r2h11y, p->i2h11y, p->h2f1f2y, p->ih2f1f2y,
        p->h2f1f2y, p->ih2f1f2y);
    temp += S2v2F12(p->czz, p->r1h1z, p->i1h1z,
        p->r1h1z, p->i1h1z, p->r1h2z, p->i1h2z,
        p->r1h2z, p->i1h2z, p->r2h11z, p->i2h11z,
        p->r2h11z, p->i2h11z, p->h2f1f2z, p->ih2f1f2z,
        p->h2f1f2z, p->ih2f1f2z);
    temp += S2v2F12(p->cxy, p->r1h1x, p->i1h1x,
        p->r1h1y, p->i1h1y, p->r1h2x, p->i1h2x,
        p->r1h2y, p->i1h2y, p->r2h11x, p->i2h11x,
        p->r2h11y, p->i2h11y, p->h2f1f2x, p->ih2f1f2x,
        p->h2f1f2y, p->ih2f1f2y);
    temp += S2v2F12(p->cyz, p->r1h1y, p->i1h1y,
        p->r1h1z, p->i1h1z, p->r1h2y, p->i1h2y,
        p->r1h2z, p->i1h2z, p->r2h11y, p->i2h11y,
        p->r2h11z, p->i2h11z, p->h2f1f2y, p->ih2f1f2y,
        p->h2f1f2z, p->ih2f1f2z);
    temp += S2v2F12(p->cxz, p->r1h1x, p->i1h1x,
        p->r1h1z, p->i1h1z, p->r1h2x, p->i1h2x,
        p->r1h2z, p->i1h2z, p->r2h11x, p->i2h11x,
        p->r2h11z, p->i2h11z, p->h2f1f2x, p->ih2f1f2x,
        p->h2f1f2z, p->ih2f1f2z);
    temp += S3v2F12(p->cxxx, p->r1h1x,
        p->i1h1x, p->r1h1x, p->i1h1x, p->r1h1x, p->i1h1x,
        p->r1h2x, p->i1h2x, p->r1h2x, p->i1h2x, p->r1h2x, p->i1h2x)
        + S3v2F12(p->cyyy,p->r1h1y,
        p->i1h1y, p->r1h1y, p->i1h1y, p->r1h1y, p->i1h1y,
        p->r1h2y, p->i1h2y, p->r1h2y, p->i1h2y, p->r1h2y, p->i1h2y);
    temp += S3v2F12(p->czzz, p->r1h1z,
        p->i1h1z, p->r1h1z, p->i1h1z, p->r1h1z, p->i1h1z,
        p->r1h2z, p->i1h2z, p->r1h2z, p->i1h2z, p->r1h2z, p->i1h2z)
        + S3v2F12(p->cxxy,p->r1h1x,
        p->i1h1x, p->r1h1x, p->i1h1x, p->r1h1y, p->i1h1y,
        p->r1h2x, p->i1h2x, p->r1h2x, p->i1h2x, p->r1h2y, p->i1h2y);
    temp += S3v2F12(p->cxxz,p->r1h1x,
        p->i1h1x, p->r1h1x, p->i1h1x, p->r1h1z, p->i1h1z,
        p->r1h2x, p->i1h2x, p->r1h2x, p->i1h2x, p->r1h2z, p->i1h2z)
        + S3v2F12(p->cxyy,p->r1h1x,
        p->i1h1x, p->r1h1y, p->i1h1y, p->r1h1y, p->i1h1y,
        p->r1h2x, p->i1h2x, p->r1h2y, p->i1h2y, p->r1h2y, p->i1h2y);
    temp += S3v2F12(p->cyyz, p->r1h1y,
        p->i1h1y, p->r1h1y, p->i1h1y, p->r1h1z, p->i1h1z,
        p->r1h2y, p->i1h2y, p->r1h2y, p->i1h2y, p->r1h2z, p->i1h2z)
        + S3v2F12(p->cxzz, p->r1h1x,
        p->i1h1x, p->r1h1z, p->i1h1z, p->r1h1z, p->i1h1z,
        p->r1h2x, p->i1h2x, p->r1h2z, p->i1h2z, p->r1h2z, p->i1h2z);
    temp += S3v2F12(p->cyzz, p->r1h1y,
        p->i1h1y, p->r1h1z, p->i1h1z, p->r1h1z, p->i1h1z,
        p->r1h2y, p->i1h2y, p->r1h2z, p->i1h2z, p->r1h2z, p->i1h2z)
        + S3v2F12(p->cxyz, p->r1h1x,
        p->i1h1x, p->r1h1y, p->i1h1y, p->r1h1z, p->i1h1z,
        p->r1h2x, p->i1h2x, p->r1h2y, p->i1h2y, p->r1h2z, p->i1h2z);
                        
    return (temp/3.);
    /* divided by 3 to get kernel (otherwise we get 3*kernel) */
}


double
DFi2F12(p)
DpassStr* p;

{
    double temp;

    temp = S2i2F12(p->cxx, p->r1h1x, p->i1h1x,
        p->r1h1x, p->i1h1x, p->r1h2x, p->i1h2x,
        p->r1h2x, p->i1h2x, p->r2h11x, p->i2h11x,
        p->r2h11x, p->i2h11x, p->h2f1f2x, p->ih2f1f2x,
        p->h2f1f2x, p->ih2f1f2x);
    temp += S2i2F12(p->cyy, p->r1h1y, p->i1h1y,
        p->r1h1y, p->i1h1y, p->r1h2y, p->i1h2y,
        p->r1h2y, p->i1h2y, p->r2h11y, p->i2h11y,
        p->r2h11y, p->i2h11y, p->h2f1f2y, p->ih2f1f2y,
        p->h2f1f2y, p->ih2f1f2y);
    temp += S2i2F12(p->czz, p->r1h1z, p->i1h1z,
        p->r1h1z, p->i1h1z, p->r1h2z, p->i1h2z,
        p->r1h2z, p->i1h2z, p->r2h11z, p->i2h11z,
        p->r2h11z, p->i2h11z, p->h2f1f2z, p->ih2f1f2z,
        p->h2f1f2z, p->ih2f1f2z);
    temp += S2i2F12(p->cxy, p->r1h1x, p->i1h1x,
        p->r1h1y, p->i1h1y, p->r1h2x, p->i1h2x,
        p->r1h2y, p->i1h2y, p->r2h11x, p->i2h11x,
        p->r2h11y, p->i2h11y, p->h2f1f2x, p->ih2f1f2x,
        p->h2f1f2y, p->ih2f1f2y);
    temp += S2i2F12(p->cyz, p->r1h1y, p->i1h1y,
        p->r1h1z, p->i1h1z, p->r1h2y, p->i1h2y,
        p->r1h2z, p->i1h2z, p->r2h11y, p->i2h11y,
        p->r2h11z, p->i2h11z, p->h2f1f2y, p->ih2f1f2y,
        p->h2f1f2z,p->ih2f1f2z);
    temp += S2i2F12(p->cxz, p->r1h1x, p->i1h1x,
        p->r1h1z, p->i1h1z, p->r1h2x, p->i1h2x,
        p->r1h2z, p->i1h2z, p->r2h11x, p->i2h11x,
        p->r2h11z, p->i2h11z, p->h2f1f2x, p->ih2f1f2x,
        p->h2f1f2z, p->ih2f1f2z);
    temp += S3i2F12(p->cxxx, p->r1h1x,
        p->i1h1x, p->r1h1x, p->i1h1x, p->r1h1x, p->i1h1x,
        p->r1h2x, p->i1h2x, p->r1h2x, p->i1h2x, p->r1h2x, p->i1h2x);
    temp += S3i2F12(p->cyyy, p->r1h1y,
        p->i1h1y, p->r1h1y, p->i1h1y, p->r1h1y, p->i1h1y,
        p->r1h2y, p->i1h2y, p->r1h2y, p->i1h2y, p->r1h2y, p->i1h2y)
        + S3i2F12(p->czzz, p->r1h1z,
        p->i1h1z, p->r1h1z, p->i1h1z, p->r1h1z, p->i1h1z,
        p->r1h2z, p->i1h2z, p->r1h2z, p->i1h2z, p->r1h2z, p->i1h2z);
    temp += S3i2F12(p->cxxy, p->r1h1x,
        p->i1h1x, p->r1h1x, p->i1h1x, p->r1h1y, p->i1h1y,
        p->r1h2x, p->i1h2x, p->r1h2x, p->i1h2x, p->r1h2y, p->i1h2y)
        + S3i2F12(p->cxxz, p->r1h1x,
        p->i1h1x, p->r1h1x, p->i1h1x, p->r1h1z, p->i1h1z,
        p->r1h2x, p->i1h2x, p->r1h2x, p->i1h2x, p->r1h2z, p->i1h2z);
    temp += S3i2F12(p->cxyy, p->r1h1x,
        p->i1h1x, p->r1h1y, p->i1h1y, p->r1h1y, p->i1h1y,
        p->r1h2x, p->i1h2x, p->r1h2y, p->i1h2y, p->r1h2y, p->i1h2y)
        + S3i2F12(p->cyyz, p->r1h1y,
        p->i1h1y, p->r1h1y, p->i1h1y, p->r1h1z, p->i1h1z,
        p->r1h2y, p->i1h2y, p->r1h2y, p->i1h2y, p->r1h2z, p->i1h2z);
    temp += S3i2F12(p->cxzz, p->r1h1x,
        p->i1h1x, p->r1h1z, p->i1h1z, p->r1h1z, p->i1h1z,
        p->r1h2x, p->i1h2x, p->r1h2z, p->i1h2z, p->r1h2z, p->i1h2z);
    temp += S3i2F12(p->cyzz, p->r1h1y,
        p->i1h1y, p->r1h1z, p->i1h1z, p->r1h1z, p->i1h1z,
        p->r1h2y, p->i1h2y, p->r1h2z, p->i1h2z, p->r1h2z, p->i1h2z)
        + S3i2F12(p->cxyz, p->r1h1x,
        p->i1h1x, p->r1h1y, p->i1h1y, p->r1h1z, p->i1h1z,
        p->r1h2x, p->i1h2x, p->r1h2y, p->i1h2y, p->r1h2z, p->i1h2z);
                        
    return (temp/3.);
    /* divided by 3 to get kernel (otherwise we get 3*kernel) */
}


double
D1n2F1(cxx,r1h1x,i1h1x)

double cxx;
double r1h1x,i1h1x;
{
    double temp;
    
    temp = S2v2F1(cxx,r1h1x,i1h1x,r1h1x,i1h1x);
    return (temp);
}


double
D1n3F1(cxx,cxxx,r1h1x,i1h1x,r2h11x,i2h11x)

double cxx,cxxx;
double r1h1x,i1h1x;
double r2h11x,i2h11x;
{
    double temp;

    temp = S2v3F1(cxx,r1h1x,i1h1x,r1h1x,i1h1x,r2h11x,i2h11x,r2h11x,i2h11x)
          + S3v3F1(cxxx,r1h1x,i1h1x,r1h1x,i1h1x,r1h1x,i1h1x);
    return (temp);
}


double
D1nF12(cxx,r1h1x,i1h1x,r1h2x,i1h2x)

double cxx;
double r1h1x,i1h1x;
double r1h2x,i1h2x;
{
    double temp;

    temp = S2vF12(cxx,r1h1x,i1h1x,r1h1x,i1h1x,r1h2x,i1h2x,r1h2x,i1h2x);
    return (0.5*temp);
}


double
D1n2F12(cxx,cxxx,r1h1x,i1h1x,r1h2x,i1h2x,r2h11x,i2h11x,h2f1f2x,ih2f1f2x)

double cxx,cxxx;
double r1h1x,i1h1x;
double r1h2x,i1h2x;
double r2h11x,i2h11x;
double h2f1f2x,ih2f1f2x;
{
    double temp;

    temp = S2v2F12(cxx,r1h1x,i1h1x,r1h1x,i1h1x,r1h2x,i1h2x,r1h2x,i1h2x,
        r2h11x,i2h11x,r2h11x,i2h11x,h2f1f2x,ih2f1f2x,h2f1f2x,ih2f1f2x)
    + S3v2F12(cxxx,r1h1x,i1h1x,r1h1x,i1h1x,r1h1x,i1h1x,
        r1h2x,i1h2x,r1h2x,i1h2x,r1h2x,i1h2x);
                        
    return (temp/3.);
    /* divided by 3 to get kernel (otherwise we get 3*kernel) */
}


double
D1i2F1(cxx,r1h1x,i1h1x)

double cxx;
double r1h1x,i1h1x;
{
    double temp;
    
    temp = S2i2F1(cxx,r1h1x,i1h1x,r1h1x,i1h1x);
    return (temp);
}


double
D1i3F1(cxx,cxxx,r1h1x,i1h1x,r2h11x,i2h11x)

double cxx,cxxx;
double r1h1x,i1h1x;
double r2h11x,i2h11x;
{
    double temp;

    temp = S2i3F1(cxx,r1h1x,i1h1x,r1h1x,i1h1x,r2h11x,i2h11x,r2h11x,i2h11x)
        + S3i3F1(cxxx,r1h1x,i1h1x,r1h1x,i1h1x,r1h1x,i1h1x);

    return (temp);
}


double
D1iF12(cxx,r1h1x,i1h1x,r1h2x,i1h2x)

double cxx;
double r1h1x,i1h1x;
double r1h2x,i1h2x;
{
    double temp;

    temp = S2iF12(cxx,r1h1x,i1h1x,r1h1x,i1h1x,r1h2x,i1h2x,r1h2x,i1h2x);
    return (0.5*temp);
}


double
D1i2F12(cxx,cxxx,r1h1x,i1h1x,r1h2x,i1h2x,r2h11x,i2h11x,h2f1f2x,ih2f1f2x)

double cxx,cxxx;
double r1h1x,i1h1x;
double r1h2x,i1h2x;
double r2h11x,i2h11x;
double h2f1f2x,ih2f1f2x;
{
    double temp;

    temp = S2i2F12(cxx,r1h1x,i1h1x,r1h1x,i1h1x,r1h2x,i1h2x,r1h2x,i1h2x,
        r2h11x,i2h11x,r2h11x,i2h11x,h2f1f2x,ih2f1f2x,h2f1f2x,ih2f1f2x)
    + S3i2F12(cxxx,r1h1x,i1h1x,r1h1x,i1h1x,r1h1x,i1h1x,
        r1h2x,i1h2x,r1h2x,i1h2x,r1h2x,i1h2x);
                        
    return (temp/3.);
    /* divided by 3 to get kernel (otherwise we get 3*kernel) */
}