Herman Sletmoen

Thanks! And thanks @goerz for guiding me :)

Does it have anything to do with #1305 where [this](https://github.com/JuliaSymbolics/Symbolics.jl/pull/1305/files#diff-b3b88cb2464df04f837776d6b612595cf6b7f6c9daee47aaa07f18b29549a48cR200) now returns a LaTeXString instead of a Symbol?

It hopefully works now, but I think we have room for improvement 😄 I suggest to keep the discussion in #1305.

This sounds somewhat familiar to what I am doing. I have cosmological models instead of chemical plants, and different "submodels" for gravity and particles instead of pump types. See e.g....

I would really appreciate a solution to this as well. I think it would be very helpful to extend the documentation page https://docs.sciml.ai/ModelingToolkit/stable/examples/remake/ with e.g. - How to optimize only...

Also cross-posted in https://github.com/lesgourg/class_public/issues/645.

Thanks for the reply! I disabled reionization in both codes just to keep the comparison minimal. If I change CAMB to use ``` reionization = T re_use_optical_depth = T re_optical_depth...

Disabling the Limber approximation makes a difference (I find it a little ambiguous whether CAMB/CLASS uses the Limber approximation for lensing-polarization in addition to lensing-lensing), but does not resolve the...

> I think the correlation coefficient here is tiny, i.e. |phi-phi/sqrt(EE*phiphi)|

[The current implementation](https://github.com/JuliaSymbolics/Symbolics.jl/blob/cb995c17081aea845aaabf78686ef5d10ed96253/src/taylor.jl) (added by me in https://github.com/JuliaSymbolics/Symbolics.jl/pull/1298) is very barebones and far from optimal. I think the best approach for improvements would be make Symbolics composable with a special-purpose...