1# $Header: /var/cvs/mbdyn/mbdyn/mbdyn-1.0/modules/module-mds/mds,v 1.8 2017/01/12 14:54:57 masarati Exp $
2# MBDyn (C) is a multibody analysis code.
3# http://www.mbdyn.org
4#
5# Copyright (C) 1996-2017
6#
7# Pierangelo Masarati	<masarati@aero.polimi.it>
8# Paolo Mantegazza	<mantegazza@aero.polimi.it>
9#
10# Dipartimento di Ingegneria Aerospaziale - Politecnico di Milano
11# via La Masa, 34 - 20156 Milano, Italy
12# http://www.aero.polimi.it
13#
14# Changing this copyright notice is forbidden.
15#
16# This program is free software; you can redistribute it and/or modify
17# it under the terms of the GNU General Public License as published by
18# the Free Software Foundation (version 2 of the License).
19#
20#
21# This program is distributed in the hope that it will be useful,
22# but WITHOUT ANY WARRANTY; without even the implied warranty of
23# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
24# GNU General Public License for more details.
25#
26# You should have received a copy of the GNU General Public License
27# along with this program; if not, write to the Free Software
28# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
29#
30
31begin: data;
32	problem: initial value;
33end: data;
34
35begin: initial value;
36	initial time: 0.;
37	final time: 10.;
38	time step: 1.e-2;
39
40	tolerance: 1.e-6;
41	max iterations: 10;
42
43	method: ms, .6;
44
45	# derivatives coefficient: 1.e-12;
46	derivatives tolerance: 1.e-3;
47	derivatives max iterations: 10;
48
49	# linear solver: umfpack, cc;
50
51	output: iterations;
52	#output: residual;
53	#output: solution;
54	#output: jacobian;
55end: initial value;
56
57begin: control data;
58	loadable elements: 1;
59end: control data;
60
61set: real M = 1.;
62set: real OMEGA = 2*pi;
63set: real XI = .1;
64set: real K = M*OMEGA^2;
65set: real D = 2*XI*M*OMEGA;
66
67begin: elements;
68	module load: "libmodule-mds";
69	user defined: 1, mds,
70		M, D, K,
71			# sine, 0., .2*OMEGA, 1., forever, 0.,
72			sine, 0., 10.*OMEGA, 1., one, 0.,
73		output, yes;
74end: elements;
75
76# vim:ft=mbd
77