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 = 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(
            "I_xx",
            functions.MomentOfInertia,
            direction1=[1.0, 0.0, 0.0],
            direction2=[1.0, 0.0, 0.0],
            aboutCM=True,
        )
        problem.addFunction(
            "I_xy",
            functions.MomentOfInertia,
            direction1=[1.0, 0.0, 0.0],
            direction2=[0.0, 1.0, 0.0],
            aboutCM=True,
        )
        problem.addFunction(
            "I_xz",
            functions.MomentOfInertia,
            direction1=[1.0, 0.0, 0.0],
            direction2=[0.0, 0.0, 1.0],
            aboutCM=True,
        )
        problem.addFunction(
            "I_yy",
            functions.MomentOfInertia,
            direction1=[0.0, 1.0, 0.0],
            direction2=[0.0, 1.0, 0.0],
            aboutCM=True,
        )
        problem.addFunction(
            "I_yz",
            functions.MomentOfInertia,
            direction1=[0.0, 1.0, 0.0],
            direction2=[0.0, 0.0, 1.0],
            aboutCM=True,
        )
        problem.addFunction(
            "I_zz",
            functions.MomentOfInertia,
            direction1=[0.0, 0.0, 1.0],
            direction2=[0.0, 0.0, 1.0],
            aboutCM=True,
        )
    return tacs_probs, fea_assembler