<div dir="ltr"><div><br></div><div>The problem is that the node families come in two flavors: with endpoints and without endpoints.  When you left the node family as `DEFAULT`, it knew to choose "with endpoints" for continuous nodes, but by overriding the node family it didn't.  I made an MR and fixed that logic.  You can fix it before the MR gets merged from the command line with "-potential_petscdualspace_lagrange_node_endpoints 1"</div></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">On Mon, Jul 7, 2025 at 3:04 PM Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div>Try running this:</div><div><br></div><div>make -f ./gmakefile test search="snes_tutorials-ex13_2d_p1_conv" EXTRA_OPTIONS="-potential_petscspace_degree 3 -potential_petscdualspace_lagrange_node_type equispaced -petscds_view ::ascii_info_detail"</div><div><br></div><div>You see that</div><div><br></div><div># >           Dual basis vector 0<br># >               Quadrature on a unknown of order -1 on 1 points (dim 2)<br># >               p0 (-0.333333, -0.333333) +1<br># >           Dual basis vector 1<br># >               Quadrature on a unknown of order -1 on 1 points (dim 2)<br># >               p0 (-0.25, -1) +1<br># >           Dual basis vector 2<br># >               Quadrature on a unknown of order -1 on 1 points (dim 2)<br># >               p0 (+0.25, -1) +1<br># >           Dual basis vector 3<br># >               Quadrature on a unknown of order -1 on 1 points (dim 2)<br># >               p0 (+0.25, -0.25) +1<br># >           Dual basis vector 4<br># >               Quadrature on a unknown of order -1 on 1 points (dim 2)<br># >               p0 (-0.25, +0.25) +1<br># >           Dual basis vector 5<br># >               Quadrature on a unknown of order -1 on 1 points (dim 2)<br># >               p0 (-1, +0.25) +1<br># >           Dual basis vector 6<br># >               Quadrature on a unknown of order -1 on 1 points (dim 2)<br># >               p0 (-1, -0.25) +1<br># >           Dual basis vector 7<br># >               Quadrature on a unknown of order -1 on 1 points (dim 2)<br># >               p0 (-1, -1) +1<br># >           Dual basis vector 8<br># >               Quadrature on a unknown of order -1 on 1 points (dim 2)<br># >               p0 (+1, -1) +1<br># >           Dual basis vector 9<br># >               Quadrature on a unknown of order -1 on 1 points (dim 2)<br># >               p0 (-1, +1) +1<br></div><div><br></div><div>but I think those should be +- 1/3, not 1/4. What is going on?</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><span class="gmail_signature_prefix">-- </span><br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="https://urldefense.us/v3/__http://www.cse.buffalo.edu/*knepley/__;fg!!G_uCfscf7eWS!cmxlEUvaAY8SqXfiOuSkCpfTppdZ87z5AnqoBNY8TFup0nd5r-RL44RvjByulFlB3Fm9amh6iNF8hZ6MBYv_M9MQ$" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>
</blockquote></div>