def _epd_function(coeffs, fsv, fdv, fkv, xcc, ycc, bgv, bge, iha, izd):
return (coeffs[<NUM_LIT:0>]*fsv*fsv +<EOL>coeffs[<NUM_LIT:1>]*fsv +<EOL>coeffs[<NUM_LIT:2>]*fdv*fdv +<EOL>coeffs[<NUM_LIT:3>]*fdv +<EOL>coeffs[<NUM_LIT:4>]*fkv*fkv +<EOL>coeffs[<NUM_LIT:5>]*fkv +<EOL>coeffs[<NUM_LIT:6>] +<EOL>coeffs[<NUM_LIT:7>]*fsv*fdv +<EOL>coeffs[<NUM_LIT:8>]*fsv*fkv +<EOL>coeffs[<NUM_LIT:9>]*fdv*fkv +<EOL>coeffs[<NUM_LIT:10>]*np.sin(<NUM_LIT:2>*pi_value*xcc) +<EOL>coeffs[<NUM_LIT:11>]*np.cos(<NUM_LIT:2>*pi_value*xcc) +<EOL>coeffs[<NUM_LIT:12>]*np.sin(<NUM_LIT:2>*pi_value*ycc) +<EOL>coeffs[<NUM_LIT>]*np.cos(<NUM_LIT:2>*pi_value*ycc) +<EOL>coeffs[<NUM_LIT>]*np.sin(<NUM_LIT:4>*pi_value*xcc) +<EOL>coeffs[<NUM_LIT:15>]*np.cos(<NUM_LIT:4>*pi_value*xcc) +<EOL>coeffs[<NUM_LIT:16>]*np.sin(<NUM_LIT:4>*pi_value*ycc) +<EOL>coeffs[<NUM_LIT>]*np.cos(<NUM_LIT:4>*pi_value*ycc) +<EOL>coeffs[<NUM_LIT>]*bgv +<EOL>coeffs[<NUM_LIT>]*bge +<EOL>coeffs[<NUM_LIT:20>]*iha +<EOL>coeffs[<NUM_LIT>]*izd)<EOL>
This is the EPD function to fit using a smoothed mag-series.