I am trying to evaluate the covariant derivative of a tensor, where both the metric and the tensor depends on some general scalar function \beta[r]. I understand that I should use
DefScalarFunction[\[Beta]]
to declare a general scalar function. I tried doing that and it gave me the correct Ricci scalar of the manifold by RicciScalar[cd][] (which depends on \beta[r] and \beta'[r]). However, when I tried to evaluate some contracted covariant derivatives for a tensor (labelled as h1) that also has \beta[r] in it
h1 = CTensor[{{-2 \[Beta][r] (1 - \[Epsilon] \[Beta][r]), 0}, {0, 0}}, {-flat, -flat}];
cd[-b][cd[-a][h1[a, b]]] // Simplify
the final answer has some terms that look like "D_1[r]" and "D_1D_1[r]", which I think stand for partial derivatives. Since index 1 is r, I would have thought that xAct will automatically output "1" and "0" for those quantities.
I have attached a minimal notebook as an example. The scalar function \beta[r] is defined in line 6 and the covariant derivatives are evaluated in line 12.
Thank you in advance for any help!
Vincent