Ignore:
Timestamp:
Mar 10, 2007, 2:08:45 PM (16 years ago)
Author:
Argimiro Resende Secchi
Message:

Removed warning messages of sample_car.mso.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/newlanguage/sample/miscellaneous/sample_car.mso

    r176 r198  
    4949
    5050        EQUATIONS
    51         diff([xl, yl, xr, yr]) = q;
     51        diff([xl, yl, xr, yr])*'s' = q;
    5252       
    53         eps^2*M/2 * diff(q(1)) = (l0-ll)*xl/ll
     53        eps^2*M/2 * diff(q(1))*'s' = (l0-ll)*xl/ll
    5454                +lambda(1)*xb +2*lambda(1)*(xl-xr);
    5555       
    56         eps^2*M/2 * diff(q(2)) = (l0-ll)*yl/ll
     56        eps^2*M/2 * diff(q(2))*'s' = (l0-ll)*yl/ll
    5757                +lambda(1)*yb +2*lambda(2)*(yr-yl) - eps^2*M/2;
    5858
    5959       
    60         eps^2*M/2 * diff(q(3)) = (l0-lr)*(xr-xb)/lr
     60        eps^2*M/2 * diff(q(3))*'s' = (l0-lr)*(xr-xb)/lr
    6161                -2*lambda(2)*(xl-xr);
    6262       
    63         eps^2*M/2 * diff(q(4)) = (l0-lr)*(yr-yb)/lr
     63        eps^2*M/2 * diff(q(4))*'s' = (l0-lr)*(yr-yb)/lr
    6464                -2*lambda(2)*(yl-yr) - eps^2*M/2;
    6565       
     
    7070        xb = sqrt(l^2-yb^2);
    7171       
    72         yb = r*sin(w*time);
     72        yb = r*sin(w*time*'rad/s');
    7373       
    7474        ll = sqrt(xl^2 + yl^2);
     
    8686        TimeStep = 0.01;
    8787        TimeEnd = 3;
    88         #integration = "index0"; # "original";
    89         #DAESolver = "mebdf";
    90 
     88        Integration = "index0";# "index0"; # "original";
     89        DAESolver = "mebdf";
     90        NLASolver = "sundials"; # "nlasolver";
    9191        RelativeAccuracy = 1e-5;
    9292        AbsoluteAccuracy = 1e-5;
    93         #indVarAccuracy = 1e-3;
    9493end
Note: See TracChangeset for help on using the changeset viewer.