Exemple #1
0
 void
jacinc1_(fint *M, fint *N, fint *NO, fint *NZ, fint *JP,
	fint *JI, real *X, real *L, real *U,
	real *Lrhs, real *Urhs, real *Inf, fint *Objdir)
{
	int n;
	cgrad *gr, **grp;
	ASL *asl = cur_ASL;

	mnnzchk_ASL(asl, M, N, *NZ, "jacinc");
	*Inf = Infinity;
	if ((n = n_con)) {
		LUcopy_ASL(n, Lrhs, Urhs, LUrhs);
		grp = Cgrad + n;
		for(; n > 0; --n)
			for(gr = *--grp; gr; gr = gr->next) {
				JI[gr->goff] = n;
				JP[gr->varno] = gr->goff + 1;
				}
		JP[c_vars] = nzjac + 1;
		}
	LUcopy_ASL(c_vars, L, U, LUv);
	memcpy(X, X0, asl->i.n_var0 * sizeof(real));
	for(n = n_obj; --n >= 0; )
		Objdir[n] = objtype[n];
	}
Exemple #2
0
 void
jacinc_(fint *M, fint *N, fint *NZ, register fint *JP,
	register short *JI, real *X, real *L, real *U,
	real *Lrhs, real *Urhs, real *Inf)
{
	int n;
	cgrad *gr, **grp;
	ASL *asl = cur_ASL;

	mnnzchk_ASL(asl, M, N, NZ, "jacinc");
	*Inf = Infinity;
	if ((n = n_con)) {
		LUcopy_ASL(n, Lrhs, Urhs, LUrhs);
		grp = Cgrad + n;
		for(; n > 0; --n)
			for(gr = *--grp; gr; gr = gr->next) {
				JI[gr->goff] = (short)n;
				JP[gr->varno] = gr->goff + 1;
				}
		JP[c_vars] = nzjac + 1;
		}
	LUcopy_ASL(c_vars, L, U, LUv);
	memcpy(X, X0, asl->i.n_var0 * sizeof(real));
	}