<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"><<a href="mailto:comer.duncan@gmail.com" target="_blank">comer.duncan@gmail.com</a>></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'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 'passed' (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"><<a href="mailto:ian.hinder@aei.mpg.de" target="_blank">ian.hinder@aei.mpg.de</a>></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 <<a href="mailto:comer.duncan@gmail.com" target="_blank">comer.duncan@gmail.com</a>> wrote:</div><br><blockquote type="cite"><div dir="ltr">Hi Ian,<div><br></div><div>Yes, I have a recurring error which I can'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 '[' and a ']' . 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'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'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 -> one,</div><div> Null,</div><div> rho1 -> 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["KrancThorn`]"; should be Get["KrancThorn`"]; *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 ") 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'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>