The Newman-Penrose curvature scalars are components of the trace-free Ricci tensor and Weyl tensor with respect to a null tetrad [L, N, M, barM] for a Lorentz metric g on a 4-dimensional manifold.

If $S$ is the trace-free Ricci tensor and $R$ is the trace of the Ricci tensor, then the Newman-Penrose Ricci scalars are:

If $W$ is the Weyl tensor, then the Weyl scalars are:

The minus sign appears in these formulas because there is a sign difference in the definition of the curvature between what is used in the Newman-Penrose formalism and what is used in the DifferentialGeometry CurvatureTensor command.

The command NPCurvatureScalars returns a pair of tables NPRicci, NPWeyl. The table NPRicci contains the values of the NP Ricci scalars and has indices "Phi00", "Phi01", "Phi02", "Phi11","Phi12", "Phi22". The table NPWeyl contains the values of the NP Weyl scalars and has indices "Psi0", "Psi1", "Psi2", "Psi3", "Psi4".

With output = ["RicciScalars"] only the table of Ricci scalars is computed.  With output = ["AllRicciScalars"] the values of the conjugate Ricci scalars are included. With output = ["WeylScalars"] only the table of Weyl scalars is computed. With output = ["AllWeylScalars"] the values of the conjugate Weyl scalars are included. When the output list includes "AllRicciScalars" or  "AllWeylScalars", the tetrad has to be passed as a list of vectors or as the name of an initialized frame created from a null tetrad.

If the curvature tensor $R$ of the metric $g$ has been previously calculated, then passing $R$ to the command NPCurvatureScalars will avoid the need to recalculate it.  For this calling sequence, any rank 4 tensor $R$ with indices ["cov_bas", "cov_bas", "cov_bas", "cov_bas"] can be used.

If the Ricci tensor $S$ of the metric $g$ has been previously calculated and only the Ricci scalars are to be computed, then passing the value of $S$ to the command NPCurvatureScalars will avoid the need to recalculate it. For this calling sequence, any symmetric, covariant rank 2 tensor $S$ can be used.

With either of the first two calling sequences, the Ricci and Weyl scalars are computed directly from the above formulas.

With the third calling sequence the Ricci and Weyl scalars are computed from the NP Ricci identities using the Ricci and Weyl scalars. These identities require knowledge of the null tetrad.  If the current frame is a coordinate frame, then the null tetrad must be specified as the second argument.  If the current frame is an anholonomic frame, then that frame is taken to be the null tetrad.

This command is part of the DifferentialGeometry:-Tensor package, and so can be used in the form NPCurvatureScalars(...) only after executing the commands with(DifferentialGeometry); with(Tensor); in that order.  It can always be used in the long form DifferentialGeometry:-Tensor:-NPCurvatureScalars.

$\mathrm{with}\left(\mathrm{DifferentialGeometry}\right)\:$$\mathrm{with}\left(\mathrm{Tools}\right)\:$$\mathrm{with}\left(\mathrm{Tensor}\right)\:$

$\mathrm{DGsetup}\left(\left[t\,x\,y\,z\right]\,Q\right)$

$\mathrm{frame\; name:\; Q}$

$g≔\mathrm{evalDG}\left(\frac{1}{x^2}\mathrm{dt}\&t\mathrm{dt}-\mathrm{dx}\&t\mathrm{dx}-\mathrm{dy}\&t\mathrm{dy}-\mathrm{dz}\&t\mathrm{dz}\right)$

$g≔\frac{1}{x^2}\mathrm{dt}\mathrm{dt}-\mathrm{dx}\mathrm{dx}-\mathrm{dy}\mathrm{dy}-\mathrm{dz}\mathrm{dz}$

$\mathrm{NTetrad}≔\mathrm{evalDG}\left(\left[\frac{1}{2}{2}^{\frac{1}{2}}x\mathrm{D\_t}+\frac{1}{2}{2}^{\frac{1}{2}}\mathrm{D\_z}\,32^{\frac{1}{2}}x\mathrm{D\_t}+22^{\frac{1}{2}}\mathrm{D\_x}+2^{\frac{1}{2}}\mathrm{D\_y}+22^{\frac{1}{2}}\mathrm{D\_z}\,\left(1+\frac{1}{2}\mathrm{I}\right)2^{\frac{1}{2}}x\mathrm{D\_t}+\frac{1}{2}{2}^{\frac{1}{2}}\mathrm{D\_x}+\frac{1}{2}\mathrm{I}{2}^{\frac{1}{2}}\mathrm{D\_y}+\left(1+\frac{1}{2}\mathrm{I}\right)2^{\frac{1}{2}}\mathrm{D\_z}\,\left(1-\frac{1}{2}\mathrm{I}\right)2^{\frac{1}{2}}x\mathrm{D\_t}+\frac{1}{2}{2}^{\frac{1}{2}}\mathrm{D\_x}-\frac{1}{2}\mathrm{I}{2}^{\frac{1}{2}}\mathrm{D\_y}+\left(1-\frac{1}{2}\mathrm{I}\right)2^{\frac{1}{2}}\mathrm{D\_z}\right]\right)$

$\mathrm{NTetrad}≔\left[\frac{\sqrt{2}x}{2}\mathrm{D\_t}+\frac{\sqrt{2}}{2}\mathrm{D\_z}\,3\sqrt{2}x\mathrm{D\_t}+2\sqrt{2}\mathrm{D\_x}+\sqrt{2}\mathrm{D\_y}+2\sqrt{2}\mathrm{D\_z}\,\left(1+\frac{\mathrm{I}}{2}\right)\sqrt{2}x\mathrm{D\_t}+\frac{\sqrt{2}}{2}\mathrm{D\_x}+\frac{\mathrm{I}}{2}\sqrt{2}\mathrm{D\_y}+\left(1+\frac{\mathrm{I}}{2}\right)\sqrt{2}\mathrm{D\_z}\,\left(1-\frac{\mathrm{I}}{2}\right)\sqrt{2}x\mathrm{D\_t}+\frac{\sqrt{2}}{2}\mathrm{D\_x}-\frac{\mathrm{I}}{2}\sqrt{2}\mathrm{D\_y}+\left(1-\frac{\mathrm{I}}{2}\right)\sqrt{2}\mathrm{D\_z}\right]$

$\mathrm{GRQuery}\left(\mathrm{NTetrad}\,g\,''NullTetrad''\right)$

$\mathrm{true}$

$\mathrm{RicciScalars},\mathrm{WeylScalars}≔\mathrm{NPCurvatureScalars}\left(\mathrm{NTetrad}\right)$

$\mathrm{RicciScalars},\mathrm{WeylScalars}≔table\left(\left[''Lambda''=-\frac{1}{6x^2}\,''Phi12''=\frac{4+3\mathrm{I}}{x^2}\,''Phi22''=\frac{10}{x^2}\,''Phi02''=\frac{1+2\mathrm{I}}{x^2}\,''Phi11''=\frac{5}{2x^2}\,''Phi00''=\frac{1}{2x^2}\,''Phi01''=\frac{1+\frac{\mathrm{I}}{2}}{x^2}\right]\right),table\left(\left[''Psi4''=\frac{\mathrm{-6}-8\mathrm{I}}{x^2}\,''Psi3''=\frac{-5\mathrm{I}}{x^2}\,''Psi1''=\frac{1-\frac{\mathrm{I}}{2}}{x^2}\,''Psi0''=\frac{1}{2x^2}\,''Psi2''=\frac{\frac{4}{3}-2\mathrm{I}}{x^2}\right]\right)$

$\mathrm{RicciScalars}\left[''Phi02''\right]$

$\frac{1+2\mathrm{I}}{x^2}$

$\mathrm{WeylScalars}\left[''Psi3''\right]$

$\frac{-5\mathrm{I}}{x^2}$

$\mathrm{NPCurvatureScalars}\left(\mathrm{NTetrad}\,\mathrm{output}=\left[''AllRicciScalars''\right]\right)$

$table\left(\left[''Phi10''=\frac{1-\frac{\mathrm{I}}{2}}{x^2}\,''Lambda''=-\frac{1}{6x^2}\,''Phi12''=\frac{4+3\mathrm{I}}{x^2}\,''Phi22''=\frac{10}{x^2}\,''Phi02''=\frac{1+2\mathrm{I}}{x^2}\,''Phi11''=\frac{5}{2x^2}\,''Phi21''=\frac{4-3\mathrm{I}}{x^2}\,''Phi00''=\frac{1}{2x^2}\,''Phi01''=\frac{1+\frac{\mathrm{I}}{2}}{x^2}\,''Phi20''=\frac{1-2\mathrm{I}}{x^2}\right]\right)$

$C≔\mathrm{CurvatureTensor}\left(g\right)$

$C≔-\frac{2}{x^2}\mathrm{D\_t}\mathrm{dx}\mathrm{dt}\mathrm{dx}+\frac{2}{x^2}\mathrm{D\_t}\mathrm{dx}\mathrm{dx}\mathrm{dt}-\frac{2}{x^4}\mathrm{D\_x}\mathrm{dt}\mathrm{dt}\mathrm{dx}+\frac{2}{x^4}\mathrm{D\_x}\mathrm{dt}\mathrm{dx}\mathrm{dt}$

$\mathrm{NPCurvatureScalars}\left(\mathrm{NTetrad}\,g\,C\,\mathrm{output}=\left[''WeylScalars''\right]\right)$

$table\left(\left[''Psi4''=\frac{\mathrm{-6}-8\mathrm{I}}{x^2}\,''Psi3''=\frac{-5\mathrm{I}}{x^2}\,''Psi1''=\frac{1-\frac{\mathrm{I}}{2}}{x^2}\,''Psi0''=\frac{1}{2x^2}\,''Psi2''=\frac{\frac{4}{3}-2\mathrm{I}}{x^2}\right]\right)$

$R≔\mathrm{evalDG}\left(a\mathrm{dt}\&s\mathrm{dx}+b\mathrm{dy}\&t\mathrm{dy}\right)$

$R≔\frac{a}{2}\mathrm{dt}\mathrm{dx}+\frac{a}{2}\mathrm{dx}\mathrm{dt}+b\mathrm{dy}\mathrm{dy}$

$\mathrm{NPCurvatureScalars}\left(\mathrm{NTetrad}\,g\,R\,\mathrm{output}=\left[''RicciScalars''\right]\right)$

$table\left(\left[''Lambda''=\frac{b}{24}\,''Phi12''=\frac{\mathrm{I}b}{2}+\frac{7xa}{4}+\frac{\mathrm{I}xa}{2}\,''Phi22''=b+6xa\,''Phi02''=-\frac{b}{4}+\frac{xa}{2}+\frac{\mathrm{I}xa}{4}\,''Phi11''=\frac{b}{8}+\frac{xa}{2}\,''Phi00''=0\,''Phi01''=\frac{xa}{8}\right]\right)$

$C≔\mathrm{CurvatureTensor}\left(g\right)$

$C≔-\frac{2}{x^2}\mathrm{D\_t}\mathrm{dx}\mathrm{dt}\mathrm{dx}+\frac{2}{x^2}\mathrm{D\_t}\mathrm{dx}\mathrm{dx}\mathrm{dt}-\frac{2}{x^4}\mathrm{D\_x}\mathrm{dt}\mathrm{dt}\mathrm{dx}+\frac{2}{x^4}\mathrm{D\_x}\mathrm{dt}\mathrm{dx}\mathrm{dt}$

$S≔\left(-\frac{1}{2}\right)\&mult\mathrm{evalDG}\left(\mathrm{RicciTensor}\left(C\right)-\frac{1}{4}g\mathrm{RicciScalar}\left(g\,C\right)\right)$

$S≔-\frac{1}{2x^4}\mathrm{dt}\mathrm{dt}+\frac{1}{2x^2}\mathrm{dx}\mathrm{dx}-\frac{1}{2x^2}\mathrm{dy}\mathrm{dy}-\frac{1}{2x^2}\mathrm{dz}\mathrm{dz}$

$L,N,M,\mathrm{barM}≔\mathrm{op}\left(\mathrm{NTetrad}\right)\:$

$\mathrm{RicciScalars}\left[''Phi00''\right]=-\mathrm{Hook}\left(\left[L\,L\right]\,S\right)$

$\frac{1}{2x^2}=\frac{1}{2x^2}$

$\mathrm{RicciScalars}\left[''Phi01''\right]=-\mathrm{Hook}\left(\left[L\,M\right]\,S\right)$

$\frac{1+\frac{\mathrm{I}}{2}}{x^2}=\frac{1+\frac{\mathrm{I}}{2}}{x^2}$

$\mathrm{RicciScalars}\left[''Phi01''\right]=-\mathrm{Hook}\left(\left[L\,M\right]\,S\right)$

$\frac{1\+\frac{1}{2}\mathrm{I}}{x^2}\=\frac{1\+\frac{1}{2}\mathrm{I}}{x^2}$

$\mathrm{RicciScalars}\left[''Phi11''\right]=-\mathrm{Hook}\left(\left[L\,N\right]\,S\right)$

$\frac{5}{2x^2}=\frac{5}{2x^2}$

$\mathrm{RicciScalars}\left[''Phi12''\right]=-\mathrm{Hook}\left(\left[N\,M\right]\,S\right)$

$\frac{4+3\mathrm{I}}{x^2}=\frac{4+3\mathrm{I}}{x^2}$

$\mathrm{RicciScalars}\left[''Phi22''\right]=-\mathrm{Hook}\left(\left[N\,N\right]\,S\right)$

$\frac{10}{x^2}=\frac{10}{x^2}$

$\mathrm{FD}≔\mathrm{FrameData}\left(\mathrm{NTetrad}\,\mathrm{NP}\right)$

$\mathrm{FD}≔\left[\left[\mathrm{E1}\,\mathrm{E2}\right]=-\frac{6\sqrt{2}\mathrm{E1}}{x}-\frac{\sqrt{2}\mathrm{E2}}{x}+\frac{\left(2-\mathrm{I}\right)\sqrt{2}\mathrm{E3}}{x}+\frac{\left(2+\mathrm{I}\right)\sqrt{2}\mathrm{E4}}{x}\,\left[\mathrm{E1}\,\mathrm{E3}\right]=-\frac{3\sqrt{2}\mathrm{E1}}{2x}-\frac{\sqrt{2}\mathrm{E2}}{4x}+\frac{\left(\frac{1}{2}-\frac{\mathrm{I}}{4}\right)\sqrt{2}\mathrm{E3}}{x}+\frac{\left(\frac{1}{2}+\frac{\mathrm{I}}{4}\right)\sqrt{2}\mathrm{E4}}{x}\,\left[\mathrm{E1}\,\mathrm{E4}\right]=-\frac{3\sqrt{2}\mathrm{E1}}{2x}-\frac{\sqrt{2}\mathrm{E2}}{4x}+\frac{\left(\frac{1}{2}-\frac{\mathrm{I}}{4}\right)\sqrt{2}\mathrm{E3}}{x}+\frac{\left(\frac{1}{2}+\frac{\mathrm{I}}{4}\right)\sqrt{2}\mathrm{E4}}{x}\,\left[\mathrm{E2}\,\mathrm{E3}\right]=\frac{\left(3+6\mathrm{I}\right)\sqrt{2}\mathrm{E1}}{x}+\frac{\left(\frac{1}{2}+\mathrm{I}\right)\sqrt{2}\mathrm{E2}}{x}-\frac{\left(2+\frac{3\mathrm{I}}{2}\right)\sqrt{2}\mathrm{E3}}{x}-\frac{5\mathrm{I}\sqrt{2}\mathrm{E4}}{2x}\,\left[\mathrm{E2}\,\mathrm{E4}\right]=\frac{\left(3-6\mathrm{I}\right)\sqrt{2}\mathrm{E1}}{x}+\frac{\left(\frac{1}{2}-\mathrm{I}\right)\sqrt{2}\mathrm{E2}}{x}+\frac{5\mathrm{I}\sqrt{2}\mathrm{E3}}{2x}+\frac{\left(\mathrm{-2}+\frac{3\mathrm{I}}{2}\right)\sqrt{2}\mathrm{E4}}{x}\,\left[\mathrm{E3}\,\mathrm{E4}\right]=-\frac{3\mathrm{I}\sqrt{2}\mathrm{E1}}{x}-\frac{\mathrm{I}\sqrt{2}\mathrm{E2}}{2x}+\frac{\left(\frac{1}{2}+\mathrm{I}\right)\sqrt{2}\mathrm{E3}}{x}+\frac{\left(-\frac{1}{2}+\mathrm{I}\right)\sqrt{2}\mathrm{E4}}{x}\right]$

$\mathrm{DGsetup}\left(\mathrm{FD}\right)$

$\mathrm{frame\; name:\; NP}$

$h≔\mathrm{evalDG}\left(2\mathrm{\Θ1}\&s\mathrm{\Θ2}-2\mathrm{\Θ3}\&s\mathrm{\Θ4}\right)$

$h≔\mathrm{\Θ1}\mathrm{\Θ2}+\mathrm{\Θ2}\mathrm{\Θ1}-\mathrm{\Θ3}\mathrm{\Θ4}-\mathrm{\Θ4}\mathrm{\Θ3}$

$\mathrm{NPCurvatureScalars}\left(\mathrm{NP}\,\mathrm{output}=\left[''WeylScalars''\right]\right)$

$table\left(\left[''Psi4''=\frac{\mathrm{-6}-8\mathrm{I}}{x^2}\,''Psi3''=\frac{-5\mathrm{I}}{x^2}\,''Psi1''=\frac{1-\frac{\mathrm{I}}{2}}{x^2}\,''Psi0''=\frac{1}{2x^2}\,''Psi2''=\frac{\frac{4}{3}-2\mathrm{I}}{x^2}\right]\right)$

$W≔\mathrm{WeylTensor}\left(h\right)\:$

$\mathrm{WeylScalars}\left[''Psi0''\right]=-\mathrm{Hook}\left(\left[\mathrm{E1}\,\mathrm{E3}\,\mathrm{E1}\,\mathrm{E3}\right]\,W\right)$

$\frac{1}{2x^2}=\frac{1}{2x^2}$

$\mathrm{WeylScalars}\left[''Psi0''\right]=-\mathrm{DGinfo}\left(W\,''CoefficientList''\,\left[\left[1\,3\,1\,3\right]\right]\right)\left[1\right]$

$\frac{1}{2x^2}=\frac{1}{2x^2}$

$\mathrm{WeylScalars}\left[''Psi1''\right]=-\mathrm{DGinfo}\left(W\,''CoefficientList''\,\left[\left[1\,2\,1\,3\right]\right]\right)\left[1\right]$

$\frac{1-\frac{\mathrm{I}}{2}}{x^2}=\frac{1-\frac{\mathrm{I}}{2}}{x^2}$

$\mathrm{WeylScalars}\left[''Psi2''\right]=-\mathrm{DGinfo}\left(W\,''CoefficientList''\,\left[\left[1\,3\,4\,2\right]\right]\right)\left[1\right]$

$\frac{\frac{4}{3}-2\mathrm{I}}{x^2}=\frac{\frac{4}{3}-2\mathrm{I}}{x^2}$

$\mathrm{WeylScalars}\left[''Psi3''\right]=-\mathrm{DGinfo}\left(W\,''CoefficientList''\,\left[\left[1\,2\,4\,2\right]\right]\right)\left[1\right]$

$\frac{-5\mathrm{I}}{x^2}=\frac{-5\mathrm{I}}{x^2}$

$\mathrm{WeylScalars}\left[''Psi4''\right]=-\mathrm{DGinfo}\left(W\,''CoefficientList''\,\left[\left[2\,4\,2\,4\right]\right]\right)\left[1\right]$

$\frac{\mathrm{-6}-8\mathrm{I}}{x^2}=\frac{\mathrm{-6}-8\mathrm{I}}{x^2}$

$\mathrm{SpinCoeff}≔\mathrm{NPSpinCoefficients}\left(\mathrm{NTetrad}\right)$

$\mathrm{SpinCoeff}≔table\left(\left[''lambda''=\frac{-\frac{5\mathrm{I}}{2}\sqrt{2}}{x}\,''sigma''=\frac{\left(\frac{1}{2}+\frac{\mathrm{I}}{4}\right)\sqrt{2}}{x}\,''gamma''=\frac{3\sqrt{2}}{x}-\frac{3\mathrm{I}\sqrt{2}}{2x}\,''tau''=\frac{3\sqrt{2}}{2x}\,''kappa''=\frac{\sqrt{2}}{4x}\,''pi''=\frac{\left(\frac{1}{2}-\mathrm{I}\right)\sqrt{2}}{x}\,''rho''=\frac{\left(\frac{1}{2}-\frac{\mathrm{I}}{4}\right)\sqrt{2}}{x}\,''epsilon''=\frac{\sqrt{2}}{2x}-\frac{\mathrm{I}\sqrt{2}}{4x}\,''alpha''=\frac{\left(\frac{3}{4}-\mathrm{I}\right)\sqrt{2}}{x}\,''nu''=\frac{\left(3-6\mathrm{I}\right)\sqrt{2}}{x}\,''beta''=\frac{5\sqrt{2}}{4x}\,''mu''=\frac{\left(2-\frac{3\mathrm{I}}{2}\right)\sqrt{2}}{x}\right]\right)$

$\mathrm{NPCurvatureScalars}\left(\mathrm{SpinCoeff}\,\mathrm{NTetrad}\,\mathrm{output}=\left[''RicciScalars''\right]\right)$

$table\left(\left[''Lambda''=-\frac{1}{6x^2}\,''Phi12''=\frac{4+3\mathrm{I}}{x^2}\,''Phi22''=\frac{10}{x^2}\,''Phi02''=\frac{1+2\mathrm{I}}{x^2}\,''Phi11''=\frac{5}{2x^2}\,''Phi00''=\frac{1}{2x^2}\,''Phi01''=\frac{1+\frac{\mathrm{I}}{2}}{x^2}\right]\right)$

$\mathrm{ChangeFrame}\left(\mathrm{NP}\right)$

$Q$

$\mathrm{NPCurvatureScalars}\left(\mathrm{SpinCoeff}\,\mathrm{output}=\left[''RicciScalars''\right]\right)$

$table\left(\left[''Lambda''=-\frac{1}{6x^2}\,''Phi12''=\frac{4+3\mathrm{I}}{x^2}\,''Phi22''=\frac{10}{x^2}\,''Phi02''=\frac{1+2\mathrm{I}}{x^2}\,''Phi11''=\frac{5}{2x^2}\,''Phi00''=\frac{1}{2x^2}\,''Phi01''=\frac{1+\frac{\mathrm{I}}{2}}{x^2}\right]\right)$

$\mathrm{DGsetup}\left(\left[t\,x\,y\,z\right]\,\left[\mathrm{z1}\,\mathrm{z2}\,\mathrm{w1}\,\mathrm{w2}\right]\,\mathrm{Spin}\right)$

$\mathrm{frame\; name:\; Spin}$

$\mathrm{W0}≔\mathrm{evalDG}\left(\left(\left(\mathrm{dz1}\&t\mathrm{dz1}\right)\&t\mathrm{dz1}\right)\&t\mathrm{dz1}\right)$

$\mathrm{W0}≔\mathrm{dz1}\mathrm{dz1}\mathrm{dz1}\mathrm{dz1}$

$\mathrm{W1}≔\mathrm{evalDG}\left(\left(\left(\mathrm{dz1}\&t\mathrm{dz1}\right)\&t\mathrm{dz1}\right)\&t\mathrm{dz2}+\left(\left(\mathrm{dz1}\&t\mathrm{dz1}\right)\&t\mathrm{dz2}\right)\&t\mathrm{dz1}+\left(\left(\mathrm{dz1}\&t\mathrm{dz2}\right)\&t\mathrm{dz1}\right)\&t\mathrm{dz1}+\left(\left(\mathrm{dz2}\&t\mathrm{dz1}\right)\&t\mathrm{dz1}\right)\&t\mathrm{dz1}\right)$

$\mathrm{W1}≔\mathrm{dz1}\mathrm{dz1}\mathrm{dz1}\mathrm{dz2}+\mathrm{dz1}\mathrm{dz1}\mathrm{dz2}\mathrm{dz1}+\mathrm{dz1}\mathrm{dz2}\mathrm{dz1}\mathrm{dz1}+\mathrm{dz2}\mathrm{dz1}\mathrm{dz1}\mathrm{dz1}$

$\mathrm{W2}≔\mathrm{evalDG}\left(\left(\left(\mathrm{dz1}\&t\mathrm{dz1}\right)\&t\mathrm{dz2}\right)\&t\mathrm{dz2}+\left(\left(\mathrm{dz1}\&t\mathrm{dz2}\right)\&t\mathrm{dz1}\right)\&t\mathrm{dz2}+\left(\left(\mathrm{dz1}\&t\mathrm{dz2}\right)\&t\mathrm{dz2}\right)\&t\mathrm{dz1}+\left(\left(\mathrm{dz2}\&t\mathrm{dz1}\right)\&t\mathrm{dz1}\right)\&t\mathrm{dz2}+\left(\left(\mathrm{dz2}\&t\mathrm{dz1}\right)\&t\mathrm{dz2}\right)\&t\mathrm{dz1}+\left(\left(\mathrm{dz2}\&t\mathrm{dz2}\right)\&t\mathrm{dz1}\right)\&t\mathrm{dz1}\right)$

$\mathrm{W2}≔\mathrm{dz1}\mathrm{dz1}\mathrm{dz2}\mathrm{dz2}+\mathrm{dz1}\mathrm{dz2}\mathrm{dz1}\mathrm{dz2}+\mathrm{dz1}\mathrm{dz2}\mathrm{dz2}\mathrm{dz1}+\mathrm{dz2}\mathrm{dz1}\mathrm{dz1}\mathrm{dz2}+\mathrm{dz2}\mathrm{dz1}\mathrm{dz2}\mathrm{dz1}+\mathrm{dz2}\mathrm{dz2}\mathrm{dz1}\mathrm{dz1}$

$\mathrm{W3}≔\mathrm{evalDG}\left(\left(\left(\mathrm{dz1}\&t\mathrm{dz2}\right)\&t\mathrm{dz2}\right)\&t\mathrm{dz2}+\left(\left(\mathrm{dz2}\&t\mathrm{dz1}\right)\&t\mathrm{dz2}\right)\&t\mathrm{dz2}+\left(\left(\mathrm{dz2}\&t\mathrm{dz2}\right)\&t\mathrm{dz1}\right)\&t\mathrm{dz2}+\left(\left(\mathrm{dz2}\&t\mathrm{dz2}\right)\&t\mathrm{dz2}\right)\&t\mathrm{dz1}\right)$

$\mathrm{W3}≔\mathrm{dz1}\mathrm{dz2}\mathrm{dz2}\mathrm{dz2}+\mathrm{dz2}\mathrm{dz1}\mathrm{dz2}\mathrm{dz2}+\mathrm{dz2}\mathrm{dz2}\mathrm{dz1}\mathrm{dz2}+\mathrm{dz2}\mathrm{dz2}\mathrm{dz2}\mathrm{dz1}$

$\mathrm{W4}≔\mathrm{DifferentialGeometry}:-\mathrm{evalDG}\left(\left(\left(\mathrm{dz2}\&t\mathrm{dz2}\right)\&t\mathrm{dz2}\right)\&t\mathrm{dz2}\right)$

$\mathrm{W4}≔\mathrm{dz2}\mathrm{dz2}\mathrm{dz2}\mathrm{dz2}$

$W≔\mathrm{evalDG}\left(\mathrm{A0}\mathrm{W0}+\mathrm{A1}\mathrm{W1}+\mathrm{A2}\mathrm{W2}+\mathrm{A3}\mathrm{W3}+\mathrm{A4}\mathrm{W4}\right)\:$

$\mathrm{NPCurvatureScalars}\left(W\right)$

$table\left(\left[''Psi4''=\mathrm{A4}\,''Psi3''=\mathrm{A3}\,''Psi1''=\mathrm{A1}\,''Psi0''=\mathrm{A0}\,''Psi2''=\mathrm{A2}\right]\right)$