1. Categories
2. Sphere
3. coreshellmicrogel (SASfit)
4. sasfit_coreshellmicrogel.c

coreshellmicrogel (SASfit) - sasfit_coreshellmicrogel.c

    ///////////////////////////////////////////////////
//    This is automatically genearted file       //
//    by sasfit_convert.py                       //
//    Some editting might be required            //
///////////////////////////////////////////////////

double Iq( double q, double W_CORE,  double SIGMA_CORE,  double W_SH,
           double SIGMA_SHIN,  double D,  double SIGMA_OUT,  double ETA_CORE,
           double ETA_SHELL,  double ETA_SOL);
double Fq( double q,  double W_CORE,  double SIGMA_CORE,  double W_SH,
           double SIGMA_SHIN,  double D,  double SIGMA_OUT,  double ETA_CORE,
           double ETA_SHELL,  double ETA_SOL);
double form_volume(  double W_CORE,  double SIGMA_CORE,  double W_SH,
                     double SIGMA_SHIN, double D, double SIGMA_OUT);

#define R_CORE		W_CORE+SIGMA_CORE
#define R_OUT		R_CORE+D+SIGMA_SHIN+W_SH+SIGMA_OUT
#define R_SHIN		R_CORE+D
double ampl(double q, double rad, double sigma, double contrast)
{
    double u,v,w,t,s,Vol;
    u = q*fabs(rad);
    v=q*(fabs(rad)+fabs(sigma));
    w=q*(fabs(rad)-fabs(sigma));
    t=sas_pow_5(q)*sas_pow_2(sigma);
    s=sas_pow_4(q)*sas_pow_2(sigma);
    Vol=4.*M_PI*(sas_pow_3(rad)/3.+rad*sas_pow_2(sigma)/6.0);

    if (Vol==0) return 0;
    if (sigma*q<1e-3)
    {
        if (q==0)
        {
            return Vol * contrast;
        }
        else
        {
            return Vol
                   * contrast / (sas_pow_3(rad)/3.+rad*sas_pow_2(sigma)/6.0) *
                   (					 (sin(u)-u*cos(u))/sas_pow_3(q)
                                        + sas_pow_2(sigma) * (sin(u)+u*cos(u))/(12.*q)
                                        - sas_pow_4(sigma) * ((1./360.)*q*u*cos(u)+1./120.*q*sin(u))
                   );
        }
    }
    else
    {
        if (q==0)
        {
            return Vol * contrast;
        }
        else
        {
            return Vol
                   * contrast / (sas_pow_3(rad)/3.+rad*sas_pow_2(sigma)/6.0)
                   * ( (rad/sas_pow_2(sigma)+1/sigma) * cos(v)/sas_pow_4(q)
                       +(rad/sas_pow_2(sigma)-1/sigma) * cos(w)/sas_pow_4(q)
                       -3*sin(v)/t
                       -3*sin(w)/t
                       -2*rad*cos(u)/s
                       +6*sin(u)/t
                     );
        }
    }
}
double Iq( double q, double W_CORE,  double SIGMA_CORE,  double W_SH,
           double SIGMA_SHIN,  double D,  double SIGMA_OUT,  double ETA_CORE,
           double ETA_SHELL,  double ETA_SOL)
{
    return sas_pow_2(Fq(q,W_CORE, SIGMA_CORE, W_SH, SIGMA_SHIN, D, SIGMA_OUT,
                        ETA_CORE, ETA_SHELL, ETA_SOL));
}
double Fq( double q,  double W_CORE,  double SIGMA_CORE,  double W_SH,
           double SIGMA_SHIN,  double D,  double SIGMA_OUT,  double ETA_CORE,
           double ETA_SHELL,  double ETA_SOL)
{
    return		ampl(q,R_OUT,SIGMA_OUT,(ETA_SHELL-ETA_SOL))
                -	ampl(q,R_SHIN,SIGMA_SHIN,(ETA_SHELL-ETA_SOL))
                +	ampl(q,R_CORE,SIGMA_CORE,(ETA_CORE -ETA_SOL));
}
double form_volume(  double W_CORE,  double SIGMA_CORE,  double W_SH,
                     double SIGMA_SHIN, double D, double SIGMA_OUT)
{
       return 4.*M_PI*sas_pow_3(R_OUT+SIGMA_OUT)/3.;
}

Back to Model Download