| 1 | //STEADYSTATE class definition
|
|---|
| 2 | //
|
|---|
| 3 | // Usage:
|
|---|
| 4 | // steadystate=new steadystate();
|
|---|
| 5 |
|
|---|
| 6 | function steadystate (){
|
|---|
| 7 | //methods
|
|---|
| 8 | this.setdefaultparameters = function(){// {{{
|
|---|
| 9 |
|
|---|
| 10 | //maximum of steady state iterations
|
|---|
| 11 | this.maxiter=100;
|
|---|
| 12 |
|
|---|
| 13 | //Relative tolerance for the steadystate convertgence
|
|---|
| 14 | this.reltol=0.01;
|
|---|
| 15 |
|
|---|
| 16 | //default output
|
|---|
| 17 | this.requested_outputs=['default'];
|
|---|
| 18 |
|
|---|
| 19 |
|
|---|
| 20 | }// }}}
|
|---|
| 21 | this.disp= function(){// {{{
|
|---|
| 22 |
|
|---|
| 23 | console.log(sprintf(' steadystate solution parameters:'));
|
|---|
| 24 |
|
|---|
| 25 | fielddisplay(this,'reltol','relative tolerance criterion');
|
|---|
| 26 | fielddisplay(this,'maxiter','maximum number of iterations');
|
|---|
| 27 | fielddisplay(this,'requested_outputs','additional requested outputs');
|
|---|
| 28 |
|
|---|
| 29 | }// }}}
|
|---|
| 30 | this.checkconsistency = function(md,solution,analyses) {// {{{
|
|---|
| 31 |
|
|---|
| 32 | //Early return
|
|---|
| 33 | if (solution!=SteadystateSolutionEnum()) return;
|
|---|
| 34 |
|
|---|
| 35 | if (md.timestepping.time_step!=0){
|
|---|
| 36 | md.checkmessage('for a steadystate computation, timestepping.time_step must be zero.');
|
|---|
| 37 | }
|
|---|
| 38 | checkfield(md,'fieldname','steadystate.requested_outputs','stringrow',1);
|
|---|
| 39 |
|
|---|
| 40 | if (isNaN(md.stressbalance.reltol)){
|
|---|
| 41 | md.checkmessage('for a steadystate computation, stressbalance.reltol (relative convergence criterion) must be defined!');
|
|---|
| 42 | }
|
|---|
| 43 | } // }}}
|
|---|
| 44 | //properties
|
|---|
| 45 | // {{{
|
|---|
| 46 |
|
|---|
| 47 | this.reltol = 0;
|
|---|
| 48 | this.maxiter = 0;
|
|---|
| 49 | this.requested_outputs = [];
|
|---|
| 50 |
|
|---|
| 51 | this.setdefaultparameters();
|
|---|
| 52 | //}}}
|
|---|
| 53 | }
|
|---|