OK, let me elaborate. As you said, SPICE (or anything based on numerical integration) has issues with functions which are discontinuous or have undefined derivatives at certain points. To that effect, convergence and simulation speed can be improved by smoothing out the transfer function. One could theoretically achieve any degree of simulation speed improvement by "ironing out" the transfer function with a progressively heavier hand. However, it appears to me that the transfer function of an NE-2 has an actual discontinuity in it (or at least something that looks an awful lot like a discontinuity), and therefore that if you speed up simulation too much, you throw the essence of the device out with the bath water. So my question was which side are you aiming for? Are you trying to make a highly accurate model that runs as fast as it can under the circumstances, or the fastest model which is as accurate as it can be under the circumstances?