Overview
Comment:The solution vector (_ui) is no longer re-zeroed each step.
Retaining the previous value appears to give dramatically better performance.

Messages to the console are now buffered through clog instead of using cerr.

Downloads: Tarball | ZIP archive
Timelines: family | ancestors | descendants | both | origin/master | trunk
Files: files | file ages | folders
SHA3-256: eb84d8e84bc7c3468edbd09108f190e3867511431202b71e7cd653b3f794f88f
User & Date: geraint@users.sourceforge.net on 2002-07-15 12:50:51.000
Other Links: branch diff | manifest | tags
Context
2002-07-15
17:00:36
Writes simulation time to shell when solution converges. check-in: 14e7fa093e user: geraint@users.sourceforge.net tags: origin/master, trunk
12:50:51
The solution vector (_ui) is no longer re-zeroed each step.
Retaining the previous value appears to give dramatically better performance.

Messages to the console are now buffered through clog instead of using cerr. check-in: eb84d8e84b user: geraint@users.sourceforge.net tags: origin/master, trunk

2002-07-14
19:02:51
Improved handling of new plot windows - multiplot and other options affect the new window instead of always controlling the first. check-in: ae9d5867bb user: geraint@users.sourceforge.net tags: origin/master, trunk
Changes
Unified Diff Ignore Whitespace Patch
Modified mttroot/mtt/lib/cc/mtt_AlgebraicSolver.cc from [ab49594ff5] to [021b3b02c0].
︙︙
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23

{
  if (_nyz > 0)
    {
      _x = x;
      _uui.insert(u,0);
      _t = t;
      _par = par;
      _ui = ColumnVector(_nyz,0.0);
      Solve();
      _uui.insert(_ui,_nu);
    }
  else
    {
      _uui = u;
    }








<








9
10
11
12
13
14
15

16
17
18
19
20
21
22

{
  if (_nyz > 0)
    {
      _x = x;
      _uui.insert(u,0);
      _t = t;
      _par = par;

      Solve();
      _uui.insert(_ui,_nu);
    }
  else
    {
      _uui = u;
    }
︙︙
Modified mttroot/mtt/lib/cc/mtt_Hybrd_Solver.cc from [478c947ac8] to [9ffbf825ef].
︙︙
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66

      break;
    default:
      unknown_errors++;
      break;
    }
  if (1 != info)
    {
      std::cerr
	<< " time " << _t << " \t"
	<< " converge (" << convergences << ") "
	<< " limit (" << limit_errors << ")"
	<< " (max error = " << std::abs (eval(_ui).max()) << ")"
	<< " other (" << input_errors + user_errors + progress_errors + unknown_errors << ") "
	<< std::endl;
    }








|








52
53
54
55
56
57
58
59
60
61
62
63
64
65
66

      break;
    default:
      unknown_errors++;
      break;
    }
  if (1 != info)
    {
      std::clog
	<< " time " << _t << " \t"
	<< " converge (" << convergences << ") "
	<< " limit (" << limit_errors << ")"
	<< " (max error = " << std::abs (eval(_ui).max()) << ")"
	<< " other (" << input_errors + user_errors + progress_errors + unknown_errors << ") "
	<< std::endl;
    }
︙︙
Modified mttroot/mtt/lib/cc/mtt_Solver.cc from [effe0746c4] to [bc6c28056d].
︙︙
8
9
10
11
12
13
14

15
16

		     const int nyz)
{
  _np  = npar; 
  _nu  = nu;
  _nx  = nx;
  _ny  = ny;
  _nyz = nyz;

  _uui = ColumnVector (_nu+_nyz);
};








>



8
9
10
11
12
13
14
15
16
17

		     const int nyz)
{
  _np  = npar; 
  _nu  = nu;
  _nx  = nx;
  _ny  = ny;
  _nyz = nyz;
  _ui  = ColumnVector (_nyz,0.0);
  _uui = ColumnVector (_nu+_nyz);
};
MTT: Model Transformation Tools
[ GitHub | SourceHut | Sourceforge | Fossil RSS ]