/* [+MEQ MatlabEQuilibrium Toolbox+] Swiss Plasma Center EPFL Lausanne 2022. All rights reserved. */
# include "meq.h"
# ifdef SINGLE
# define FCT   sbint
# else
# define FCT   dbint
# endif

void FCT(FLT * Fi, FLT * Fx, int * k, FLT * c, int ni, int nz) {
 
 int i;
 for (i=ni; i--; k++, c+=4) {
  *Fi++ = Fx[(*k)     ] * c[0] +
          Fx[(*k)+nz  ] * c[1] +
          Fx[(*k)+nz+1] * c[2] +
          Fx[(*k)+1   ] * c[3];
 };

}
