A common technique used in old analog function generators was a two step process. Convert the square wave into a triangle/sawtooth waveform using an integrator (easy with an op amp), and follow that with a PieceWise Linear wave shaping circuit, often using a multiple diode shaping circuit, where diodes are used to progressively shape the triangle into a near sine wave.
The quality of a piecewise linear sine shaper output is critically dependent on the input amplitude.
Unfortunately, the effective gain of the square to triangle integrator is inversely proportional to frequency. The traditional way to handle that was to gang the adjustment of the integrator components with that of the squarewave oscillator, so the integrator gain tracked with the oscillator frequency, but that's no good for an external signal.
One is left with the problem of designing an AGC circuit that can handle a 1000:1 range of input levels while maintaining near constant output. Even if you do forward correction with a F to V circuit and an analog multiplier that's still tricky, and you'll probably need a PID servo loop to trim the gain.

Some improvement can be gained at the expense of considerable complexity by switching integrator capacitors by frequency decade (or even simply hi/lo), with some hysterisis so it doesn't jitter.
Any deviation from an exact 50% duty cycle will cause the integrator output to 'walk' until it rails and distorts, so you'll also need another PID servo loop to inject a centering bias.
Both PID loops will need tuning to handle the maximum possible rate of change of input frequency.
That's an Eurocard protoboard or two of OPAMPs and precision passives, with a considerable price tag. Its not surprising we are pushing the O.P towards digital methods.