def METHOD_NAME(self, comm):
    """
    Setup mesh and pytacs object for problem we will be testing.
    """
    # Overwrite default check values
    if self.dtype == complex:
        self.rtol = 1e-6
        self.atol = 1e-6
        self.dh = 1e-50
    else:
        self.rtol = 2e-1
        self.atol = 1e-3
        self.dh = 1e-6
    # Instantiate FEA Assembler
    struct_options = {}
    fea_assembler = pytacs.pyTACS(bdf_file, comm, options=struct_options)
    # Set up constitutive objects and elements
    fea_assembler.initialize()
    # Read in forces from BDF and create tacs struct problems
    tacs_probs = fea_assembler.createTACSProbsFromBDF()
    # Convert from dict to list
    tacs_probs = list(tacs_probs.values())
    # Set convergence to be tight for test
    for problem in tacs_probs:
        problem.setOption("L2Convergence", 1e-20)
        problem.setOption("L2ConvergenceRel", 1e-20)
    # Add Functions
    for problem in tacs_probs:
        problem.addFunction("mass", functions.StructuralMass)
        problem.addFunction("compliance", functions.Compliance)
        problem.addFunction(
            "x_disp",
            functions.KSDisplacement,
            ksWeight=ksweight,
            direction=[10.0, 0.0, 0.0],
        )
        problem.addFunction(
            "y_disp",
            functions.KSDisplacement,
            ksWeight=ksweight,
            direction=[0.0, 10.0, 0.0],
        )
        problem.addFunction(
            "z_disp",
            functions.KSDisplacement,
            ksWeight=ksweight,
            direction=[0.0, 0.0, 10.0],
        )
        problem.addFunction(
            "-x_disp",
            functions.KSDisplacement,
            ksWeight=ksweight,
            direction=[-10.0, 0.0, 0.0],
        )
        problem.addFunction(
            "-y_disp",
            functions.KSDisplacement,
            ksWeight=ksweight,
            direction=[0.0, -10.0, 0.0],
        )
        problem.addFunction(
            "-z_disp",
            functions.KSDisplacement,
            ksWeight=ksweight,
            direction=[0.0, 0.0, -10.0],
        )
    constr = fea_assembler.createVolumeConstraint("volume")
    constr.addConstraint("con")
    tacs_probs.append(constr)
    return tacs_probs, fea_assembler