On Fri, Feb 15, 2013 at 11:52 AM, Anders Logg <[email protected]> wrote: > On Fri, Feb 15, 2013 at 11:10:59AM +0100, Roland Siegbert wrote: > My suggestion would be to edit the CBC.Twist code to make it do what > you want. If you find a good solution and a nice interface, it might > be merged back in. OK, for anybody interested in the result: I added a create_solver() routine and modified solve() of the CBCProblem class derived CBC.Twist class Hyperelasticity, in a way that doesn't break the other demos and gave me the possibility to do this using release.py written by Harish Narayanan:
<code> # Setup and solve problem problem = Release() problem.create_solver() bc = DirichletBC(problem.solver.vector.sub(0), Constant(0.0), "x[0] == 1.0") problem.solver.bcu.append(bc) problem.parameters["solver_parameters"]["save_solution"] = True print problem problem.solve() interactive() </code> As can be seen, I also added vector to the instance attributes of the CBCSolver derived class in cbc.Twist (CG1...). I'll write a patch for review, after I understand how launchpad works in that way. ;-) Roland _______________________________________________ Mailing list: https://launchpad.net/~fenics Post to : [email protected] Unsubscribe : https://launchpad.net/~fenics More help : https://help.launchpad.net/ListHelp
