El objetivo principal de este artículo es analizar la eficiencia de los solucionadores generales para problemas parabólicos con operadores elípticos de potencia fraccionaria. Tales esquemas discretos pueden ser utilizados en casos de operadores elípticos no constantes, mallas espaciales no uniformes y dominios espaciales generales. Se demuestran resultados de estabilidad para todos los algoritmos y la precisión de las aproximaciones obtenidas se estima resolviendo problemas de prueba bien conocidos. Se presenta una modificación del esquema de división de segundo orden, que combina el método de división para resolver localmente el subproblema no lineal y el algoritmo AAA para resolver el subproblema de difusión no local. Se presentan y analizan los resultados de experimentos computacionales.