157 void mass_spring_system(int nx, int nu, double *Ac, double *Bc)  in mass_spring_system()  argument
168 		Ac[ii] = 0.0;  in mass_spring_system()
170 		Ac[ii+nx*(pp+ii)] = 1.0;  in mass_spring_system()
172 		Ac[pp+ii+nx*(ii)] = -2.0;  in mass_spring_system()
174 		Ac[pp+ii+nx*(1+ii)] = 1.0;  in mass_spring_system()
176 		Ac[1+pp+ii+nx*(ii)] = 1.0;  in mass_spring_system()
192 	double *Ac;  member
215 	ls->Ac = (double *) c_ptr;  in d_create_linear_system()
229 	double *Ac = ls->Ac;  in d_cvt_colmaj_to_linear_system()  local
232 		Ac[ii] = A[ii];  in d_cvt_colmaj_to_linear_system()
246 	double *Ac = ls->Ac;  in d_expl_linear_ode()  local
252 			xdot[ii] += Ac[ii+nx*jj] * x[jj];  in d_expl_linear_ode()
267 	double *Ac = ls->Ac;  in d_expl_linear_vde()  local
275 				xdot[ii+nx*kk] += Ac[ii+nx*jj] * x[jj+nx*kk];  in d_expl_linear_vde()
295 	double *Ac = ls->Ac;  in d_res_impl_linear_ode()  local
301 			res[ii] += Ac[ii+nx*jj] * x[jj];  in d_res_impl_linear_ode()
316 	double *Ac = ls->Ac;  in d_res_impl_linear_vde()  local
324 				res[ii+nx*kk] += Ac[ii+nx*jj] * x[jj+nx*kk];  in d_res_impl_linear_vde()
343 	double *Ac = ls->Ac;  in d_jac_impl_linear_ode()  local
346 			jac[ii+nx*jj] = Ac[ii+nx*jj];  in d_jac_impl_linear_ode()
368 	double *Ac; d_zeros(&Ac, nx, nx);  in main()  local
371 	mass_spring_system(nx, nu, Ac, Bc);  in main()
373 	d_print_mat(nx, nx, Ac, nx);  in main()
383 	d_cvt_colmaj_to_linear_system(Ac, Bc, &ls);  in main()
385 	d_print_mat(nx, nx, ls.Ac, nx);  in main()
413 		A[ii] = Ts*Ac[ii];  in main()
421 	dgesv_3l(nx, nu, Ac, nx, ipiv, B, nx, &info);  in main()
625 	free(Ac);  in main()