For your first question, yes you can definitely do this! You can use a couple changepoint kernels. Where you have your parametric model defined, y = c x^{c_2}, use that as your mean function (c and c_2 can be unknown parameters that are given priors) and set your covariance function to return all zeros.
RE your second question, maybe check out PyMC-BART?