Hi, I am trying to make a Kranc file to build a thorn to make a Brill wave using CT_Cosmology. I am bogged down with errors in the script and would appreciate some help in getting the thing to work. I am attaching the script and the output err file. Basically what I am trying to do is to build a single kranc file to handle the specification of the initial metric and extrinsic curvature and have the Hamiltonian constraint solved using the CT_Cosmology thorn CT_Multilevel. I have been trying to copy existing thorns but have not been able to iron out the problems.
Thanks for the help! Comer On Wed, Aug 12, 2015 at 11:17 AM, Comer Duncan <[email protected]> wrote: > Here are the two files. Sorry. > > Comer > > On Wed, Aug 12, 2015 at 11:16 AM, Comer Duncan <[email protected]> > wrote: > >> Hi Ian, >> >> 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? >> >> Thanks. >> >> Comer >> >> On Tue, Aug 11, 2015 at 10:45 AM, Ian Hinder <[email protected]> >> wrote: >> >>> >>> On 11 Aug 2015, at 15:07, Comer Duncan <[email protected]> wrote: >>> >>> Hi Ian, >>> >>> 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! >>> >>> >>> Hi Comer, >>> >>> Observations in the order I have them: >>> >>> >>> - 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 >>> >>> >>> epsi -> one, >>> Null, >>> rho1 -> Sqrt[x*x + y*y], >>> >>> in your calculation. Just move this to outside the calculation. >>> >>> 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. >>> >>> 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. >>> >>> 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. >>> >>> 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 >>> >>> You might want to take a look at >>> >>> >>> https://bitbucket.org/einsteintoolkit/einsteinanalysis/src/master/WeylScal4/m/WeylScal4.m >>> >>> 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. >>> >>> -- >>> Ian Hinder >>> http://members.aei.mpg.de/ianhin >>> >>> >> >
CT_BrillAnalytic.m
Description: Binary data
CT_BrillAnalytic.err
Description: Binary data
_______________________________________________ Users mailing list [email protected] http://lists.einsteintoolkit.org/mailman/listinfo/users
