docs: def assign_line_results(grid, line_data):
      """
      Write results obtained from PF to graph
  
      Parameters
      -----------
      grid: ding0.network
      line_data: pandas.DataFrame
          DataFrame containing active/reactive at nodes obtained from PF analysis
      """
  
      package_path = ding0.__path__[0]
  
      edges = [edge for edge in grid.graph_edges()
               if (edge['adj_nodes'][0] in grid._graph.nodes() and not isinstance(
              edge['adj_nodes'][0], LVLoadAreaCentreDing0))
               and (
               edge['adj_nodes'][1] in grid._graph.nodes() and not isinstance(
                   edge['adj_nodes'][1], LVLoadAreaCentreDing0))]
  
      decimal_places = 6
      for edge in edges:
          s_res = [
              round(sqrt(
                  max(abs(line_data.loc["MV_{0}_lin_{1}".format(grid.id_db, edge[
                      'branch'].id_db), 'p0'][0]),
                      abs(line_data.loc["MV_{0}_lin_{1}".format(grid.id_db, edge[
                          'branch'].id_db), 'p1'][0])) ** 2 +
                  max(abs(line_data.loc["MV_{0}_lin_{1}".format(grid.id_db, edge[
                      'branch'].id_db), 'q0'][0]),
                      abs(line_data.loc["MV_{0}_lin_{1}".format(grid.id_db, edge[
                          'branch'].id_db), 'q1'][0])) ** 2),decimal_places),
              round(sqrt(
                  max(abs(line_data.loc["MV_{0}_lin_{1}".format(grid.id_db, edge[
                      'branch'].id_db), 'p0'][1]),
                      abs(line_data.loc["MV_{0}_lin_{1}".format(grid.id_db, edge[
                          'branch'].id_db), 'p1'][1])) ** 2 +
                  max(abs(line_data.loc["MV_{0}_lin_{1}".format(grid.id_db, edge[
                      'branch'].id_db), 'q0'][1]),
                      abs(line_data.loc["MV_{0}_lin_{1}".format(grid.id_db, edge[
                          'branch'].id_db), 'q1'][1])) ** 2),decimal_places)]
  
          edge['branch'].s_res = s_res