python - how to compare string in one line with string in next line? -

-

i have file 16,000 lines in it. of them have same format. here simple example if it:

atom    139  c1  dppc   18      17.250  58.420  10.850  1.00  0.00  <...>  atom    189  c1  dppc   19      23.050  20.800  11.000  1.00  0.00

i need check if lines contains string dppc , has identifier 18 forms 50 line block before identifier switches 19, etc.

so now, have following code:

cnt = 0 open('test_file.pdb') f1:     open('out','a') f2:         lines = f1.readlines()         i, line in enumerate(lines):              if "dppc" in line:                    = line.strip()[22:26]                    if a[i] == [i+1]:                        cnt = cnt + 1                    elif a[i] != a[i+1]:                        cnt = 0

and here stuck. found examples how compare subsequent lines similar approach did not work here. still cannot figure out how compare value of a in line[i] value of a in line[i+1].

try (explanations in comments).

data = """atom    139  c1  dppc   18      17.250  58.420  10.850  1.00  0.00 atom    139  c1  dppc   18      17.250  58.420  10.850  1.00  0.00 atom    139  c1  dppc   18      17.250  58.420  10.850  1.00  0.00 atom    139  c1  dppc   18      17.250  58.420  10.850  1.00  0.00 atom    139  c1  dppc   18      17.250  58.420  10.850  1.00  0.00 atom    139  c1  dppc   18      17.250  58.420  10.850  1.00  0.00 atom    139  c1  dppc   18      17.250  58.420  10.850  1.00  0.00 atom    139  c1  dppc   18      17.250  58.420  10.850  1.00  0.00 atom    139  c1  dppc   18      17.250  58.420  10.850  1.00  0.00 atom    189  c1  dppc   19      23.050  20.800  11.000  1.00  0.00 atom    189  c1  dppc   19      23.050  20.800  11.000  1.00  0.00 atom    189  c1  dppc   19      23.050  20.800  11.000  1.00  0.00 atom    189  c1  dppc   19      23.050  20.800  11.000  1.00  0.00 atom    189  c1  dppc   19      23.050  20.800  11.000  1.00  0.00 atom    189  c1  dppc   19      23.050  20.800  11.000  1.00  0.00 atom    189  c1  dppc   19      23.050  20.800  11.000  1.00  0.00 atom    189  c1  dppc   19      23.050  20.800  11.000  1.00  0.00 atom    189  c1  dppc   19      23.050  20.800  11.000  1.00  0.00 atom    189  c1  dppc   19      23.050  20.800  11.000  1.00  0.00"""  # last code seen in 5th column. code = none  # count of lines of current code. count = 0  line in data.split("\n"):     # 5th column.     c = line.split()[4]      # code in 5th column changed.     if c != code:         # if aren't @ start of file, print count         # code ended.         if code:             print("{}: {}".format(code, count))          # rember new code.         code = c      # count line     count = count + 1  # print count last code. print("{}: {}".format(code, count))

output:

18: 9 19: 19

Comments

Post a Comment

Popular posts from this blog

sequelize.js - Sequelize group by with association includes id -

-
so when requesting group sequelize follows: return models.workingcalendar .findall({ attributes: [ 'workingcalendar.periodid', 'workingcalendar.date', 'period.name' ], include: [ { model: models.period, attributes: [] } ], where: { getsudoid: currentgetsudo.id, unitplantid: unitplantid }, group: ['workingcalendar.periodid', 'workingcalendar.date', 'period.name'], }); sequelize run query: select [workingcalendar].[id], [workingcalendar].[periodid], [workingcalendar].[date], [period].[name] [workingcalendars] [workingcalendar] left outer join [periods] [period] on [wor...
Read more

c++ - Migration from QScriptEngine to QJSEngine -

-
i'm migrating qscriptengine code on qjsengine. now, have: class pars { public: static qscriptvalue printmainlog(qscriptcontext* c, qscriptengine* e); }; qscriptvalue pars::printmainlog(qscriptcontext* c, qscriptengine* e) { //some actions return e->globalobject().property(""); } ... qscriptengine engine; ... engine.globalobject().setproperty("printlog",engine.newfunction(pars::printmainlog)); so, user can put printlog("what ever"); in application in , example, qlineedit , function pars::printmainlog evalute. is there way qjsengine? so, user put same printlog("what ever");? way find here , user should put logger.printlog("what ever"); logger class inherited qobject printlog slot.
Read more

delphi - Take screenshot in webcam using VFrames in Console Application -

-
i'm doing program take picture of webcam using delphi xe2 , vframes achieve this, problem have figured out, in graphic application works fine, when use unit in console application, returns me error saying first chance exception @ $76b6b727. exception class eaccessviolation message 'access violation @ address 004a271b in module 'console.exe'. read of address 00000260'. process console.exe (3676) my unit : unit webcam; interface uses sysutils, windows, vcl.imaging.jpeg, vcl.graphics, vsample, vframes, classes; type twebcam = class private procedure newvideoframeevent(sender: tobject; width, height: integer; dataptr: pointer); public constructor create; destructor destroy; override; procedure capture_webcam(take_name: string); end; var web_image: tvideoimage; name_screen: string; implementation constructor twebcam.create; begin inherited create; end; destructor twebcam.destroy; begin inherited destroy; end; ...
Read more