Hi,

I'm afraid that I'll need some specific help. I've tried recasting my code to run in interactive mode but am getting the following error:

I can't see how to resolve this issue from the linear_elasticity_interactive.py case.

I've written the respective variables as:

I've also tried writing:

David

I'm afraid that I'll need some specific help. I've tried recasting my code to run in interactive mode but am getting the following error:

`ValueError: wrong arguments shapes for "dw_laplace.1.Omega(coef.val, s, T)" term! (see above)`

I can't see how to resolve this issue from the linear_elasticity_interactive.py case.

I've written the respective variables as:

`coef = Material('coef', val=[2.0])`

T = FieldVariable('T', 'unknown', field)

s = FieldVariable('s', 'test', field, primary_var_name='T')

I've also tried writing:

coef = Material('coef', values={'.val' : 0.0})

as per its2D_interactive.py, but with the same effect.

Thanks,David

On Tuesday, 6 December 2016 15:26:44 UTC+1, David Jessop wrote:Hi,

I have a non-linear form of the Poisson equation where the "diffusion" coefficient depends on the derivatives of the state variable (see image). There are no Dirichlet-type boundary conditions (i.e. no ebcs) so I really need to specify a good initial guess at the solution.

My questions are:

1) How are initial guesses passed to the solver? From the documentation, it doesn't seem like initial guesses can be directly passed to the nls.newton solver as an argument (though they are an argument of the subroutine __call__).

2) In what form should the initial guess vector/array be written? Can a function be passed? Do numerical values have be be specified at the mesh nodes?

Thanks for any help on this matter.

D