The ebcs appear to be given correctly. They should be applied automatically, when using simple.py. In your eigenvalue problem solver you should take care of that by calling problem.time_update(), as in schroedinger_app.py, line 104.

Does it help?

r.

----- Reply message -----
From: "steve" <s...@spvi.com>
To: <sfep...@googlegroups.com>
Subject: Cylindrical coordinates
Date: Fri, Jun 8, 2012 14:54


I've simplified this to a case the works a little better (no crash) but still not there.

I've started with the gmsh tutorial file 1: t1.geo and modified it a bit

------------------------------------------
lc = 0.009;

Point(1) = {0, 0, 0, lc};
Point(2) = {.1, 0,  0, lc} ;
Point(3) = {.1, .3, 0, lc} ;
Point(4) = {0,  .3, 0, lc} ;

Line(1) = {1,2} ;
Line(2) = {3,2} ;
Line(3) = {3,4} ;
Line(4) = {4,1} ;

Line Loop(5) = {4,1,-2,3} ;

Plane Surface(6) = {5} ;

Physical Point(1) = {1,2} ;
Physical Point(2) = {3,4} ;

Physical Surface(7) = {6};
-------------------------------

I'd like to set the voltage along line 1 to 0V, and along line 3 to 150V say. I tried 'Physical Line', but that added '2_2' elements to the vtk file that sfepy didn't like. So.. I tried 'Physical Point' and that seemed to allow my code to run, but its still not setting up the region correctly. I used:

regions = {
    'Omega' : ('all', {}),
    'Omega_Source' : ('nodes of group 1', {}),
    'Omega_Target' : ('nodes of group 2', {}),
}

which ran "OK", but didn't set the voltages with:

ebcs = {
    'u1' : ('Omega_Source', {'u.0' : 0.0}),
    'u2' : ('Omega_Target', {'u.0' : 150.0}),
}

I'll keep looking. Thanks for any insight.

-steve



--
You received this message because you are subscribed to the Google Groups "sfepy-devel" group.
To view this discussion on the web visit https://groups.google.com/d/msg/sfepy-devel/-/XDGjDglwIKwJ.
To post to this group, send email to sfepy...@googlegroups.com.
To unsubscribe from this group, send email to sfepy-devel...@googlegroups.com.
For more options, visit this group at http://groups.google.com/group/sfepy-devel?hl=en.