OK! Good luck with this... see output for a few modes. Thanks for all the help!<div><br></div><div>Now.. I have Yet Another Problem... gah. Back to the electrostatic problem with Dirichlet BCs.&nbsp;</div><div><br></div><div>This is what I'd love:</div><div><br></div><div>Build a .geo file with various points, lines and surfaces.&nbsp;</div><div><br></div><div>In the geo file specify that certain lines are tagged with a specific ID.</div><div><br></div><div>In the problem description set up boundary conditions by specifying the value of the potential on all the nodes that are located on those lines.</div><div><br></div><div>Is there any way to do that? Everything I've tried so far as failed. I'm trying something like:</div><div><br></div><div>...</div><div>// many lines skipped...</div><div><br></div><div><div>Line(17) = {9, 8};</div><div>Line(18) = {10, 11};</div><div>Line(19) = {11, 12};</div><div>Line(20) = {12, 9};</div><div>Line(21) = {12, 8};</div><div>Line(22) = {12, 13};</div><div><br></div><div>Line Loop(1) = {1, 2, 3, 4};</div><div>Line Loop(2) = {5,6,7,-3};</div><div>Line Loop(3) = {14,-6,-9};</div><div>Line Loop(5) = {11,12,13,-8};</div><div>Line Loop(6) = {-13, -21, 22, -14};</div><div>Line Loop(8) = {-17, -20, 21};</div><div>Line Loop(9) = {18, 19, 20, -16};</div><div><br></div><div>Plane Surface(1) = {1};</div><div>Plane Surface(2) = {2};</div><div>Plane Surface(3) = {3};</div><div>Plane Surface(5) = {5};</div><div>Plane Surface(6) = {6};</div><div>Plane Surface(8) = {8};</div><div>Plane Surface(9) = {9};</div><div><br></div><div>Physical Surface(7) = {1,2,3,5,6,8,9};</div><div><br></div><div>Physical Line(1) = {7,4};</div><div>Physical Line(2) = {5,9,8};</div><div>Physical Line(3) = {12,17,16};</div><div>Physical Line(4) = {19};</div><div>Physical Line(5) = {1};</div><div><br></div><div>I'm thinking that by making "Physical Lines" I can refer to them in the problem description by id. Is that right?</div><div><br></div><div>Anyway.. when I include these in the .geo file I get:</div><div><br></div><div><div>dfr-slusher:sfepy steve$ python simple.py es-test.py</div><div>sfepy: left over: ['verbose', '__builtins__', 'get_r_dependence', '__doc__', '__name__', 'Newregions', 'data_dir', 'nm', '__package__', '_filename', '__file__']</div><div>sfepy: reading mesh (/Users/steve/Development/sfepy/tmp/es-lens.vtk)...</div><div>sfepy: ...done in 0.04 s</div><div>Traceback (most recent call last):</div><div>&nbsp; File "simple.py", line 129, in &lt;module&gt;</div><div>&nbsp; &nbsp; main()</div><div>&nbsp; File "simple.py", line 122, in main</div><div>&nbsp; &nbsp; app = PDESolverApp(conf, options, output_prefix)</div><div>&nbsp; File "/Users/steve/Development/sfepy/sfepy/applications/pde_solver_app.py", line 304, in __init__</div><div>&nbsp; &nbsp; **kwargs)</div><div>&nbsp; File "/Users/steve/Development/sfepy/sfepy/fem/problemDef.py", line 82, in from_conf</div><div>&nbsp; &nbsp; domain = Domain(mesh.name, mesh)</div><div>&nbsp; File "/Users/steve/Development/sfepy/sfepy/fem/domain.py", line 194, in __init__</div><div>&nbsp; &nbsp; gel = GeometryElement(desc)</div><div>&nbsp; File "/Users/steve/Development/sfepy/sfepy/fem/geometry_element.py", line 132, in __init__</div><div>&nbsp; &nbsp; gd = geometry_data[name]</div><div>KeyError: '2_2'</div></div><div><br></div><div>I tried a post-mortem debug.. and indeed 'desc' is '2_2', but there is no matching key in geometry_data. ;-(</div><div><br></div><div>If this is documented, I missed it somehow. Any hints?</div><div><br></div><div>thanks!</div><div>-steve</div><div><br></div><br>On Thursday, June 7, 2012 8:01:52 AM UTC-6, steve wrote:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;">Hi Robert,<div><br><br>On Thursday, June 7, 2012 4:39:02 AM UTC-6, Robert Cimrman wrote:<blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex">On 06/07/2012 07:09 AM, steve wrote:
<br><br>&gt; The trick was:
<br>&gt;
<br>&gt; &nbsp; &nbsp; &nbsp;ebc_1 = {'name':'dummy', 'region':'Surface', 'dofs':{}}
<br>
<br>or, literally, just:
<br>
<br>ebcs = {}
<br></blockquote><div><br></div><div>Dang.. you're right. I thought *sure* I tried that and failed.&nbsp;</div><blockquote class="gmail_quote" style="margin:0;margin-left:0.8ex;border-left:1px #ccc solid;padding-left:1ex">
<br>- the short syntax form... :)
<br>
<br><br>There might be some spurios effects due to mesh. What happens if you refine?&nbsp;</blockquote><div><br></div><div>Demonstrates my ignorance... not sure what that means.... I'll look it up.</div><div><br></div><div>-steve</div><div>&nbsp;</div></div></blockquote></div>