← Back to branch summary

~cbc-core/cbc.solve/main

« back to all changes in this revision

Viewing changes to cbc/swing/new_dualproblem.py

  • Committer: Gabriel Balaban
  • Date: 2012-08-28 20:52:20 UTC
  • Revision ID: gabrib@math.uio.no-20120828205220-i0q978uktc2l9ylo
"FE" method set to default for Dual solver

Show diffs side-by-side

added added

removed removed

Lines of Context:
cbc/swing/new_dualproblem.py
15
15
                      Z_F,  Y_F,  X_F,  Z_S,  Y_S,  Z_M,  Y_M,
16
16
                      Z_F0, Y_F0, X_F0, Z_S0, Y_S0, Z_M0, Y_M0,
17
17
                      U_F0, P_F0, U_S0, P_S0, U_M0,
18
 
                      U_F1, P_F1, U_S1, P_S1, U_M1,parameters):
 
18
                      U_F1, P_F1, U_S1, P_S1, U_M1, parameters):
19
19
    """
20
20
    Return bilinear and linear forms for a time step
21
21
    method - FE is forward Euler, BE is backward Euler, CG1 is midpoint rule
22
22
    """
23
23
 
24
24
    # Choose method here
25
 
    method = "CG1"
 
25
    method = parameters["dualsolver"]["timestepping"]
 
26
    parameters["FSINewtonSolver"]["fluid_domain_time_discretization"] = \
 
27
    parameters["dualsolver"]["fluid_domain_time_discretization"]                                                                  
26
28
 
27
29
    info_blue("Creating dual forms")
28
30