<div dir="ltr">Here are the two files. Sorry.<div><br></div><div>Comer</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Wed, Aug 12, 2015 at 11:16 AM, Comer Duncan <span dir="ltr">&lt;<a href="mailto:comer.duncan@gmail.com" target="_blank">comer.duncan@gmail.com</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Ian,<div><br></div><div>I&#39;ve made some changes taking info from the WeylScalar thorn stuff, but apparently not really well enough.  I am attaching the BrillAnalytic.err and BrillAnalytic.m files.  The error I get indicates that there is a big mismatch between expected argument number (2) and &#39;passed&#39; (17)!.  I can not tell where this is coming from from the error.  Can you take a look, please?</div><div><br></div><div>Thanks.</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Comer</div></font></span></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Aug 11, 2015 at 10:45 AM, Ian Hinder <span dir="ltr">&lt;<a href="mailto:ian.hinder@aei.mpg.de" target="_blank">ian.hinder@aei.mpg.de</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><br><div><span><div>On 11 Aug 2015, at 15:07, Comer Duncan &lt;<a href="mailto:comer.duncan@gmail.com" target="_blank">comer.duncan@gmail.com</a>&gt; wrote:</div><br><blockquote type="cite"><div dir="ltr">Hi Ian,<div><br></div><div>Yes, I have a recurring error which I can&#39;t track down. In fact there are probably several things wrong with the script.  I am attaching the script BrillAnalytic.m and BrillAnalytic.err.  I can not seem to locate a missing paren or bracket or brace.  vi tells me that there is a mismatch of 1 in the occurrence of a &#39;[&#39; and a &#39;]&#39; .  Also I am not at all sure that I have configured the calculation of TrRicci correctly.  On one hand I specify gxx,gxy,gxz, etc but never associate these expressions with g[la,lb].  I guess I don&#39;t understand well enough the working of Kranc.  Thanks for taking a look!</div></div></blockquote><div><br></div></span><div>Hi Comer,</div><div><br></div><div>Observations in the order I have them:</div><div><br></div><div><ul><li>You can&#39;t define q[rho1_,zcyl_] in the way that you have, inside the calculation.  The mathematica expression f[x_] := y has the effect of defining the function f, but returns Null.  So you effectively have </li></ul><div><br></div><div><div>    epsi -&gt; one,</div><div>    Null,</div><div>    rho1 -&gt; Sqrt[x*x + y*y],</div></div><div><br></div><div>in your calculation.  Just move this to outside the calculation.</div><div><br></div><div>Get[&quot;KrancThorn`]&quot;; should be Get[&quot;KrancThorn`&quot;];  *that* was a subtle mistake!  I commented out the whole of the rest of the file before I found it! The issue is that the closing parenthesis is inside the string (before the &quot;) so Mathematica thinks the whole of the rest of the file is inside the Get[], and when it reaches the end of the file, it complains that the Get hasn&#39;t been terminated.  This is your main error.</div><div><br></div><div>1/detg detgExpr is indeed 1.  However, detg is calculated from g only once, rather than in every component of gu, so it is more efficient.</div><div><br></div><div>You need to define g as a tensor.  Anything which is used as a tensor in the .m file needs to be defined as a tensor with DefineTensor.  This is the cause of the RecursionLimit errors, which happen in MatrixInverse.  Kranc could in fact detect this; I am thinking about how to add such a check.</div><div><br></div><div>The standard 3+1 metric variables are defined in the ADMBase thorn, and are called gxx, gxy, etc.  Kranc needs tensors to be called g11, g12, etc, so we usually just redefine these components: g11=gxx; g21=gxy; g22=gyy; g31=gxz; g32=gyz; g33=gzz</div></div><div><br></div><div>You might want to take a look at </div><div><br></div><div><a href="https://bitbucket.org/einsteintoolkit/einsteinanalysis/src/master/WeylScal4/m/WeylScal4.m" target="_blank">https://bitbucket.org/einsteintoolkit/einsteinanalysis/src/master/WeylScal4/m/WeylScal4.m</a></div><div><br></div><div>which is a thorn for computing the Weyl scalars from the ADMBase variables, which is similar to what you are trying to do in TrRicciCalc.</div></div><span><br><div>
<div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div style="color:rgb(0,0,0);letter-spacing:normal;text-align:start;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;word-wrap:break-word"><div>-- </div><div>Ian Hinder</div><div><a href="http://members.aei.mpg.de/ianhin" target="_blank">http://members.aei.mpg.de/ianhin</a></div></div></div></div></div>
</div>
<br></span></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>