Пример #1
0
two_division *divide_network_in_two\
(mod_matrix *mod_mat, eigen_pair *leading_eigen_pair, int use_improve) {
	elem_vector *s;
	two_division *result;
	int i;
	if (IS_POSITIVE(leading_eigen_pair->value)) {
		if ((s = get_s_vector_for(leading_eigen_pair->vector)) == NULL) {
			return NULL;
		}
		if ((result = allocate_two_division(s)) == NULL) {
			free_elem_vector(s);
			return NULL;
		}
		if (use_improve != 0) {
			if(improve_network_division(mod_mat, result) == 0) {
				free_two_division(result); /*also frees s !!! */
				return NULL;
			}
		}
		if (calculate_modularity_score(mod_mat, s, &result->quality) == 0) {
			free_two_division(result); /*also frees s !!! */
			return NULL;
		}
		if (IS_POSITIVE(result->quality)) {
			return result;
		}
		free_two_division(result); /*also frees s !!! */
	}
	/*return the trivial division*/
	if ((s = allocate_elem_vector(mod_mat->A_g->n)) == NULL) {
		return NULL;
	}
	if ((result = allocate_two_division(s)) == NULL) {
		free_elem_vector(s);
		return NULL;
	}
	for(i=0; i<result->s_vector->n; i++) {
		result->s_vector->values[i] = 1;
	}
	return result;
}
Пример #2
0
/* get the s vector corresponding to 'vector'
 * meaning s[i] = sign(vector[i])
 */
elem_vector *get_s_vector_for(elem_vector *vector) {
	elem_vector *s;
	int i;
	if ((s = allocate_elem_vector(vector->n)) == NULL) {
		return NULL;
	}
	for(i=0; i<s->n; i++) {
		if (IS_POSITIVE(vector->values[i])) {
			s->values[i] = 1;
		} else {
			s->values[i] = -1;
		}
	}
	return s;
}
Пример #3
0
 bool OBResidue::GetAminoAcidProperty(int property) const
 {
   switch(property)
     {
     case OBAminoAcidProperty::ACIDIC:
       return IS_ACIDIC(_aakey)      != 0;
     case OBAminoAcidProperty::ACYCLIC:
       return IS_ACYCLIC(_aakey)     != 0;
     case OBAminoAcidProperty::ALIPHATIC:
       return IS_ALIPHATIC(_aakey)   != 0;
     case OBAminoAcidProperty::AROMATIC:
       return IS_AROMATIC(_aakey)    != 0;
     case OBAminoAcidProperty::BASIC:
       return IS_BASIC(_aakey)       != 0;
     case OBAminoAcidProperty::BURIED:
       return IS_BURIED(_aakey)      != 0;
     case OBAminoAcidProperty::CHARGED:
       return IS_CHARGED(_aakey)     != 0;
     case OBAminoAcidProperty::CYCLIC:
       return IS_CYCLIC(_aakey)      != 0;
     case OBAminoAcidProperty::HYDROPHOBIC:
       return IS_HYDROPHOBIC(_aakey) != 0;
     case OBAminoAcidProperty::LARGE:
       return IS_LARGE(_aakey)       != 0;
     case OBAminoAcidProperty::MEDIUM:
       return IS_MEDIUM(_aakey)      != 0;
     case OBAminoAcidProperty::NEGATIVE:
       return IS_NEGATIVE(_aakey)    != 0;
     case OBAminoAcidProperty::NEUTRAL:
       return IS_NEUTRAL(_aakey)     != 0;
     case OBAminoAcidProperty::POLAR:
       return IS_POLAR(_aakey)       != 0;
     case OBAminoAcidProperty::POSITIVE:
       return IS_POSITIVE(_aakey)    != 0;
     case OBAminoAcidProperty::SMALL:
       return IS_SMALL(_aakey)       != 0;
     case OBAminoAcidProperty::SURFACE:
       return IS_SURFACE(_aakey)     != 0;
     default:
       return false;
     }
 }
Пример #4
0
void frictionContact3D_AlartCurnierABGenerated(
  double rn,
  double rt1,
  double rt2,
  double un,
  double ut1,
  double ut2,
  double mu,
  double rhon,
  double rhot1,
  double rhot2,
  double *result)
{
    double x1;
    int x2;
    x1=-rhon*un + rn;
    x2=x1 <= 0;
    int x3;
    x3=x1 > 0;
    double x4;
    double x5;
    double x6;
    double x7;
    double x8;
    double x9;
    double x10;
    int x11;
    int x12;
    int x13;
    int x14;
    int x15;
    x4=rhot1*ut1;
    x5=-rt1 + x4;
    x6=rhot2*ut2;
    x7=-rt2 + x6;
    x8=x5*x5 + x7*x7;
    x9=(assert(IS_POSITIVE(x8)), sqrt(x8));
    x10=mu*x1;
    x11=x9 <= x10;
    x12=x9 <= 0;
    x13=x9 > x10;
    x14=x9 > 0;
    x15=(x11 || x2) && (x11 || x12 || x2) && (x11 || x13 || x2) && (x11 || x2 || x3) && (x12 || x2 || x3) && (x11 || x12 || x13 || x2) && (x11 || x12 || x2 || x3) && (x11 || x13 || x14 || x2) && (x11 || x14 || x2 || x3);
    int x20;
    int x21;
    double x16;
    double x17;
    double x18;
    double x19;
    x20=x13 && x3 || x14 && x2;
    x21=x20 && x3;
    int x22;
    x22=x11 && x3 || x12 && x2;
    int x23;
    x23=x2 && x20;
    /* Assignment result[0, 0]=Piecewise((0, x2), (rhon, x3)) */
    if (x2)
    {
        DEBUG_PRINT("Case (x2) is True.\n");

        /* Assignment result[0, 0]=0 */
        result[0] = 0;
    }
    else if (x3)
    {
        DEBUG_PRINT("Case (x3) is True.\n");

        /* Assignment result[0, 0]=rhon */
        result[0] = rhon;
    }

    /* Assignment result[1, 0]=Piecewise((0, x15), (rhon*x19, x21)) */
    if (x15)
    {
        DEBUG_PRINT("Case (x15) is True.\n");

        /* Assignment result[1, 0]=0 */
        result[1] = 0;
    }
    else if (x21)
    {
        DEBUG_PRINT("Case (x21) is True.\n");
        x16=rt1 - x4;
        x17=1.0/(assert(IS_NOT_ZERO(x9)), x9);
        x18=mu*x17;
        x19=x16*x18;

        /* Assignment result[1, 0]=rhon*x19 */
        result[1] = rhon*x19;
    }

    /* Assignment result[2, 0]=Piecewise((0, x15), (rhon*x31, x21)) */
    double x30;
    double x31;if (x15)
    {
        DEBUG_PRINT("Case (x15) is True.\n");

        /* Assignment result[2, 0]=0 */
        result[2] = 0;
    }
    else if (x21)
    {
        DEBUG_PRINT("Case (x21) is True.\n");
        x17=1.0/(assert(IS_NOT_ZERO(x9)), x9);
        x18=mu*x17;
        x30=rt2 - x6;
        x31=x18*x30;

        /* Assignment result[2, 0]=rhon*x31 */
        result[2] = rhon*x31;
    }

    /* Assignment result[0, 1]=0 */
    result[3] = 0;

    /* Assignment result[1, 1]=Piecewise((rhot1, x22), (0, x23), (rhot1*x24 + rhot1*x27, x21)) */
    double x24;
    double x25;
    double x26;
    double x27;if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");

        /* Assignment result[1, 1]=rhot1 */
        result[4] = rhot1;
    }
    else if (x23)
    {
        DEBUG_PRINT("Case (x23) is True.\n");

        /* Assignment result[1, 1]=0 */
        result[4] = 0;
    }
    else if (x21)
    {
        DEBUG_PRINT("Case (x21) is True.\n");
        x16=rt1 - x4;
        x17=1.0/(assert(IS_NOT_ZERO(x9)), x9);
        x24=x10*x17;
        x25=pow(x8, -3.0/2.0);
        x26=mu*x1*x16*x25;
        x27=x26*x5;

        /* Assignment result[1, 1]=rhot1*x24 + rhot1*x27 */
        result[4] = rhot1*x24 + rhot1*x27;
    }

    /* Assignment result[2, 1]=Piecewise((0, x15), (rhot1*x33, x21)) */
    double x32;
    double x33;if (x15)
    {
        DEBUG_PRINT("Case (x15) is True.\n");

        /* Assignment result[2, 1]=0 */
        result[5] = 0;
    }
    else if (x21)
    {
        DEBUG_PRINT("Case (x21) is True.\n");
        x25=pow(x8, -3.0/2.0);
        x30=rt2 - x6;
        x32=mu*x1*x25*x30;
        x33=x32*x5;

        /* Assignment result[2, 1]=rhot1*x33 */
        result[5] = rhot1*x33;
    }

    /* Assignment result[0, 2]=0 */
    result[6] = 0;

    /* Assignment result[1, 2]=Piecewise((0, x15), (rhot2*x28, x21)) */
    double x28;if (x15)
    {
        DEBUG_PRINT("Case (x15) is True.\n");

        /* Assignment result[1, 2]=0 */
        result[7] = 0;
    }
    else if (x21)
    {
        DEBUG_PRINT("Case (x21) is True.\n");
        x16=rt1 - x4;
        x25=pow(x8, -3.0/2.0);
        x26=mu*x1*x16*x25;
        x28=x26*x7;

        /* Assignment result[1, 2]=rhot2*x28 */
        result[7] = rhot2*x28;
    }

    /* Assignment result[2, 2]=Piecewise((rhot2, x22), (0, x23), (rhot2*x24 + rhot2*x34, x21)) */
    double x34;if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");

        /* Assignment result[2, 2]=rhot2 */
        result[8] = rhot2;
    }
    else if (x23)
    {
        DEBUG_PRINT("Case (x23) is True.\n");

        /* Assignment result[2, 2]=0 */
        result[8] = 0;
    }
    else if (x21)
    {
        DEBUG_PRINT("Case (x21) is True.\n");
        x17=1.0/(assert(IS_NOT_ZERO(x9)), x9);
        x24=x10*x17;
        x25=pow(x8, -3.0/2.0);
        x30=rt2 - x6;
        x32=mu*x1*x25*x30;
        x34=x32*x7;

        /* Assignment result[2, 2]=rhot2*x24 + rhot2*x34 */
        result[8] = rhot2*x24 + rhot2*x34;
    }

    /* Assignment result[0, 3]=Piecewise((1, x2), (0, x3)) */
    if (x2)
    {
        DEBUG_PRINT("Case (x2) is True.\n");

        /* Assignment result[0, 3]=1 */
        result[9] = 1;
    }
    else if (x3)
    {
        DEBUG_PRINT("Case (x3) is True.\n");

        /* Assignment result[0, 3]=0 */
        result[9] = 0;
    }

    /* Assignment result[1, 3]=Piecewise((0, x15), (-x19, x21)) */
    if (x15)
    {
        DEBUG_PRINT("Case (x15) is True.\n");

        /* Assignment result[1, 3]=0 */
        result[10] = 0;
    }
    else if (x21)
    {
        DEBUG_PRINT("Case (x21) is True.\n");
        x16=rt1 - x4;
        x17=1.0/(assert(IS_NOT_ZERO(x9)), x9);
        x18=mu*x17;
        x19=x16*x18;

        /* Assignment result[1, 3]=-x19 */
        result[10] = -x19;
    }

    /* Assignment result[2, 3]=Piecewise((0, x15), (-x31, x21)) */
    if (x15)
    {
        DEBUG_PRINT("Case (x15) is True.\n");

        /* Assignment result[2, 3]=0 */
        result[11] = 0;
    }
    else if (x21)
    {
        DEBUG_PRINT("Case (x21) is True.\n");
        x17=1.0/(assert(IS_NOT_ZERO(x9)), x9);
        x18=mu*x17;
        x30=rt2 - x6;
        x31=x18*x30;

        /* Assignment result[2, 3]=-x31 */
        result[11] = -x31;
    }

    /* Assignment result[0, 4]=0 */
    result[12] = 0;

    /* Assignment result[1, 4]=Piecewise((0, x22), (1, x23), (-x27 + x29, x21)) */
    double x29;if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");

        /* Assignment result[1, 4]=0 */
        result[13] = 0;
    }
    else if (x23)
    {
        DEBUG_PRINT("Case (x23) is True.\n");

        /* Assignment result[1, 4]=1 */
        result[13] = 1;
    }
    else if (x21)
    {
        DEBUG_PRINT("Case (x21) is True.\n");
        x16=rt1 - x4;
        x17=1.0/(assert(IS_NOT_ZERO(x9)), x9);
        x24=x10*x17;
        x25=pow(x8, -3.0/2.0);
        x26=mu*x1*x16*x25;
        x27=x26*x5;
        x29=-x24 + 1;

        /* Assignment result[1, 4]=-x27 + x29 */
        result[13] = -x27 + x29;
    }

    /* Assignment result[2, 4]=Piecewise((0, x15), (-x33, x21)) */
    if (x15)
    {
        DEBUG_PRINT("Case (x15) is True.\n");

        /* Assignment result[2, 4]=0 */
        result[14] = 0;
    }
    else if (x21)
    {
        DEBUG_PRINT("Case (x21) is True.\n");
        x25=pow(x8, -3.0/2.0);
        x30=rt2 - x6;
        x32=mu*x1*x25*x30;
        x33=x32*x5;

        /* Assignment result[2, 4]=-x33 */
        result[14] = -x33;
    }

    /* Assignment result[0, 5]=0 */
    result[15] = 0;

    /* Assignment result[1, 5]=Piecewise((0, x15), (-x28, x21)) */
    if (x15)
    {
        DEBUG_PRINT("Case (x15) is True.\n");

        /* Assignment result[1, 5]=0 */
        result[16] = 0;
    }
    else if (x21)
    {
        DEBUG_PRINT("Case (x21) is True.\n");
        x16=rt1 - x4;
        x25=pow(x8, -3.0/2.0);
        x26=mu*x1*x16*x25;
        x28=x26*x7;

        /* Assignment result[1, 5]=-x28 */
        result[16] = -x28;
    }

    /* Assignment result[2, 5]=Piecewise((0, x22), (1, x23), (x29 - x34, x21)) */
    if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");

        /* Assignment result[2, 5]=0 */
        result[17] = 0;
    }
    else if (x23)
    {
        DEBUG_PRINT("Case (x23) is True.\n");

        /* Assignment result[2, 5]=1 */
        result[17] = 1;
    }
    else if (x21)
    {
        DEBUG_PRINT("Case (x21) is True.\n");
        x17=1.0/(assert(IS_NOT_ZERO(x9)), x9);
        x24=x10*x17;
        x25=pow(x8, -3.0/2.0);
        x29=-x24 + 1;
        x30=rt2 - x6;
        x32=mu*x1*x25*x30;
        x34=x32*x7;

        /* Assignment result[2, 5]=x29 - x34 */
        result[17] = x29 - x34;
    }
}
Пример #5
0
void frictionContact3D_AlartCurnierFABGenerated(
  double rn,
  double rt1,
  double rt2,
  double un,
  double ut1,
  double ut2,
  double mu,
  double rhon,
  double rhot1,
  double rhot2,
  double *result)
{
    double x1;
    double x2;
    int x3;
    x1=rhon*un;
    x2=rn - x1;
    x3=x2 <= 0;
    int x4;
    x4=x2 > 0;
    double x5;
    double x6;
    double x7;
    double x8;
    double x9;
    double x10;
    int x11;
    double x12;
    int x13;
    int x14;
    x5=rhot1*ut1;
    x6=-rt1 + x5;
    x7=rhot2*ut2;
    x8=-rt2 + x7;
    x9=x6*x6 + x8*x8;
    x10=(assert(IS_POSITIVE(x9)), sqrt(x9));
    x11=x10 <= 0;
    x12=mu*x2;
    x13=x10 <= x12;
    x14=x11 && x3 || x13 && x4;
    int x15;
    int x16;
    int x17;
    int x18;
    x15=x10 > 0;
    x16=x10 > x12;
    x17=x15 && x3 || x16 && x4;
    x18=x17 && x3;
    int x22;
    double x19;
    double x20;
    double x21;
    x22=x17 && x4;
    int x23;
    x23=(x13 || x3) && (x11 || x13 || x3) && (x11 || x3 || x4) && (x13 || x16 || x3) && (x13 || x3 || x4) && (x11 || x13 || x16 || x3) && (x11 || x13 || x3 || x4) && (x13 || x15 || x16 || x3) && (x13 || x15 || x3 || x4);
    /* Assignment result[0, 0]=Piecewise((rn, x3), (x1, x4)) */
    if (x3)
    {
        DEBUG_PRINT("Case (x3) is True.\n");

        /* Assignment result[0, 0]=rn */
        result[0] = rn;
    }
    else if (x4)
    {
        DEBUG_PRINT("Case (x4) is True.\n");

        /* Assignment result[0, 0]=x1 */
        result[0] = x1;
    }

    /* Assignment result[1, 0]=Piecewise((x5, x14), (rt1, x18), (rt1 - x19*x21, x22)) */
    if (x14)
    {
        DEBUG_PRINT("Case (x14) is True.\n");

        /* Assignment result[1, 0]=x5 */
        result[1] = x5;
    }
    else if (x18)
    {
        DEBUG_PRINT("Case (x18) is True.\n");

        /* Assignment result[1, 0]=rt1 */
        result[1] = rt1;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x19=rt1 - x5;
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x21=x12*x20;

        /* Assignment result[1, 0]=rt1 - x19*x21 */
        result[1] = rt1 - x19*x21;
    }

    /* Assignment result[2, 0]=Piecewise((x7, x14), (rt2, x18), (rt2 - x21*x31, x22)) */
    double x31;if (x14)
    {
        DEBUG_PRINT("Case (x14) is True.\n");

        /* Assignment result[2, 0]=x7 */
        result[2] = x7;
    }
    else if (x18)
    {
        DEBUG_PRINT("Case (x18) is True.\n");

        /* Assignment result[2, 0]=rt2 */
        result[2] = rt2;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x21=x12*x20;
        x31=rt2 - x7;

        /* Assignment result[2, 0]=rt2 - x21*x31 */
        result[2] = rt2 - x21*x31;
    }

    /* Assignment result[0, 1]=Piecewise((0, x3), (rhon, x4)) */
    if (x3)
    {
        DEBUG_PRINT("Case (x3) is True.\n");

        /* Assignment result[0, 1]=0 */
        result[3] = 0;
    }
    else if (x4)
    {
        DEBUG_PRINT("Case (x4) is True.\n");

        /* Assignment result[0, 1]=rhon */
        result[3] = rhon;
    }

    /* Assignment result[1, 1]=Piecewise((0, x23), (rhon*x25, x22)) */
    double x24;
    double x25;if (x23)
    {
        DEBUG_PRINT("Case (x23) is True.\n");

        /* Assignment result[1, 1]=0 */
        result[4] = 0;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x19=rt1 - x5;
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x24=mu*x20;
        x25=x19*x24;

        /* Assignment result[1, 1]=rhon*x25 */
        result[4] = rhon*x25;
    }

    /* Assignment result[2, 1]=Piecewise((0, x23), (rhon*x32, x22)) */
    double x32;if (x23)
    {
        DEBUG_PRINT("Case (x23) is True.\n");

        /* Assignment result[2, 1]=0 */
        result[5] = 0;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x24=mu*x20;
        x31=rt2 - x7;
        x32=x24*x31;

        /* Assignment result[2, 1]=rhon*x32 */
        result[5] = rhon*x32;
    }

    /* Assignment result[0, 2]=0 */
    result[6] = 0;

    /* Assignment result[1, 2]=Piecewise((rhot1, x14), (0, x18), (rhot1*x21 + rhot1*x28, x22)) */
    double x26;
    double x27;
    double x28;if (x14)
    {
        DEBUG_PRINT("Case (x14) is True.\n");

        /* Assignment result[1, 2]=rhot1 */
        result[7] = rhot1;
    }
    else if (x18)
    {
        DEBUG_PRINT("Case (x18) is True.\n");

        /* Assignment result[1, 2]=0 */
        result[7] = 0;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x19=rt1 - x5;
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x21=x12*x20;
        x26=pow(x9, -3.0/2.0);
        x27=mu*x19*x2*x26;
        x28=x27*x6;

        /* Assignment result[1, 2]=rhot1*x21 + rhot1*x28 */
        result[7] = rhot1*x21 + rhot1*x28;
    }

    /* Assignment result[2, 2]=Piecewise((0, x23), (rhot1*x34, x22)) */
    double x33;
    double x34;if (x23)
    {
        DEBUG_PRINT("Case (x23) is True.\n");

        /* Assignment result[2, 2]=0 */
        result[8] = 0;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x26=pow(x9, -3.0/2.0);
        x31=rt2 - x7;
        x33=mu*x2*x26*x31;
        x34=x33*x6;

        /* Assignment result[2, 2]=rhot1*x34 */
        result[8] = rhot1*x34;
    }

    /* Assignment result[0, 3]=0 */
    result[9] = 0;

    /* Assignment result[1, 3]=Piecewise((0, x23), (rhot2*x29, x22)) */
    double x29;if (x23)
    {
        DEBUG_PRINT("Case (x23) is True.\n");

        /* Assignment result[1, 3]=0 */
        result[10] = 0;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x19=rt1 - x5;
        x26=pow(x9, -3.0/2.0);
        x27=mu*x19*x2*x26;
        x29=x27*x8;

        /* Assignment result[1, 3]=rhot2*x29 */
        result[10] = rhot2*x29;
    }

    /* Assignment result[2, 3]=Piecewise((rhot2, x14), (0, x18), (rhot2*x21 + rhot2*x35, x22)) */
    double x35;if (x14)
    {
        DEBUG_PRINT("Case (x14) is True.\n");

        /* Assignment result[2, 3]=rhot2 */
        result[11] = rhot2;
    }
    else if (x18)
    {
        DEBUG_PRINT("Case (x18) is True.\n");

        /* Assignment result[2, 3]=0 */
        result[11] = 0;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x21=x12*x20;
        x26=pow(x9, -3.0/2.0);
        x31=rt2 - x7;
        x33=mu*x2*x26*x31;
        x35=x33*x8;

        /* Assignment result[2, 3]=rhot2*x21 + rhot2*x35 */
        result[11] = rhot2*x21 + rhot2*x35;
    }

    /* Assignment result[0, 4]=Piecewise((1, x3), (0, x4)) */
    if (x3)
    {
        DEBUG_PRINT("Case (x3) is True.\n");

        /* Assignment result[0, 4]=1 */
        result[12] = 1;
    }
    else if (x4)
    {
        DEBUG_PRINT("Case (x4) is True.\n");

        /* Assignment result[0, 4]=0 */
        result[12] = 0;
    }

    /* Assignment result[1, 4]=Piecewise((0, x23), (-x25, x22)) */
    if (x23)
    {
        DEBUG_PRINT("Case (x23) is True.\n");

        /* Assignment result[1, 4]=0 */
        result[13] = 0;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x19=rt1 - x5;
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x24=mu*x20;
        x25=x19*x24;

        /* Assignment result[1, 4]=-x25 */
        result[13] = -x25;
    }

    /* Assignment result[2, 4]=Piecewise((0, x23), (-x32, x22)) */
    if (x23)
    {
        DEBUG_PRINT("Case (x23) is True.\n");

        /* Assignment result[2, 4]=0 */
        result[14] = 0;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x24=mu*x20;
        x31=rt2 - x7;
        x32=x24*x31;

        /* Assignment result[2, 4]=-x32 */
        result[14] = -x32;
    }

    /* Assignment result[0, 5]=0 */
    result[15] = 0;

    /* Assignment result[1, 5]=Piecewise((0, x14), (1, x18), (-x28 + x30, x22)) */
    double x30;if (x14)
    {
        DEBUG_PRINT("Case (x14) is True.\n");

        /* Assignment result[1, 5]=0 */
        result[16] = 0;
    }
    else if (x18)
    {
        DEBUG_PRINT("Case (x18) is True.\n");

        /* Assignment result[1, 5]=1 */
        result[16] = 1;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x19=rt1 - x5;
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x21=x12*x20;
        x26=pow(x9, -3.0/2.0);
        x27=mu*x19*x2*x26;
        x28=x27*x6;
        x30=-x21 + 1;

        /* Assignment result[1, 5]=-x28 + x30 */
        result[16] = -x28 + x30;
    }

    /* Assignment result[2, 5]=Piecewise((0, x23), (-x34, x22)) */
    if (x23)
    {
        DEBUG_PRINT("Case (x23) is True.\n");

        /* Assignment result[2, 5]=0 */
        result[17] = 0;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x26=pow(x9, -3.0/2.0);
        x31=rt2 - x7;
        x33=mu*x2*x26*x31;
        x34=x33*x6;

        /* Assignment result[2, 5]=-x34 */
        result[17] = -x34;
    }

    /* Assignment result[0, 6]=0 */
    result[18] = 0;

    /* Assignment result[1, 6]=Piecewise((0, x23), (-x29, x22)) */
    if (x23)
    {
        DEBUG_PRINT("Case (x23) is True.\n");

        /* Assignment result[1, 6]=0 */
        result[19] = 0;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x19=rt1 - x5;
        x26=pow(x9, -3.0/2.0);
        x27=mu*x19*x2*x26;
        x29=x27*x8;

        /* Assignment result[1, 6]=-x29 */
        result[19] = -x29;
    }

    /* Assignment result[2, 6]=Piecewise((0, x14), (1, x18), (x30 - x35, x22)) */
    if (x14)
    {
        DEBUG_PRINT("Case (x14) is True.\n");

        /* Assignment result[2, 6]=0 */
        result[20] = 0;
    }
    else if (x18)
    {
        DEBUG_PRINT("Case (x18) is True.\n");

        /* Assignment result[2, 6]=1 */
        result[20] = 1;
    }
    else if (x22)
    {
        DEBUG_PRINT("Case (x22) is True.\n");
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x21=x12*x20;
        x26=pow(x9, -3.0/2.0);
        x30=-x21 + 1;
        x31=rt2 - x7;
        x33=mu*x2*x26*x31;
        x35=x33*x8;

        /* Assignment result[2, 6]=x30 - x35 */
        result[20] = x30 - x35;
    }
}
Пример #6
0
void frictionContact3D_AlartCurnierFGenerated(
  double rn,
  double rt1,
  double rt2,
  double un,
  double ut1,
  double ut2,
  double mu,
  double rhon,
  double rhot1,
  double rhot2,
  double *result)
{
    double x1;
    double x2;
    int x3;
    x1=rhon*un;
    x2=rn - x1;
    x3=x2 <= 0;
    int x4;
    x4=x2 > 0;
    double x5;
    double x6;
    double x7;
    double x8;
    int x9;
    x5=rhot1*ut1;
    x6=rhot2*ut2;
    x7=(assert(IS_POSITIVE((-rt1 + x5)*(-rt1 + x5) + (-rt2 + x6)*(-rt2 + x6))), sqrt((-rt1 + x5)*(-rt1 + x5) + (-rt2 + x6)*(-rt2 + x6)));
    x8=mu*x2;
    x9=x3 && x7 <= 0 || x4 && x7 <= x8;
    int x10;
    int x11;
    x10=x3 && x7 > 0 || x4 && x7 > x8;
    x11=x10 && x3;
    int x13;
    double x12;
    x13=x10 && x4;
    /* Assignment result[0, 0]=Piecewise((rn, x3), (x1, x4)) */
    if (x3)
    {
        DEBUG_PRINT("Case (x3) is True.\n");

        /* Assignment result[0, 0]=rn */
        result[0] = rn;
    }
    else if (x4)
    {
        DEBUG_PRINT("Case (x4) is True.\n");

        /* Assignment result[0, 0]=x1 */
        result[0] = x1;
    }

    /* Assignment result[1, 0]=Piecewise((x5, x9), (rt1, x11), (rt1 - x12*(rt1 - x5), x13)) */
    if (x9)
    {
        DEBUG_PRINT("Case (x9) is True.\n");

        /* Assignment result[1, 0]=x5 */
        result[1] = x5;
    }
    else if (x11)
    {
        DEBUG_PRINT("Case (x11) is True.\n");

        /* Assignment result[1, 0]=rt1 */
        result[1] = rt1;
    }
    else if (x13)
    {
        DEBUG_PRINT("Case (x13) is True.\n");
        x12=(assert(IS_NOT_ZERO(x7)), mu*x2/x7);

        /* Assignment result[1, 0]=rt1 - x12*(rt1 - x5) */
        result[1] = rt1 - x12*(rt1 - x5);
    }

    /* Assignment result[2, 0]=Piecewise((x6, x9), (rt2, x11), (rt2 - x12*(rt2 - x6), x13)) */
    if (x9)
    {
        DEBUG_PRINT("Case (x9) is True.\n");

        /* Assignment result[2, 0]=x6 */
        result[2] = x6;
    }
    else if (x11)
    {
        DEBUG_PRINT("Case (x11) is True.\n");

        /* Assignment result[2, 0]=rt2 */
        result[2] = rt2;
    }
    else if (x13)
    {
        DEBUG_PRINT("Case (x13) is True.\n");
        x12=(assert(IS_NOT_ZERO(x7)), mu*x2/x7);

        /* Assignment result[2, 0]=rt2 - x12*(rt2 - x6) */
        result[2] = rt2 - x12*(rt2 - x6);
    }
}
Пример #7
0
void frictionContact3D_AlartCurnierJeanMoreauABGenerated(
  double rn,
  double rt1,
  double rt2,
  double un,
  double ut1,
  double ut2,
  double mu,
  double rhon,
  double rhot1,
  double rhot2,
  double *result)
{
    double x1;
    int x2;
    x1=-rhon*un + rn;
    x2=x1 <= 0;
    int x3;
    x3=x1 > 0;
    int x4;
    double x5;
    double x6;
    double x7;
    double x8;
    double x9;
    double x10;
    int x11;
    int x12;
    double x13;
    int x14;
    int x15;
    x4=rn <= 0;
    x5=rhot1*ut1;
    x6=-rt1 + x5;
    x7=rhot2*ut2;
    x8=-rt2 + x7;
    x9=x6*x6 + x8*x8;
    x10=(assert(IS_POSITIVE(x9)), sqrt(x9));
    x11=x10 <= 0;
    x12=rn > 0;
    x13=mu*rn;
    x14=x10 <= x13;
    x15=x11 && x4 || x12 && x14;
    int x16;
    int x17;
    int x18;
    int x19;
    x16=x10 > 0;
    x17=x10 > x13;
    x18=x12 && x17 || x16 && x4;
    x19=x18 && x4;
    int x26;
    double x20;
    double x21;
    double x22;
    double x23;
    double x24;
    double x25;
    x26=x12 && x18;
    int x27;
    x27=(x14 || x4) && (x11 || x12 || x4) && (x11 || x14 || x4) && (x12 || x14 || x4) && (x14 || x17 || x4) && (x11 || x12 || x14 || x4) && (x11 || x14 || x17 || x4) && (x12 || x14 || x16 || x4) && (x14 || x16 || x17 || x4);
    /* Assignment result[0, 0]=Piecewise((0, x2), (rhon, x3)) */
    if (x2)
    {
        DEBUG_PRINT("Case (x2) is True.\n");

        /* Assignment result[0, 0]=0 */
        result[0] = 0;
    }
    else if (x3)
    {
        DEBUG_PRINT("Case (x3) is True.\n");

        /* Assignment result[0, 0]=rhon */
        result[0] = rhon;
    }

    /* Assignment result[1, 0]=0 */
    result[1] = 0;

    /* Assignment result[2, 0]=0 */
    result[2] = 0;

    /* Assignment result[0, 1]=0 */
    result[3] = 0;

    /* Assignment result[1, 1]=Piecewise((rhot1, x15), (0, x19), (rhot1*x21 + rhot1*x25, x26)) */
    if (x15)
    {
        DEBUG_PRINT("Case (x15) is True.\n");

        /* Assignment result[1, 1]=rhot1 */
        result[4] = rhot1;
    }
    else if (x19)
    {
        DEBUG_PRINT("Case (x19) is True.\n");

        /* Assignment result[1, 1]=0 */
        result[4] = 0;
    }
    else if (x26)
    {
        DEBUG_PRINT("Case (x26) is True.\n");
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x21=x13*x20;
        x22=rt1 - x5;
        x23=pow(x9, -3.0/2.0);
        x24=mu*rn*x22*x23;
        x25=x24*x6;

        /* Assignment result[1, 1]=rhot1*x21 + rhot1*x25 */
        result[4] = rhot1*x21 + rhot1*x25;
    }

    /* Assignment result[2, 1]=Piecewise((0, x27), (rhot1*x33, x26)) */
    double x31;
    double x32;
    double x33;if (x27)
    {
        DEBUG_PRINT("Case (x27) is True.\n");

        /* Assignment result[2, 1]=0 */
        result[5] = 0;
    }
    else if (x26)
    {
        DEBUG_PRINT("Case (x26) is True.\n");
        x23=pow(x9, -3.0/2.0);
        x31=rt2 - x7;
        x32=mu*rn*x23*x31;
        x33=x32*x6;

        /* Assignment result[2, 1]=rhot1*x33 */
        result[5] = rhot1*x33;
    }

    /* Assignment result[0, 2]=0 */
    result[6] = 0;

    /* Assignment result[1, 2]=Piecewise((0, x27), (rhot2*x28, x26)) */
    double x28;if (x27)
    {
        DEBUG_PRINT("Case (x27) is True.\n");

        /* Assignment result[1, 2]=0 */
        result[7] = 0;
    }
    else if (x26)
    {
        DEBUG_PRINT("Case (x26) is True.\n");
        x22=rt1 - x5;
        x23=pow(x9, -3.0/2.0);
        x24=mu*rn*x22*x23;
        x28=x24*x8;

        /* Assignment result[1, 2]=rhot2*x28 */
        result[7] = rhot2*x28;
    }

    /* Assignment result[2, 2]=Piecewise((rhot2, x15), (0, x19), (rhot2*x21 + rhot2*x34, x26)) */
    double x34;if (x15)
    {
        DEBUG_PRINT("Case (x15) is True.\n");

        /* Assignment result[2, 2]=rhot2 */
        result[8] = rhot2;
    }
    else if (x19)
    {
        DEBUG_PRINT("Case (x19) is True.\n");

        /* Assignment result[2, 2]=0 */
        result[8] = 0;
    }
    else if (x26)
    {
        DEBUG_PRINT("Case (x26) is True.\n");
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x21=x13*x20;
        x23=pow(x9, -3.0/2.0);
        x31=rt2 - x7;
        x32=mu*rn*x23*x31;
        x34=x32*x8;

        /* Assignment result[2, 2]=rhot2*x21 + rhot2*x34 */
        result[8] = rhot2*x21 + rhot2*x34;
    }

    /* Assignment result[0, 3]=Piecewise((1, x2), (0, x3)) */
    if (x2)
    {
        DEBUG_PRINT("Case (x2) is True.\n");

        /* Assignment result[0, 3]=1 */
        result[9] = 1;
    }
    else if (x3)
    {
        DEBUG_PRINT("Case (x3) is True.\n");

        /* Assignment result[0, 3]=0 */
        result[9] = 0;
    }

    /* Assignment result[1, 3]=Piecewise((0, x27), (-x22*x29, x26)) */
    double x29;if (x27)
    {
        DEBUG_PRINT("Case (x27) is True.\n");

        /* Assignment result[1, 3]=0 */
        result[10] = 0;
    }
    else if (x26)
    {
        DEBUG_PRINT("Case (x26) is True.\n");
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x22=rt1 - x5;
        x29=mu*x20;

        /* Assignment result[1, 3]=-x22*x29 */
        result[10] = -x22*x29;
    }

    /* Assignment result[2, 3]=Piecewise((0, x27), (-x29*x31, x26)) */
    if (x27)
    {
        DEBUG_PRINT("Case (x27) is True.\n");

        /* Assignment result[2, 3]=0 */
        result[11] = 0;
    }
    else if (x26)
    {
        DEBUG_PRINT("Case (x26) is True.\n");
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x29=mu*x20;
        x31=rt2 - x7;

        /* Assignment result[2, 3]=-x29*x31 */
        result[11] = -x29*x31;
    }

    /* Assignment result[0, 4]=0 */
    result[12] = 0;

    /* Assignment result[1, 4]=Piecewise((0, x15), (1, x19), (-x25 + x30, x26)) */
    double x30;if (x15)
    {
        DEBUG_PRINT("Case (x15) is True.\n");

        /* Assignment result[1, 4]=0 */
        result[13] = 0;
    }
    else if (x19)
    {
        DEBUG_PRINT("Case (x19) is True.\n");

        /* Assignment result[1, 4]=1 */
        result[13] = 1;
    }
    else if (x26)
    {
        DEBUG_PRINT("Case (x26) is True.\n");
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x21=x13*x20;
        x22=rt1 - x5;
        x23=pow(x9, -3.0/2.0);
        x24=mu*rn*x22*x23;
        x25=x24*x6;
        x30=-x21 + 1;

        /* Assignment result[1, 4]=-x25 + x30 */
        result[13] = -x25 + x30;
    }

    /* Assignment result[2, 4]=Piecewise((0, x27), (-x33, x26)) */
    if (x27)
    {
        DEBUG_PRINT("Case (x27) is True.\n");

        /* Assignment result[2, 4]=0 */
        result[14] = 0;
    }
    else if (x26)
    {
        DEBUG_PRINT("Case (x26) is True.\n");
        x23=pow(x9, -3.0/2.0);
        x31=rt2 - x7;
        x32=mu*rn*x23*x31;
        x33=x32*x6;

        /* Assignment result[2, 4]=-x33 */
        result[14] = -x33;
    }

    /* Assignment result[0, 5]=0 */
    result[15] = 0;

    /* Assignment result[1, 5]=Piecewise((0, x27), (-x28, x26)) */
    if (x27)
    {
        DEBUG_PRINT("Case (x27) is True.\n");

        /* Assignment result[1, 5]=0 */
        result[16] = 0;
    }
    else if (x26)
    {
        DEBUG_PRINT("Case (x26) is True.\n");
        x22=rt1 - x5;
        x23=pow(x9, -3.0/2.0);
        x24=mu*rn*x22*x23;
        x28=x24*x8;

        /* Assignment result[1, 5]=-x28 */
        result[16] = -x28;
    }

    /* Assignment result[2, 5]=Piecewise((0, x15), (1, x19), (x30 - x34, x26)) */
    if (x15)
    {
        DEBUG_PRINT("Case (x15) is True.\n");

        /* Assignment result[2, 5]=0 */
        result[17] = 0;
    }
    else if (x19)
    {
        DEBUG_PRINT("Case (x19) is True.\n");

        /* Assignment result[2, 5]=1 */
        result[17] = 1;
    }
    else if (x26)
    {
        DEBUG_PRINT("Case (x26) is True.\n");
        x20=1.0/(assert(IS_NOT_ZERO(x10)), x10);
        x21=x13*x20;
        x23=pow(x9, -3.0/2.0);
        x30=-x21 + 1;
        x31=rt2 - x7;
        x32=mu*rn*x23*x31;
        x34=x32*x8;

        /* Assignment result[2, 5]=x30 - x34 */
        result[17] = x30 - x34;
    }
}
Пример #8
0
void frictionContact3D_AlartCurnierJeanMoreauFGenerated(
  double rn,
  double rt1,
  double rt2,
  double un,
  double ut1,
  double ut2,
  double mu,
  double rhon,
  double rhot1,
  double rhot2,
  double *result)
{
    double x1;
    double x2;
    x1=rhon*un;
    x2=rn - x1;
    double x3;
    int x4;
    double x5;
    double x6;
    int x7;
    double x8;
    int x9;
    x3=rhot1*ut1;
    x4=rn <= 0;
    x5=rhot2*ut2;
    x6=(assert(IS_POSITIVE((-rt1 + x3)*(-rt1 + x3) + (-rt2 + x5)*(-rt2 + x5))), sqrt((-rt1 + x3)*(-rt1 + x3) + (-rt2 + x5)*(-rt2 + x5)));
    x7=rn > 0;
    x8=mu*rn;
    x9=x4 && x6 <= 0 || x7 && x6 <= x8;
    int x10;
    int x11;
    x10=x4 && x6 > 0 || x7 && x6 > x8;
    x11=x10 && x4;
    int x13;
    double x12;
    x13=x10 && x7;
    /* Assignment result[0, 0]=Piecewise((rn, x2 <= 0), (x1, x2 > 0)) */
    if (x2 <= 0)
    {
        DEBUG_PRINT("Case (x2 <= 0) is True.\n");

        /* Assignment result[0, 0]=rn */
        result[0] = rn;
    }
    else if (x2 > 0)
    {
        DEBUG_PRINT("Case (x2 > 0) is True.\n");

        /* Assignment result[0, 0]=x1 */
        result[0] = x1;
    }

    /* Assignment result[1, 0]=Piecewise((x3, x9), (rt1, x11), (rt1 - x12*(rt1 - x3), x13)) */
    if (x9)
    {
        DEBUG_PRINT("Case (x9) is True.\n");

        /* Assignment result[1, 0]=x3 */
        result[1] = x3;
    }
    else if (x11)
    {
        DEBUG_PRINT("Case (x11) is True.\n");

        /* Assignment result[1, 0]=rt1 */
        result[1] = rt1;
    }
    else if (x13)
    {
        DEBUG_PRINT("Case (x13) is True.\n");
        x12=(assert(IS_NOT_ZERO(x6)), mu*rn/x6);

        /* Assignment result[1, 0]=rt1 - x12*(rt1 - x3) */
        result[1] = rt1 - x12*(rt1 - x3);
    }

    /* Assignment result[2, 0]=Piecewise((x5, x9), (rt2, x11), (rt2 - x12*(rt2 - x5), x13)) */
    if (x9)
    {
        DEBUG_PRINT("Case (x9) is True.\n");

        /* Assignment result[2, 0]=x5 */
        result[2] = x5;
    }
    else if (x11)
    {
        DEBUG_PRINT("Case (x11) is True.\n");

        /* Assignment result[2, 0]=rt2 */
        result[2] = rt2;
    }
    else if (x13)
    {
        DEBUG_PRINT("Case (x13) is True.\n");
        x12=(assert(IS_NOT_ZERO(x6)), mu*rn/x6);

        /* Assignment result[2, 0]=rt2 - x12*(rt2 - x5) */
        result[2] = rt2 - x12*(rt2 - x5);
    }
}
Пример #9
0
void frictionContact3D_AlartCurnierJeanMoreauFABGenerated(
  double rn,
  double rt1,
  double rt2,
  double un,
  double ut1,
  double ut2,
  double mu,
  double rhon,
  double rhot1,
  double rhot2,
  double *result)
{
    double x1;
    double x2;
    int x3;
    x1=rhon*un;
    x2=rn - x1;
    x3=x2 <= 0;
    int x4;
    x4=x2 > 0;
    double x5;
    int x6;
    double x7;
    double x8;
    double x9;
    double x10;
    double x11;
    int x12;
    int x13;
    double x14;
    int x15;
    int x16;
    x5=rhot1*ut1;
    x6=rn <= 0;
    x7=-rt1 + x5;
    x8=rhot2*ut2;
    x9=-rt2 + x8;
    x10=x7*x7 + x9*x9;
    x11=(assert(IS_POSITIVE(x10)), sqrt(x10));
    x12=x11 <= 0;
    x13=rn > 0;
    x14=mu*rn;
    x15=x11 <= x14;
    x16=x12 && x6 || x13 && x15;
    int x17;
    int x18;
    int x19;
    int x20;
    x17=x11 > 0;
    x18=x11 > x14;
    x19=x13 && x18 || x17 && x6;
    x20=x19 && x6;
    int x24;
    double x21;
    double x22;
    double x23;
    x24=x13 && x19;
    int x28;
    x28=(x15 || x6) && (x12 || x13 || x6) && (x12 || x15 || x6) && (x13 || x15 || x6) && (x15 || x18 || x6) && (x12 || x13 || x15 || x6) && (x12 || x15 || x18 || x6) && (x13 || x15 || x17 || x6) && (x15 || x17 || x18 || x6);
    /* Assignment result[0, 0]=Piecewise((rn, x3), (x1, x4)) */
    if (x3)
    {
        DEBUG_PRINT("Case (x3) is True.\n");

        /* Assignment result[0, 0]=rn */
        result[0] = rn;
    }
    else if (x4)
    {
        DEBUG_PRINT("Case (x4) is True.\n");

        /* Assignment result[0, 0]=x1 */
        result[0] = x1;
    }

    /* Assignment result[1, 0]=Piecewise((x5, x16), (rt1, x20), (rt1 - x21*x23, x24)) */
    if (x16)
    {
        DEBUG_PRINT("Case (x16) is True.\n");

        /* Assignment result[1, 0]=x5 */
        result[1] = x5;
    }
    else if (x20)
    {
        DEBUG_PRINT("Case (x20) is True.\n");

        /* Assignment result[1, 0]=rt1 */
        result[1] = rt1;
    }
    else if (x24)
    {
        DEBUG_PRINT("Case (x24) is True.\n");
        x21=rt1 - x5;
        x22=1.0/(assert(IS_NOT_ZERO(x11)), x11);
        x23=x14*x22;

        /* Assignment result[1, 0]=rt1 - x21*x23 */
        result[1] = rt1 - x21*x23;
    }

    /* Assignment result[2, 0]=Piecewise((x8, x16), (rt2, x20), (rt2 - x23*x32, x24)) */
    double x32;if (x16)
    {
        DEBUG_PRINT("Case (x16) is True.\n");

        /* Assignment result[2, 0]=x8 */
        result[2] = x8;
    }
    else if (x20)
    {
        DEBUG_PRINT("Case (x20) is True.\n");

        /* Assignment result[2, 0]=rt2 */
        result[2] = rt2;
    }
    else if (x24)
    {
        DEBUG_PRINT("Case (x24) is True.\n");
        x22=1.0/(assert(IS_NOT_ZERO(x11)), x11);
        x23=x14*x22;
        x32=rt2 - x8;

        /* Assignment result[2, 0]=rt2 - x23*x32 */
        result[2] = rt2 - x23*x32;
    }

    /* Assignment result[0, 1]=Piecewise((0, x3), (rhon, x4)) */
    if (x3)
    {
        DEBUG_PRINT("Case (x3) is True.\n");

        /* Assignment result[0, 1]=0 */
        result[3] = 0;
    }
    else if (x4)
    {
        DEBUG_PRINT("Case (x4) is True.\n");

        /* Assignment result[0, 1]=rhon */
        result[3] = rhon;
    }

    /* Assignment result[1, 1]=0 */
    result[4] = 0;

    /* Assignment result[2, 1]=0 */
    result[5] = 0;

    /* Assignment result[0, 2]=0 */
    result[6] = 0;

    /* Assignment result[1, 2]=Piecewise((rhot1, x16), (0, x20), (rhot1*x23 + rhot1*x27, x24)) */
    double x25;
    double x26;
    double x27;if (x16)
    {
        DEBUG_PRINT("Case (x16) is True.\n");

        /* Assignment result[1, 2]=rhot1 */
        result[7] = rhot1;
    }
    else if (x20)
    {
        DEBUG_PRINT("Case (x20) is True.\n");

        /* Assignment result[1, 2]=0 */
        result[7] = 0;
    }
    else if (x24)
    {
        DEBUG_PRINT("Case (x24) is True.\n");
        x21=rt1 - x5;
        x22=1.0/(assert(IS_NOT_ZERO(x11)), x11);
        x23=x14*x22;
        x25=pow(x10, -3.0/2.0);
        x26=mu*rn*x21*x25;
        x27=x26*x7;

        /* Assignment result[1, 2]=rhot1*x23 + rhot1*x27 */
        result[7] = rhot1*x23 + rhot1*x27;
    }

    /* Assignment result[2, 2]=Piecewise((0, x28), (rhot1*x34, x24)) */
    double x33;
    double x34;if (x28)
    {
        DEBUG_PRINT("Case (x28) is True.\n");

        /* Assignment result[2, 2]=0 */
        result[8] = 0;
    }
    else if (x24)
    {
        DEBUG_PRINT("Case (x24) is True.\n");
        x25=pow(x10, -3.0/2.0);
        x32=rt2 - x8;
        x33=mu*rn*x25*x32;
        x34=x33*x7;

        /* Assignment result[2, 2]=rhot1*x34 */
        result[8] = rhot1*x34;
    }

    /* Assignment result[0, 3]=0 */
    result[9] = 0;

    /* Assignment result[1, 3]=Piecewise((0, x28), (rhot2*x29, x24)) */
    double x29;if (x28)
    {
        DEBUG_PRINT("Case (x28) is True.\n");

        /* Assignment result[1, 3]=0 */
        result[10] = 0;
    }
    else if (x24)
    {
        DEBUG_PRINT("Case (x24) is True.\n");
        x21=rt1 - x5;
        x25=pow(x10, -3.0/2.0);
        x26=mu*rn*x21*x25;
        x29=x26*x9;

        /* Assignment result[1, 3]=rhot2*x29 */
        result[10] = rhot2*x29;
    }

    /* Assignment result[2, 3]=Piecewise((rhot2, x16), (0, x20), (rhot2*x23 + rhot2*x35, x24)) */
    double x35;if (x16)
    {
        DEBUG_PRINT("Case (x16) is True.\n");

        /* Assignment result[2, 3]=rhot2 */
        result[11] = rhot2;
    }
    else if (x20)
    {
        DEBUG_PRINT("Case (x20) is True.\n");

        /* Assignment result[2, 3]=0 */
        result[11] = 0;
    }
    else if (x24)
    {
        DEBUG_PRINT("Case (x24) is True.\n");
        x22=1.0/(assert(IS_NOT_ZERO(x11)), x11);
        x23=x14*x22;
        x25=pow(x10, -3.0/2.0);
        x32=rt2 - x8;
        x33=mu*rn*x25*x32;
        x35=x33*x9;

        /* Assignment result[2, 3]=rhot2*x23 + rhot2*x35 */
        result[11] = rhot2*x23 + rhot2*x35;
    }

    /* Assignment result[0, 4]=Piecewise((1, x3), (0, x4)) */
    if (x3)
    {
        DEBUG_PRINT("Case (x3) is True.\n");

        /* Assignment result[0, 4]=1 */
        result[12] = 1;
    }
    else if (x4)
    {
        DEBUG_PRINT("Case (x4) is True.\n");

        /* Assignment result[0, 4]=0 */
        result[12] = 0;
    }

    /* Assignment result[1, 4]=Piecewise((0, x28), (-x21*x30, x24)) */
    double x30;if (x28)
    {
        DEBUG_PRINT("Case (x28) is True.\n");

        /* Assignment result[1, 4]=0 */
        result[13] = 0;
    }
    else if (x24)
    {
        DEBUG_PRINT("Case (x24) is True.\n");
        x21=rt1 - x5;
        x22=1.0/(assert(IS_NOT_ZERO(x11)), x11);
        x30=mu*x22;

        /* Assignment result[1, 4]=-x21*x30 */
        result[13] = -x21*x30;
    }

    /* Assignment result[2, 4]=Piecewise((0, x28), (-x30*x32, x24)) */
    if (x28)
    {
        DEBUG_PRINT("Case (x28) is True.\n");

        /* Assignment result[2, 4]=0 */
        result[14] = 0;
    }
    else if (x24)
    {
        DEBUG_PRINT("Case (x24) is True.\n");
        x22=1.0/(assert(IS_NOT_ZERO(x11)), x11);
        x30=mu*x22;
        x32=rt2 - x8;

        /* Assignment result[2, 4]=-x30*x32 */
        result[14] = -x30*x32;
    }

    /* Assignment result[0, 5]=0 */
    result[15] = 0;

    /* Assignment result[1, 5]=Piecewise((0, x16), (1, x20), (-x27 + x31, x24)) */
    double x31;if (x16)
    {
        DEBUG_PRINT("Case (x16) is True.\n");

        /* Assignment result[1, 5]=0 */
        result[16] = 0;
    }
    else if (x20)
    {
        DEBUG_PRINT("Case (x20) is True.\n");

        /* Assignment result[1, 5]=1 */
        result[16] = 1;
    }
    else if (x24)
    {
        DEBUG_PRINT("Case (x24) is True.\n");
        x21=rt1 - x5;
        x22=1.0/(assert(IS_NOT_ZERO(x11)), x11);
        x23=x14*x22;
        x25=pow(x10, -3.0/2.0);
        x26=mu*rn*x21*x25;
        x27=x26*x7;
        x31=-x23 + 1;

        /* Assignment result[1, 5]=-x27 + x31 */
        result[16] = -x27 + x31;
    }

    /* Assignment result[2, 5]=Piecewise((0, x28), (-x34, x24)) */
    if (x28)
    {
        DEBUG_PRINT("Case (x28) is True.\n");

        /* Assignment result[2, 5]=0 */
        result[17] = 0;
    }
    else if (x24)
    {
        DEBUG_PRINT("Case (x24) is True.\n");
        x25=pow(x10, -3.0/2.0);
        x32=rt2 - x8;
        x33=mu*rn*x25*x32;
        x34=x33*x7;

        /* Assignment result[2, 5]=-x34 */
        result[17] = -x34;
    }

    /* Assignment result[0, 6]=0 */
    result[18] = 0;

    /* Assignment result[1, 6]=Piecewise((0, x28), (-x29, x24)) */
    if (x28)
    {
        DEBUG_PRINT("Case (x28) is True.\n");

        /* Assignment result[1, 6]=0 */
        result[19] = 0;
    }
    else if (x24)
    {
        DEBUG_PRINT("Case (x24) is True.\n");
        x21=rt1 - x5;
        x25=pow(x10, -3.0/2.0);
        x26=mu*rn*x21*x25;
        x29=x26*x9;

        /* Assignment result[1, 6]=-x29 */
        result[19] = -x29;
    }

    /* Assignment result[2, 6]=Piecewise((0, x16), (1, x20), (x31 - x35, x24)) */
    if (x16)
    {
        DEBUG_PRINT("Case (x16) is True.\n");

        /* Assignment result[2, 6]=0 */
        result[20] = 0;
    }
    else if (x20)
    {
        DEBUG_PRINT("Case (x20) is True.\n");

        /* Assignment result[2, 6]=1 */
        result[20] = 1;
    }
    else if (x24)
    {
        DEBUG_PRINT("Case (x24) is True.\n");
        x22=1.0/(assert(IS_NOT_ZERO(x11)), x11);
        x23=x14*x22;
        x25=pow(x10, -3.0/2.0);
        x31=-x23 + 1;
        x32=rt2 - x8;
        x33=mu*rn*x25*x32;
        x35=x33*x9;

        /* Assignment result[2, 6]=x31 - x35 */
        result[20] = x31 - x35;
    }
}
Пример #10
0
int improve_network_division(mod_matrix *mod_mat, two_division *division) {
	int *unmoved;
	int *indices;
	int i, j, i_tag, j_tag;
	elem delta_Q;
	elem *score, *improve;
	if ((score = calloc(mod_mat->A_g->n, sizeof(elem))) == NULL) {
		MEMORY_ALLOCATION_FAILURE_AT("improve_network_division: score");
		return 0;
	}
	if ((improve = calloc(mod_mat->A_g->n, sizeof(elem))) == NULL) {
		MEMORY_ALLOCATION_FAILURE_AT("improve_network_division: improve");
		free(score);
		return 0;
	}
	if ((indices = calloc(mod_mat->A_g->n, sizeof(int))) == NULL) {
		MEMORY_ALLOCATION_FAILURE_AT("improve_network_division: indices");
		free(score);
		free(improve);
		return 0;
	}
	if ((unmoved = calloc(mod_mat->A_g->n, sizeof(int))) == NULL) {
		MEMORY_ALLOCATION_FAILURE_AT("improve_network_division: unmoved");
		free(indices);
		free(score);
		free(improve);
		return 0;
	}
	do {
		/* Initialize values in unmoved */
		for(i=0; i<mod_mat->A_g->n; i++) {
			unmoved[i] = 1;
		}

		/* initialize score values */
		if (init_division_score(mod_mat, division->s_vector, score) == 0) {
			free(unmoved);
			free(indices);
			free(score);
			free(improve);
			return 0;
		}



		/* trying to find an improvement of the partition defined by s */
		for(i=0; i<mod_mat->A_g->n; i++) {

			/* find j (ptr_to_j_tag->value) with maximal score[j] */
			for(j_tag = -1, j=0; j<mod_mat->A_g->n; j++) {
				if (unmoved[j]) {
					if (j_tag==-1 || score[j_tag] < score[j]) {
						j_tag = j;
					}
				}
			}


			division->s_vector->values[j_tag] *= -1; /* S[j_tag] = -S[j_tag] */
			indices[i] = j_tag;
			improve[i] = score[j_tag];
			if (i)
				improve[i] += improve[i-1];
			unmoved[j_tag] = 0;

			/* update scores */
			update_division_score(mod_mat, division->s_vector, j_tag, score);
		}
		/* Find the maximum improvement of s and update s accordingly */
		for(i=0, i_tag=0; i<mod_mat->A_g->n; i++)
			if(improve[i_tag]<improve[i])
				i_tag = i;

		for(i=mod_mat->A_g->n-1; i>i_tag; i--)
			division->s_vector->values[indices[i]] *= -1;

		delta_Q = (i_tag==mod_mat->A_g->n-1) ? 0 : improve[i_tag];
		division->quality += delta_Q;

	} while(IS_POSITIVE(delta_Q));
	free(unmoved);
	free(indices);
	free(improve);
	free(score);
	return(1);
}