<div dir="ltr"><div><div><div><div><div><div>Hello,<br></div> I need to assemble the operators for an eigenvalue problem. <br></div>Consider for example a 1D Schrodinger equation<br><br></div>-div(grad(psi)) = E psi for x in [0, L] and psi(0) = 0 and psi(L) = 0<br>
<br></div>I understand how to get the 1D mesh into a DM object. <br><br></div><div>I was thinking to create the default section using 1 scalar field per node add boundary points.<br></div><div><br></div><div>Use DMCreateMatrix to get the operators.<br>
<br></div><div>Now if i do element by element assembly, <br>Do i need to use matsetclosure or masetvalues?<br>how are the elements with boundary nodes handled ? <br></div><div>Are the values for boundary nodes ignored ?<br>
<br></div><div>I will pass the assembled operators to slepc to solve the eigenvalue problem. <br><br></div><div>Thanks<br>Reddy<br></div><div><br></div></div><div><div><div><div><div><div><div><div><div><br><br>-- <br>-----------------------------------------------------<br>
Dharmendar Reddy Palle<br>Graduate Student<br>Microelectronics Research center,<br>University of Texas at Austin,<br>10100 Burnet Road, Bldg. 160<br>MER 2.608F, TX 78758-4445<br>e-mail: <a href="mailto:dharmareddy84@gmail.com" target="_blank">dharmareddy84@gmail.com</a><br>
Phone: +1-512-350-9082<br>United States of America.<br>Homepage: <a href="https://webspace.utexas.edu/~dpr342" target="_blank">https://webspace.utexas.edu/~dpr342</a><br>
</div></div></div></div></div></div></div></div></div></div>