ScharfetterGummel

An old sample code that solves semiconductor drift-diffusion equations.

You should be able to run it with the arguments:

This is only here as a guide of how to use a `ScharfetterGummelFaceVariable` and how to do Newton iterations. This script runs with FiPy 2.1, but it diverges with PySparse and the current is not stable with Trilinos. Caveat Programmator.

Note: the ScharfetterGummelFaceVariable was included in FiPy after this code was written.