En este artículo construimos solucionadores paralelos y analizamos la eficiencia y precisión de solucionadores paralelos generales para problemas parabólicos tridimensionales con la potencia fraccionaria de operadores elípticos. El método discreto propuesto está dirigido a operadores elípticos generales no constantes, la segunda motivación para el uso de tales esquemas surge cuando son esenciales mallas espaciales no uniformes. Se requieren solucionadores paralelos para resolver los grandes sistemas obtenidos de ecuaciones lineales. Se realiza un análisis detallado de la escalabilidad para comparar la eficiencia de los algoritmos paralelos propuestos. Se presentan y analizan los resultados de experimentos computacionales.